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INTRODUCTION 


The  impetus  for  study  of  close-in  phenomena  of  buried  explosions  stems  from  a  desire  to 
understand  how  the  explosive  disturbance  generates  Sblsmic  waves  at  distances  far  from  underground 
explosions.  It  is  certain  that  the  elastic  waves  produced  are  characterized  by  features  of  the  initial 
explosion  and  by  properties  of  the  propagating  medium  as  well.  Oelineatimi  of  the  characteristics  of 
the  close-in  disturbance  provides  insight  into  the  manner  of  seismic  wave  generation  by  explosions,  a 
central  aim  of  the  Vela- Uniform  Program. 

The  close-in  region  about  a  buried  explosion  may  be  described  as  that  in  which  wave  propagation 
is  nonelastic.  This  close-in  region  may  be  subdivided  into  at  least  two  distinct  subregions  which  are 
characterized  by  the  manner  of  the  medium's  response  to  the  explosive  disturbance.  In  ihe  region 
immediately  adjacent  to  the  explosion  where  disturbance  of  the  medium  is  Intense,  wave  propagation 
is  by  means  of  shock  waves;  this  region  is  termed  hydrodynamic  or  fluid.  At  distances  farther 
removed  from  the  explosion  center,  the  shock  wave  disturbance  decays  to  such  an  extent  that  the 
medium  no  longer  responds  as  a  fluid  but  not  to  the  extent  where  the  medium’s  behavior  is  yet  elastic, 
This  nonfluid-nonelaatic  region  Is  characterized  by  many  possible  complex  phenomena  whose  details 
are  largely  unknown.  Processea  such  as  yielding,  plastic  How,  crushing  and  cracking,  phase  changes, 
anomalous  compressibility,  and  viscosity  are  likely  to  be  encountered  in  this  region.  In  spite  of  the 
vast  number  of  existing  experimental  and  theoretical  publications  oti  these  processes,  it  is  important 
to  realize  that  little  is  actually  known  about  these  (dtenomena  particularly  as  they  may  apply  to  under¬ 
ground  explosions.  However,  numerous  and  valuable  data  from  meaaurements  in  the  nonfluid- 
nonelaatic  region  about  buried  explosions  await  detailed  theoretical  description,  and  from  these,  it  is 
hoped,  we  can  glean  information. 

Significant  progress  has  been  made  during  the  past  two  or  three  yearn  in  theoretical  calculations 
of  the  wave  disturbance  from  buried  explosions  in  the  close-in  region.  Most  su  lking  is  the  success 
of  these  calculations  for  the  dtsturfoonce  in  the  hydroi(ynamic  region.  It  seems  eertsln  that  the  general 
features  of  shock  wave  propagation  in  the  hydrodynamic  region  are  correctly  deacribed  by  the  theory. 
However,  quantitative  agreement  between  tbeoretteal  predtctiona  and  data  recently  obtained  from 
meaaurements  in  the  fluid  region  is  not  In  every  case  obtained,  indicallng  that  refinements  in  theory 
are  required. 


Understandably,  much  less  progress  has  beat  made  in  detailed  calculations  of  the  wave  disturb¬ 
ance  in  the  nonfluiu-nonelastic  region.  This  lack  of  success  is  principal];,  a  result  of  our  ignorance  of 
the  medium's  properties  in  this  regicm.  For  this  regicm  a  model  of  the  medium  is  required  which 
describes  its  prq>erties  and  (iynamic  behavior  under  the  conditions  created  by  buried  explosions. 
Criteria  for  dynamic  yielding  must  be  established,  the  Influence  of  phase  changes  evaluated,  the  signi¬ 
ficance  of  viscosity  and  strain  rate  effects  determined,  and  other  possibly  significant  phenomena 
investigated  before  an  adequate  model  can  be  obtained. 

Once  a  realistic  model  is  obtained,  then  the  explosive  disturbance  may  be  theoretically  described 
throughout  both  the  fluid  and  nonfluid-nonelastic  regions  to  the  point  where  the  disturbance  is  elastic. 

At  that  region  in  space  about  the  eiqdosion  where  the  disturbance  is  described  and  where  the  medium 
responds  elastically,  the  explosion  source  characteristics  for  seismology  may  be  considered  obtained. 
At  present,  the  greatest  barrier  to  a  fuller  understanding  of  seismic  wave  generation  from  buried 
explosions  is  the  lack  of  knowledge  of  the  medium's  properties  in  the  nonfluld-nonelastic  region. 

Much  of  the  information  required  for  construction  of  an  adequate  model  can  be  obtained  from  laboratory 
experiments.  The  remainder  must  come  from  direct  measurement  during  full-scale  explosion  experi¬ 
ments. 

The  research  described  in  this  report  is  an  attempt  to  obtain  some  information  about  close-in 
phenomena  which  will  be  helpful  to  the  realisation  of  objectives  in  the  Vela-lktlform  Program.  Specific 
aims  of  this  research  have  been  the  development  of  instruments  capable  of  making  measurements  in 
the  nonelaatic  region  about  buried  explosions,  investigation  of  medium  properties  of  geologic  solids 
relevant  to  studies  of  the  close-in  region,  and  develo|wnent  of  a  theory  for  description  of  spherical 
wave  propagatioi«  from  burled  explosions. 

Instrumeida  developed  are  a  peak  ahock  preaaure  gage  for  meaaurement  of  100-  to  lOOO-kllobar 
ahoek  preasurea  in  aolida  and  long-period,  large-amplitude  diaplacement  gagea  for  meaauremeots  in 
the  nonfluld-nonelaailc  region.  All  instruments  developed  have  been  tested  In  field  experiments  with 
nuclear  explcwiona  under  conditions  for  which  they  were  designed  These  tests  hsve  been  auecessftil 
and  the  Inatruments  developed  are  eurrenUy  being  used  to  collect  dsta  from  experiments  with  burled 
nuclear  ex|doaiana. 

Medtinm  propeiliss  investigstsd  have  bean  Hugoniots  and  the  Ibigonlet  elasUe  limits  of  gselogtc 
mateHals.  fhigontois  of  voleanic  Udt  and  two  types  of  porous  soil  have  bsea  measured  iUso  the 
Wlusnct  of  perosl^  and  water  contaut  on  the  Hugoniots  of  porous  msdio  has  boon  sxamiasd  fb^oniot 


elastic  limits  have  be^i  investigated  by  means  of  a  new  technique  eiiq>lqying  quartz  piezoelectric 
crystals  to  measure  directly  the  stress-time  profiles  of  finite  amplitude  waves  in  rock  specimens.  By 
this  technique  some  qualitative  informaticm  has  been  obtained  from  measurements  of  stress  wave  pro¬ 
files  in  granite,  basalt,  and  halite. 

A  phenomenological  descr4>tlon  of  spherical  wave  prt^agation  from  buried  explosions  has  been 
developed  and  compared  with  experimental  data.  The  express  purpose  of  this  effort  has  been  to  obtain 
a  simple  analytical  description  without  resort  to  the  formidable  difficulties  of  a  more  rigorous  mathe¬ 
matical  approach  and  without  detailed  consideration  of  the  numerous  possible  complex  phenomena 
which  occur  in  the  close-in  region.  Ob.vectives  were  to  obtain,  in  terms  of  a  few  parameters  which 
characterize  gross  properties  of  media,  c]q>ressions  for  the  decay  of  wave  pressure  with  distance 
which  could  then  be  compared  with  numerous  available  data.  Calculations  have  been  compared  with 
data  from  both  nuclear  and  chemical  mqtloaions  in  granite,  halite,  volcanic  tuft  and  desert  alluvium 
over  the  complete  range  of  pressures  found  in  all  regions  from  the  hydrotiynamic  to  the  elastic.  The 
results  of  calculation  are  presently  as  accurate  as  those  of  other  methods. 


Section  I 


HYDRODYNAMIC  SHOCK  PRESSURE  GAGE 
A.  J.  Chabai,  R.  C.  Bass,  and  H.  L.  Hawk 

Background 

The  purpose  of  developing  an  instrument  that  can  measure  pressures  of  amplitude  from  0.1  to 
1  megabar  (100  thousand  to  1  million  atmospheres)  is  to  obtain  data  on  the  shock  wave  decay  with 
distance  in  the  hydrodynamic  region  about  buried  explosions  and  to  provide  information  for  the  evalua^ 
tion  of  theoretical  calculations.  Conventional  pressure  gages  can  not  operate  at  the  extreme  pressure 
levels  associated  with  shock  waves  in  solids,  nor  are  their  principles  of  operation  generally  applicable 
for  extension  to  the  high  pressures  in  question.  Indeed,  the  material  of  which  the  gage  is  constructed 
may  be  a  significant  perturbation  on  the  shock  wave  whose  pressure  is  to  be  measured.  At  shock 
pressures  of  0.1  to  1  megabar  all  materials  suffer  a  large  compression  and  so  any  transducer  will 
have  its  physical  form  drastically  altered;  in  short,  it  will  be  destroyed  by  shock  waves.  The  extreme 
pressures  and  high  accelerations  associated  with  shock  waves  present  a  difficult  environment  for  any 
transducer.  Clearly  any  transducer  in  this  environment  must  complete  its  measurement  before  it  and 
its  component  circuitry  and  cabling  are  destroyed. 

Of  the  five  hydrodynamic  quantities  -  pressure,  density,  temperature,  particle  velocity  and  shock 
wave  velocity  -  which  characterize  any  shock  wave,  the  one  most  susceptible  to  direct  measurement  is 
shock  wave  velocity.  Particle  velocity  comes  next  in  ease  of  measurement. Pressure  and  temr'orature 
behind  shock  waves  in  solids  have  never  been  directly  measured.  While  direct  measurements  of 
density  behind  shock  waves  in  solids  have  been  made  by  flash-X-ray  tecliniques  in  laboratory  experi¬ 
ments  the  method  is  not  considered  feasible  for  use  in  a  field  gage. 

The  boundary  conditions  at  a  shock  front  specifying  conservation  of  mass  and  momentum  contain 
the  four  unknown  hydrodynamic  quantities  -  pressure,  density,  particle  velocity,  and  shock  velocity. 

If  any  two  of  these  can  bo  measured,  then  the  remaining  two  quantities  are  determined  by  the  conserva¬ 
tion  equations.  In  laboratory  experiments^  on  shock  waves  in  solids,  shock  velocity  and  particle 
velocity  arc  measured;  shock  pressure  and  density  are  then  inferred  from  the  conservation  equations. 

’'From  measurements  of  free  surface  velocity  particle  velocity  is  Inferred  (see  Reference  1). 
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This  laboratory  technique  has  been  adopted  by  Lombard^  in  development  of  a  peak  shock  pressure  gage 
for  field  use.  Another  laboratory  technique^  used  to  determine  hydrodynamic  quantities  behind  shock 
waves  is  the  impedance-mismatch  method.  In  this  method  use  is  made  of  some  material  whose 
ilugoniot  is  known,  and  measurement  of  only  one  quantity,  shock  velocity,  is  required  in  order  to  infer 
pressure,  density,  and  particle  velocity.  Since  shock  velocity  is  more  easily  measured  than  particle 
velocity  the  principles  of  the  impedance-mismatch  technique  have  been  chosen  for  the  development  of 
the  peak  shock  pressure  transducer  described  in  this  report. 

Other  methods  for  measurement  of  peak  shock  pressures,  employing  new  and  different  principlesi 

3 

are  currently  being  investigated.  For  example,  Lombard  has  made  use  of  electric  current  released 
by  Lucite  under  the  influence  of  shock  pressure.  It  is  found  that  over  a  certain  range  of  shock  pressures 
the  peak  current  released  by  Lucite  is  directly  proportional  to  peak  shock  pressure. 

Research  is  also  being  conducted  to  develop  shock  pressure  gages  capable  of  measuring  pressure 
as  a  function  of  time  behind  the  incident  shock  fronts.  Significant  progress  along  these  lines  has  been 

4 

made.  Keough  and  Bernstein  have  employed  a  technique  where  the  changes  in  resistance  of  a  wire 
over  w'hich  a  shock  wave  is  passing  are  related  to  the  pressure  changes  behind  the  wave.  Manganin 
wire  appears  to  be  the  most  promising  sensor.  The  percentage  increase  in  wire  resistance  is  a  linear 
function  of  shock  front  pressure  over  the  range  investigated  (10  to  150  kilobars).  Measurement  of 
pressure-time  history  behind  a  shock  wave  can  be  made  for  times  of  the  order  of  tens  of  microseconds 
before  the  gage  is  destroyed. 

Another  gage  being  developed  for  measurement  of  shock  wave  pressure-time  profiles  consists 
of  a  disk-shaped  cell  filled  with  distilled  water  into  which  two  silver  electrodes  are  inserted.  Applying 
a  potential  to  the  electrodes  and  passing  a  shock  wave  through  the  cell  results  in  an  output  signal  which 
varies  directly  with  incident  shock  pressure  and  directly  with  the  initially  applied  potential.  No  signal 
is  observed  with  zero  applied  potential.  Pressure-time  waveforms  are  observed  for  periods  of  10  to 
20  psec  and  have  shapes  similar  to  those  obtained  with  the  manganin  resistance  wire  gage. 

While  not  all  of  the  pressure  gages  currently  under  development  have  been  extensively  tested  in 
field  experiments,  their  eventual  use  ’vill  provide  much  information  about  wave  propagation  in  the  non- 
(•lastlc  region  about  buried  explosions. 

As  a  result  of  research  activity  in  shock  pressure  transducers,  roughly  25  measurements  of  peak 
shock  pressure  have  been  successfully  performed  during  the  past  two  years  in  the  hydrodynamic  region 


about  buried  explosions.  The  measurements  were  made  using  the  impedance-mismatch  gages  and  the 
2  7  8  9 

Lombard  gage.  Shock  pressures  recorded  '  ‘  have  been  from  20  to  more  than  600  kUobars  in  granite, 
volcanic  tuff,  and  desert  alluvium.  These  pressure  measurements  have  yielded  valuable  information 
and,  together  with  independent  measurements  of  shock  wave  position  as  a  function  of  time,  provide  the 
basis  for  confidence  in  the  theoretical  description  of  shock  wave  propagation  in  the  hydrodynamic  region 
about  buried  explosions. 


Impedance- Mismatch  Shock  Pressure  Gage 

Conservation  of  momentum  and  mass  at  the  front  of  a  shock  wave  are  specified  by  the  equations^ 

P  =  p  Vu  (1.1 

o 


P^V  =  p(V-u),  (1.2] 

where  P,  p,  and  u  are  respectively  pressure,  density  and  particle  velocity  behind  the  shock  wave, 

V  is  shock  velocity,  and  p^  is  the  density  of  the  medium  into  which  the  wave  progresses.  In  Equations 
1.1  and  1.2  it  is  assumed  that  the  medium  is  at  rest,  u^  =  0,  and  that  ambient  medium  pressure,  P^, 
is  essentially  zero.  The  increase  in  internal  energy,  is  given  by 


The  Hugoniot  of  a  material  is  a  relationship  between  P  and  p  and  represents  the  locus  of  the 
thermodynamic  states  which  ctin  be  attained  in  the  material  by  mesins  of  shock  waves.  Every  material 
has  a  unique  Hugoniot.  From  laboratory  experiments^^  it  has  been  found  that  for  most  solids  shock 
velocity  is  a  linear  function  of  particle  velocity  over  a  wide  range  of  shock  pressures, 

V  =  C  +  Su  .  (1.4] 

Eliminating  velocities  in  Equations  1.1,  1.2  and  1.4,  the  Hugoniot  is  obtained  in  terms  of  the 

material  constants  p  ,  C  and  S; 

o 


p  C  p(p-p  ) 
o _ *^0 

[p^S.p(S-l)j^ 


13 


By  eliminating  shock  velocity  from  Equations  1.1  and  1.4  another  form  for  the  Hugoniot  can  be 
expressed  in  terms  of  pressure,  particle  vdlocity.  and  the  material  constants: 

P  =  p  (C+Su)u  .  (1.6) 

o 

The  principle  of  operation  of  the  impedance-mismatch  gage  is  an  extension  of  well  known^ 
laboratory  techniques  and  is  illustrated  by  the  following  example.  Consider  a  plane  shock  wave 
progressing  through  a  medium  adtose  Hugoniot  is  unknown,  and  assume  that  this  incident  wave 
encounters  two  or  more  different  materials  of  known  Hugoniot  in  its  path  (see  Figure  1. 1). 


p. 

- >  V 

Material 

1 

2 

Medium 

'’l 

Material 

P  a  0 

0 

Material 

“l 

3 

u  =  0 

0 

Figure  1,1  Shock  Wave  Incident  on  Dissimilar  Materials 

At  the  interfaces  of  the  materials  with  known  Hugoniots,  shock  waves  are  transmitted  into  these 
materials  and  waves  are  reflected  back  into  the  medium.  By  measuring  the  shock  wave  transit  times 
through  materials  2  and  3,  shock  pressure  in  each  material  is  determined  since  their  Hugoniots  are 
known.  This  information  together  with  a  measurement  of  the  incident  shock  velocity  allows  the 
incident  shock  pressure  to  be  inferred. 

The  known  Hugoniots  (Equation  1.6)  of  materials  2  and  3  of  Figure  1.1  are  illustrated  in  Figure  1.2 
together  with  the  medium  Hugoniot,  1,  considered  unknown.  From  the  shock-wave  transit-time 
measurements  in  materials  1,  2,  and  3  we  obtain  shock  velocities,  and  with  the  known  values  of 
material  densities,  shock  impedance,  p^V,  is  determined  for  each  material. 

By  Equation  1.1  it  is  seen  that  shock  impedance  is  represented  by  a  line  of  slope  p^V  passing 
through  the  origin  of  coordinates  in  Figure  1.2,  Since  the  momentum  equation  must  be  satisfied  and 
since  the  pressure  of  the  shock  wave  must  be  given  by  one  of  the  points  on  the  Hugoniot  curve,  the  inter¬ 
section  of  the  line  of  slope  p^V  with  the  medium  Hugoniot  curve  determines  the  pressure  and  particle 
velocity  of  the  shock  wave  whose  Impedance  is  p^V.  This  point  of  intersection  is  given  by  the  solution 


Equations  1.1  and  1.6  once  p^V  has  been  measured.  It  is  because  the  Hug(mi<4s  of  materials  2  mid  3 
are  known  that  pressures  and  and  corre^iondinc  particle  velocities  u^  and  u^  may  be  obtained 
from  shock  transit  time  measurements  in  these  materials. 


p 


Figure  1.2  Illustration  of  Peak  Shock  Pressure  Determination  From 
Pres  sure- Particle  Velocity  Diagram 

Now  the  transmitted  pi  assures  and  particle  velocities  established  in  materials  2  and  3  are 
determined  by  the  strength  of  the  incident  shock  wave  in  the  medium  material,  1.  Since  the  shock 
impedances  of  materials  2  and  3  differ  from  the  medium  impedance,  the  transmitted  shock  pressures 
will  be  different  from  those  of  the  incident  shock  and  waves  will  be  reflected  back  into  the  medium 
material  at  the  interfaces.  Applying  the  conservation  equations  to  the  system  of  waves  resulting  after 
the  incident  shock  passes  an  interface,  it  is  found  that  pressure  and  particle  velocities  of  a  reflected 
wave  must  lie  on  a  curve  called  the  "reflection  Hugoniot."^’^^'^*  This  curve.  Illustrated  by  Ir  for  the 
medium  material  in  Figure  1.2,  represents  the  reflected  wave  traveling  in  a  direction  opposite  from 
the  incident  wave  and  determines  the  locus  of  states.  (P,  u),  which  may  be  attained  by  reflected  waves 
starting  from  the  shocked  state.  (P^,  u^).  From  the  boundary  conditions  that  pressure  and  particle 
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velocity  must  be  ccmtinuous  across  an  interface,  It-is.  seen  that  the  intersection  of  curves  Ir  and  2  of 
figure  1.2  determines  the  values  of  pressure  and  particle  velocity  for  both  the  transmitted  and  reflected 
waves  at  the  interface  between  materials  1  and  2.  hi  this  case  it  is  seen  that  the  reflected  wave  is  a 
shock  wave  (P^  >  P^)  since  the  shock  impedance  of  material  2  is  greater  than  that  of  the  medium 
material.  Similarly,  considering  materials  1  and  3,  it  is  seen  from  Figure  1.2  that  from  this  interface 

a  rarefaction  wave  (P„  <  P.)  is  reflected  since  (p  V).  <  (p  V),  . 

3  1  o  o  o  1 

It  is  found^  experimentally  that  within  a  few  percent  the  reflection  Hugoniot,  Ir,  is  the  mirror 
image  of  the  Hugoniot  curve,  1,  about  the  line  u  =  u^  .  Consequently  the  reflection  Hugoniot  may  be 
expressed  as 

P  =  p^  [s(2uj-u)^  +  C{2uj-u)j  .  (1.7) 

The  points  (P^,  u^)  and  (P^,  u^),  establijhed  by  shock  transit  time  measurements  in  materials  2 
and  3,  are  also  points  on  the  unknown  Hugoniot  of  the  medium  (noted  by  open  circles  in  Figure  1.2), 

The  impedance-mismatch  gage  can  make  use  of  this  fact  in  obtaining  information  on  the  Hugoniot  of  the 
medium  material. 

If  materials  2  and  3  are  chosen  to  have  a  higher  and  lower  shock  impedance  than  the  medium  as 
shown  in  Figure  1.2,  then  the  incident  shock  pressure  may  be  estimated  by  Interpolation  as  that 
pressure,  determined  by  the  intersection  of  the  line  of  slope  (PqV)^  and  a  line  joining  the  points 

(Pj,  u^)  and  (Pg,  Ug)  on  the  Hugoniot  of  materials  2  and  3.  If  the  two  materials  of  known  Hugoniot  both 
have  higher  or  lower  shock  Impedances  than  the  medium,  then  an  extrapolation  procedure  may  be 
employed  to  establish  P^^  .  The  interpolated  pressure,  Pjj^j.  will,  as  seen  in  Figure  1,2,  always  be 
greater  than  the  incident  shock  pressure,  P^;  however,  it  is  found  empirically  from  examination  of  the 
interpolatim  procedure  for  various  combinations  of  materials  1,  2,  and  3  and  for  various  incident 
pressures,  P^,  that  the  difference  between  Pj^^  and  P^  resulting  from  interpolation  will,  In  almost 
every  case,  be  much  less  than  the  maximum  combined  percentage  error  (25  percent)  expected  from 
other  sources  in  the  measurement  process.  A  10  percent  error  appt^ars  to  be  a  more  realistic  error 
incurred  in  P^  by  a  straight  line  interpolation  procedure. 

Once  Pj^  and  also  u^^  are  obtained  Equation  1.7  and  the  points  (P^.  u^)  and  (P^,  u^)  may  be 
utilised  to  estimate  the  unknown  Hugoniot  of  the  medium.  Since  medium  density  is  known  and  u^^ 
has  been  obtained,  the  two  unknown  constants,  and  C^,  which  determine  the  medium  Hugoniot 
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(Equation  1.5),  may  be  evaluated.  This  procedure  gives  at  least  an  estimate  of  the  Hugoniot  in  the 
vicinity  of  the  pressure,  P^.  From  similar  measurements  of  shock  pressure  at  different  levels  of 
the  incident  pressure,  (at  different  distances  from  an  explosion),  additional  Hugoniot  points  and 
values  of  and  near  the  pressure  P^  are  obtained,  which,  when  all  taken  together,  provide  an 
estimate  of  the  in  situ  medium  Hugoniot  curve  (curve  1  of  Figure  1.2)  over  the  range  of  measured 
pressures. 

It  must  be  emphasized  that  the  impedance- mismatch  method  of  measuring  shock  pressures  is 
strictly  applicable  only  to  strong  shock  waves.  At  lower  wave  pressures  where,  for  example,  phase 
changes  may  be  induced  in  the  medium  material,  the  Hugoniot  (curve  1  in  Figure  1.2)  will  have  a  cusp 
or  point  of  inflection,  and  the  reflection  Hugoniot  (curve  Ir)  will  likely  not  be  the  mirror  image  of  the 
Hugoniot,  particularly  that  portion  of  the  reflection  Hugoniot  which  joins  the  initial  shock  point  (Pj,  u^) 
with  a  lower  impedance  point  such  as  (Pg.  Ug)  in  Figure  1.2.  Fortunately,  it  appears  from  presently 
available  data  that  even  marked  phase  changes  induced  by  shock  waves  in  geologic  materials  do  not 
produce  marked  cusps  or  deviations  from  a  smooth  (P-u)  Hugoniot  curve  (Equation  1.6)  whose  C  and  S 
values  are  evaluated  by  strong  shock  measurements.  In  addition,  that  portion  of  the  reflection 
Hugoniot  joining  points  such  as  (P^,  Uj)  and  (Pg,  Ug)  in  Figure  1.2  is  actually  an  adiabat  which  is 
absolutely  unknown  for  geologic  materials  and  which,  more  than  likely,  has  no  Inflection  points  as  does 
a  Hugoniot  with  phase  changes.  The  observation  that  a  Hugoniot  curve  in  the  P-u  plane  is  not  grossly 
altered  from  a  continuous  curve  as  given  by  Equation  1.6  by  phenomena  such  as  phase  changes,  and  the 
qualitative  belief  that  the  adiabatic  portion  of  the  reflection  Hugoniot  will  be  a  continuous  curve  and 
nearly  equal  to  the  mirror  image  of  the  Hugoniot,  tend  to  be  compensations  for  errors  incurred  in 
Pjjj^  when  an  impedance-mismatch  gage  is  used  at  pressure  levels  where  phase  changes  may  occur. 

Until  more  precise  methods  are  available  for  measuring  shock  pressures,  it  is  felt  that  the  impedance- 
mismatch  technique  will  provide  reasonable  values  of  incident  shock  pressure  even  at  pressures  where 
phenomena  such  as  phase  changes  cause  cusps  or  departures  from  a  smooth  Hugoniot  curve. 

To  test  the  impedance- mismatch  method  of  determining  incident  shock  pressures,  laboratory 
experiments  were  conducted  in  which  2024  aluminum  (whose  Hugoniot  is  known^)  was  used  to  represent 
the  medium.  Brass  and  Plexiglas  were  used  as  the  mismatch  materials  2  uid  3  of  known 
Hugoniot.  Plane  shock  waves  were  generated  in  the  aluminum  using  TNT  and  Composition- B  explosives. 
The  210-kllobBr  shock  from  TNT  and  the  350-kilobar  shock  from  Composltion-B  expected  in  aluminum 
were  measured  by  the  impedance-mismatch  method  with  errors  of  about  5  percent. 
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As  envisaged  for  field  use  in  drill  holes  radial  to  an  explosion  center  the  impedance- mismatch 
gage  will  hive  a  large  circular  disk  (Figure  1.3)  cut  from  rock  cores  taken  at  or  near  positions  in  the 
medium  where  pressure  measurements  are  desired. 


Figure  1.3  Schematic  of  Impedance- Mismatch  Gage  Assembly 

A  typical  arrangement  of  the  basic  elements  for  an  impedance-mismatch  gage  is  shown  in  Figure 
1.3.  On  the  flat  surface  of  the  large  circular  disk  of  medium  material  (8  inches  in  diameter  and  1  to  2 
inches  in  thickness)  are  mounted  three  smaller  disks  (dimensions  1-1/2  inch  diameter  and  1  cm  thick). 
The  two  disks  of  mismatch  material  and  the  small  disk  of  medium  material  have  their  centers  on  a 
circle  of  1-1/2  inch  radius  and  are  spaced  120  degrees  apart.  Dlameter-to-thlcknuss  ratio  of  the  small 
disks  must  be  about  1  or  more  in  order  to  avoid  the  influence  of  edge  effects^^  on  detectors  as  the  shock 
wave  progresses  through  tin  disks.  Shock  wave  detectors  are  placed  at  the  center  of  the  back  face  of 
each  of  the  small  disks  (numbered  1,  2  and  3  in  Figure  1.3).  In  addition,  three  shock  detectors 
(numbered  4,  S  and  6)  are  placed  on  the  back  surface  of  the  large  diak.  These  latter  three,  mounted  on 
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a  circle  of  3/4>inch  radius,  are  spaced  120  degrees  a{>art  and  located  symmetrically  to  the  small  disks 
as  shown  in  the  figure. 

Detectors  2  and  3  each  in  combination  with  any  one  or  all  of  the  dt  ^ectors  4,  5  and  6  serve  to 
determine  shock  transit  times  in  the  mismatch  materials.  Detector  1  with  4.  5  or  6  establishes  transit 
time  in  the  medium  material.  Thus,  for  example,  in  the  case  of  a  normally  incidait  plane  shock  wave, 
the  impedance  slopes,  (p  V),,  (p  VL,  and  (p  V)„  required  in  Figure  1.2  for  determination  of  medium 
shock  pressure,  may  be  obtained  from  transit  time  measurement  of  detector  combinations  1  and  4, 

2  and  5,  and  3  and  6,  respectively. 

When  the  incident  wave  is  not  normal,  correction  for  tilt  (angle  between  the  longitudinal  axis  of 
the  large  medium  disk  and  a  normal  to  the  plane  of  the  incident  shock  wave)  must  be  made.  It  is  for  the 
purpose  of  measuring  shock  tilt  that  three  detectors,  4,  5,  and  6,  are  used  on  the  large  medium  disk. 
Detectors  4,  5,  6.  and  1  form  a  set  of  four  noncoplanar  points  which  in  principle  are  sufficient  to 
determine  both  degree  of  shock  tilt  and  shock  velocity  in  the  medium  material. 

Consider,  for  example,  a  plane  wave  incident  on  the  shock  gage  and  tilted  with  respect  to  the 
longitudinal  axis  of  the  gage  of  Figure  1.3.  Suppose  that  detector  number  4  first  senses  the  shock  wave. 
Construct  a  cartesian  system  of  coordinates  with  the  center  of  detector  4  as  origin  and  with  one  of  the 
coordinate  axes  parallel  to  the  longitudinal  gage  axis.  Then  from  the  four  equations 

VAtj^  ■  AXj^costfj  +  Ztyj^costtj  +  A*j^coBOg 

'"^‘64  '  ^*64"°"“l  ^y64‘=°“"2  *  ^'64‘=°“3  «.8) 

2  2  2 
1  ■  cos  e  j  +  cos  o  2  *  “  3  • 

incident  shock  velocity.  V,  and  its  direction  cosines  coso^  coso^,  and  coso^  may  be  determined  in 
terms  of  the  known  positions  of  detectors  1.  5  and  6  with  respect  to  detector  4  and  in  terms  of  Ino 
measured  time  differences  ^^64*  ^  solution  of  U.e  equations  is  poattible  only  11  the 

four  detectors  are  not  coplanar.  B  is  assumed  that  a  field  ga;{e  can  be  located  with  some  degree  of 
accuracy  so  that  shock  wave  tUt  will  not  be  too  groat  and  so  that,  possihle  shock  wavs  refraction  effects 
will  bo  negligible. 
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For  experiments  in  which  drill  holes  radial  to  an  explosion  ceiter  are  not  feasible  but  in  which 
vertical  drill  holes  to  the  depth  of  an  explosion  center  are  available,  the  impedance- mismatch  gage  can 
be  a  rectangular  slab  with  elem^ts  mounted  on  the  back  face  in  a  manner  similar  to  that  shown  in 
Figure  1.3.  Whenever  a  gage  is  placed  in  a  vertical  drill  hole,  orientation  of  the  gage  will  be  required. 
Field  experience  with  orientation  has  shown  that,  by  means  of  gyroscopes  attached  to  gages,  an  axis 
normal  to  the  gage  or  slab  face  may  be  aligned  in  a  1000-foot-deep  hole  to  within  1  degree  with  a 
horizontal  line  passing  through  the  explosion  center. 

hi  final  form  for  field  installation  the  gage  will  be  mounted  Ln  a  canister  which  protects  the  working 
elements  of  the  assembly  and  to  which  the  signal  cables  and  gage  lowering  equipment  may  be  attached. 
The  front  surface  of  the  medium  disk  must,  of  course,  be  completely  exposed  so  as  not  to  impose  a 
perturbatimi  on  the  pressure  measurement.  Bonding  the  medium  disk  to  the  Jg  situ  material  is  discussed 
below. 

Shock  Wave  Detectors 

Several  types  of  shock  detector  elements  have  been  considered  for  use  with  the  impedance- 

mismatch  gage.  The  properties  of  both  piezoelectric  ceramic  detectors  and  shorting  pin  detectors  have 
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been  examined.  In  addition,  a  shorting  'circular  switch*  type  of  detector  has  been  investigated.  Up 

to  shock  pressures  of  about  1  megabar,  piezoelectric  ceramics  of  lead-zirconate-titanate  "re  found  to 

•8 

be  by  far  the  best:  however,  any  detector  element  of  fast  rise  time  (-*-10  sec)  r<nd  output  of  10  or 
more  volts  would  serve  as  well. 

to  Figure  1.4  are  shown  the  componmts  of  a  PZT  oetectf^r  which  has  proven  to  be  extremely 
reliable.  A  small  wafer*  of  PZT,  0.120  inch  in  diameter  bv  0,020  inch  in  thickness  is  sandwiched  into 
a  small  brass  housing  onto  which  a  standard  Microdot  coaxial  cable  connector  is  attached.  For  use 
with  an  impedance  gage  the  assembly  is  simply  cemented  into  position,  usually  with  an  epoxy  mixture. 

It  has  been  found  that  the  PZT  detector  of  Figure  1.4  releases  considerable  amounts  of  charge 
over  a  relatively  wide  range  of  pressiu'ea.  For  example,  at  a  shock  pressure  of  300  kilobars  about  100 
volts  are  generated  by  the  PZT  wafer  across  a  30-ohm  resistive  load.  As  shock  pressure  is  Increased, 
however,  the  ou^mt  signal  progressively  decreases.  Precise  experiments  by  Reynolds  and  Seay^^ 
have  shown  a  marked  decrease  in  charge  per  unit  area  released  by  PZT  as  shock  pressure  increases 
from  20  to  ISO  kilobrnrs.  Our  InvsMigations  of  the  PZT  detector  shown  in  Figure  1.4  indicate  that 

*Obi«iasd  from  Oevite  Corp. .  Zlsetronic  Cwnponents  Div. ,  3405  Perkins  Ave. .  Cleveland  I.  Ohio. 
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output  voltages  from  PZT  are  still  sufficiently  high  and  usable  in  a  pressure  gage  up  to  1  -megabar 
pressures.  At  1-megabar  pressure,  the  output  is  about  10  volts  across  a  50-ohm  load.  At  slwck 
pressures  in  excess  of  1  megabar  PZT  ceramics  are  not  adequate  as  dcstectors  for  field  gages.  For 
shock  wave  detection  at  pressures  in  excess  of  1  megabar,  we  have  used  the  0.032-inch  coaxial  self¬ 
shorting  pins  manufactured  by  Edgerton,  Germeshausen  and  Grier,  b)c. ;  however,  the  shorting  pins  do 
not  have  the  same  reliability  that  PZT  ceramics  have  at  lower  pressures,  bivestigations  of  different 
types  of  piezoelectric  ceramics  other  than  PZT  are  being  made  in  the  hope  of  finding  one  which  will 
operate  at  pressures  above  1  megabar. 


CUTAWAY  VIEW  OF  STANDARD  BRASS  INNER  PZT  BRASS  OUTER 

MICRODOT  COAXIAL  ELECTRODE  CERAMIC  ELECTRODE 

CABLE  CONNECTOR  (CLEVTIE  PZT-4) 


SCALE  -  lO/l  APPROX. 


Figure  1.4  Assembly  Drawing  of  PZT  Shock  Wave  Detector 


Delay  Time  Coding  Mixer  Circuit 

When  a  shock  wave  is  ni^arly  normally  incident  on  the  gage  of  Figure  1.3.  output  signals  from 
detectors  4,  5,  and  e  will  be  nearly  coincident,  Also,  depending  on  tl»e  incident  shock  pressure  and  the 
material  of  the  medium  and  mismatch  disks,  signals  from  detectors  I,  2,  and  3  may  occur  nearly 
aimullaneously.  In  addition,  since  the  degree  of  tilt  of  the  shock  wave  at  the  pressure  gage  cannot  he 
known  before  a  measurement  is  msde.  one  doee  iwt  know,  for  example,  which  of  detectors  4,  S.  or  6 
will  give  sn  output  signal  first  Td  Identify  each  detector  signal  and  to  dlidingutsh  unambiguously  those 
signals  which  occur  nearly  elmnitaneouely,  a  Aday  time  coding  mixer  elreuU  has  been  developed. 
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An  electronic  circuit  (Figure  l.S),  developed  for  use  with  an  impedance- mismatch  pressure  gage, 
is  employed  to  identify  detector  output  signals  and  to  allow  the  times  at  which  a  shock  wave  activates  the 
detectors  to  be  measured  with  some  precision. 


Figure  1.5.  Circuit  Flow  Diagram  ter  Time-Delay  Coding  MUer  Circuitry 

The  time  delay  coding  mixer  fTIX^M)  circuit  performs  the  function  of  delaying  in  time  several  of 
the  output  signals  from  detectors  and  mixing  these  signals  in  appropriate  conthinations  for  recording. 
Delaying  signals  by  prescribed  amounts  serves  the  purpose  of  identifying  signals  with  particular 
detectors  and  of  increasing  the  separation  in  lime  of  those  signals  which  occur  nearly  simultaneously 
(two  output  signals  are  considered  simultaneous  when  they  are  seftarated  In  time  by  an  amount  which  is 
less  than  the  pulse  width  of  either  signal).  Provisions  are  made  in  the  TDCM  for  receiving  pulses  from 
a  pulse  generator  for  presluH  calibration  of  time  delays  and  for  receiving  a  sero  time  signal. 

From  Figure  1.5  it  is  seen  that  output  aignais  from  detectors  4.  5,  f,  and  )  -  which  measure  shock 
velocity  and  lUl  of  the  shock  wave  In  the  medium  •  are  delayed  and  mixed  in  the  TDCM  circuit  and  fed  out 
on  a  single  line,  ehannet  A.  to  the  recording  station.  IHgnals  from  detectors  t,  4.  and  k  •  which 
measure  shock  velocity  in  one  mlsmatrh  material  •  are  delayed,  mixed,  and  fed  out  on  channel  R. 
SimUarly,  aivwin  from  detectors  3,  4,  and  S  on  chnnnol  C  dMermine  ahech  veioelly  in  the  other  miamalch 


material.  There  is  some  redundancy  of  signals  on  channels  B  and  C  used  to  obtain  velocity  in  the  mis¬ 
match  material.  This  redundancy  has  been  specifically  designed  into  the  circuit  as  a  safeguard  against 
the  possible  failure  of  one  or  two  of  the  PZT  detectors,  4,  5,  and  6. 

I'  igure  1. 6  is  a  block  diagram  of  the  TDCM.  Its  principal  circuit  for  reception  of  signals  from  the 
detectors  inyolves  for  each  detector  an  "and"  gate,  a  26  psec  blanking  monostable  multivibrator  (MMV), 
a  delay  MMV  and  inverter  (for  detectors  4,  5,  and  6  only),  a  pulse-shaping  MMV,  an  "or"  gate,  and  an 
emitter  follower.  The  subsidiary  and  optional  enabling  circuit  consists  of  a  120  psec  delay  MMV,  an 
inverter,  and  an  enabling  MMV  for  operation  of  the  "and"  gates. 


Figure  l.C  Block  Diagram  of  Time-Delay  Coding  Mixer  Circuit 

The  optional  enabling  circuit  has  been  designed  for  use  with  pressure  gages  that  will  be  placed 
very  close  to  explosions.  Its  purpose  is  to  blank  out  any  early  extraneous  signals  from  an  explosion 
and  prevent  them  from  triggering  the  TDCM  and  recording  scopes.  When  used,  the  enabling  circuit, 
as  shown,  provides  a  500-psec  -long  signal  which  activates  the  "and"  gates  and  makes  them  receptive 
to  detector  signals  from  120  psec  to  620  psec  after  zero  time.  The  time  at  which  the  "and"  gates  are 
activated  and  the  duration  of  activation  can  be  varied  to  suit  the  requirement  of  a  particular  pressure 
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gage.  At  shock  pressure  gage  positions  where  it  is  felt  that  extraneous  signals  will  not  be  troublesome, 
the  enabling  circuit  is  not  used,  the  "and"  gates  are  bypassed,  and  detector  signals  are  presented 
directly  to  the  26  psec  delay  MMV's, 

In  the  main  TDCM  circuit,  the  26  psec  blanking  MMV's  are  used  to  eliminate  noise  and  undesired 
signals,  such  as  cable  breaks  which  may  occur  after  the  measurement  signals  are  obtained  from 
detectors. 

The  delay  MMV's  take  the  nearly  simultaneous  signals  from  detectors  4,  5,  and  6  and  delay  each 
signal  by  a  different  amount  so  that  it  may  be  identified  with  the  detector  which  generated  it.  The  p'artic 
ular  values  of  delay  times  shown  in  Figure  1.4  were  more  or  less  arbitrarily  chosen  and  can  be  varied 
according  to  estimated  shock  wave  transit  tiines  in  the  gage  and  to  meet  conditions  required  for  record¬ 
ing. 

From  inverters  or  directly  from  blanking  MMV's,  signals  are  passed  to  the  pulse-shaping  MMV's 
after  first  being  differentiated.  Output  from  a  pulse-shaping  MMV  is  a  single  positive  pulse  of  0.2  psec 
duration  and  of  about  8-volt  amplitude  (see  Figure  1.7),  , 


Figure  1.7  Output  of  Pulse  Shaping  MMV,  2  Volts/cm,  0.1  psec/cm 

Figure  1.8  shows  the  tyjjical  operation  performed  by  the  TDCM  on  the  signal  from  detector 
number  6. 

Output  signals  from  detectors  1,  4,  5,  and  6  -  after  being  delayed,  inverted,  differentiated,  and 
shaped  -  are  mixed  by  the  "or"  gate  and  then  passed  to  an  emitter  follower  whose  output  is  channel  A. 

In  Figure  1.9  are  shown  output  displays  which  arc  expected  from  channel  A  during  calibration  and  for  a 
normal  incidence  shock  wave  with  a  velocity  of  2.5  mm/psec  in  a  medium  disk  with  a  thickness  of  1  cm. 


Figure  1.8  Typical  Operation  of  TDCM,  Signal  Output  From  Detector  No.  6 


A  prototype  TDCM  circuit  has  been  built  and  its  characteristics  and  electronic  operation  evaluated. 
Figures  1.10(a)  through  l.lO(o)  present  detailed  circuit  diagrams  of  the  TDCM.  Figures  1.11(a)  and  1.11(b) 
show  two  views  of  the  prototype  circuit.  Channel  A,  B,  and  C  outputs  resulting  from  calibration  signals 
applied  simultaneously  to  simulate  detector  signals  1  through  G  are  shown  in  Figure  1.12.  Stability  tests 
indicate  that  delay  times  produced  by  tlic  TDCM  are  constant  to  better  than  0,02  /.'sec  over  periods  of  more 
than  30  minutes. 
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NOTES 

1  resistance  values  are  in  ohms  unless  otherwise  5  the  followins  transistors  are  snigis  type; 

SPECiFiEtL  O'  THRU  Q6,  QB  THRU  012.  QlS  THRU  022,  023 

026,  027,  OaCL  031,  Q3T  thru  QAL  ALL  OTHERS  ARE 

2  capacitance  values  ask  in  UICROMICROFARaDS  2M7A4  TYPE. 

UNLESS  otherwise  specified. 


Figure  1.10(d)  Enabling  Circuit  Diagram 


calibrate 


NOTES 


1  resistance  values  are  in  ohms  unless  otherwise 
SPECiFiEa 

2  CAPACITANCE  VALUES  ARE  IN  MICROMICROFARaOS 
UNLESS  otherwise  SPECIFlEa 


transistors  are  2Niei3  type: 

Q22.  023 

8N*4  0TH«4  are 


Figure  1, 10(e)  Calibration  Circuit  Diagram 


2G 


ll(a>  Front  View  of  Prototype  TDCM  Figure  1.11(b)  Rear  View  of  Prototype  TDCM 


Figure  1.12  Calibration  Tracta  of  Channela  A.  B  and  C. 
1  volt/cm.  2  fisec/cm 
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A  plane  wave  explosive  experiment  was  conducted  to  test  the  TDCM  operation  at  shock  pressure 
levels  simUar  to  those  expected  in  the  field.  Composition-B  explosive  was  used  to  generate  a  450-kilobar 
shock  wave  in  brass.  On  the  back  surface  of  tlie  brass  plate  were  placed  three  shock  detectors  (corre¬ 
sponding  to  4,  5,  and  6  of  Figure  1.3).  In  addition  a  1  cm  thick  disk  of  brass  with  a  shock  detector  on  its 
back  surface  (corresponding  to  1  of  Figure  1.3)  was  mounted  on  the  back  of  the  brass  plate.  The  test 
then  was  an  evaluation  of  the  channel  A  portion  of  the  TDCM.  This  test  was  successful  and  the  TDCM 
circuit  functioned  properly. 

Performance  of  the  TDCM  circuit  appears  adequate  for  field  use.  Modifications  of  the  basic 
circuit  can  be  made  to  accommodate  specific  applications  without  undue  difficulties. 

With  the  TDCM  circuit  the  impedance-mismatch  gage  can  operate  as  described,  yielding  shock 
velocity  in  the  medium,  degree  of  tilt  of  the  shock  wave,  shock  pressure  in  the  medium,  and  information 
on  the  medium  Hugoniot  in  the  vicinity  of  the  medium  shock  pressure. 

Bonding  Gage  to  Medium 

A  practical  consideration  which  has  not  received  sufficient  attention  to  date  is  that  of  bonding  the 

pressure  gage  to  the  medium  at  its  emplacement  position.  A  grout  mixture  is  desired  which  will  be 

convenient  to  use  and  will  make  a  good  bond  or  contact  between  the  medium  and  the  gage  in  its  drill  hole, 

More  important,  a  grout  mixture  is  required  whose  Hugoniot  matches  that  of  the  medium.  If  the  grout 

has  a  shock  impedance  much  different  from  that  of  the  medium  then  a  significant  perturbation  may 

be  made  on  the  pressure  measurements.  The  preparation  of  grout  mixtures  with  Hugoniots  matching 

those  of  a  given  rock  medium,  at  least  over  a  small  range  of  pressures  where  a  measurement  is  con- 

17 

templated,  is  felt  to  be  a  relatively  easy  task.  An  analysis  of  the  synthesis  of  Hugoniots  of  multi- 
component  materials  has  indicated  how  the  grout  preparation  problem  may  be  approached. 

Shock  Pressure  Measurements  About  Nuclear  Explosions 

The  impedance-mismatch  gage  as  described  above  has  never  been  used  to  measure  shock  pressures 
in  a  field  experiment.  The  principal  reasons  for  this  have  been  (1)  that  no  adequate  electronic  circuitry, 
such  as  the  TDCM.  was  available  for  proper  recording  of  pressure  gage  signals,  and  (2)  that  no  oppor¬ 
tunity  has  arisen  where  an  impedance- mismatch  gage  could  be  placed  near  an  underground  explosion  for 
a  meanini^  evaluation  of  its  performance.  Now  that  the  impedance- mismatch  gage  has  been  developed, 
it  is  readly  for  field  trials.  Measurements  of  peak  ahock  pressure  wUl  be  attempted  at  the  1-,  0.6-,  and 
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0.4-megabar  pressure  levels  using  impedance-mismatch  gages  in  a  scheduled  nuclear  explosion  (Shoal)  in 
granite. 

As  a  result  of  research  conducted  on  the  impedance-mismatch  gage,  simplified  shock  pressure  gages 
have  been  fabricated,  tested,  and  used  in  field  experiments.  Shock  pressures  created  by  buried  nuclear 
explosions  have  been  successfully  measured  with  these  gages.  However,  in  each  case  where  a  simplified 
gage  was  used,  no  absolute  measure  -  in  the  sense  of  that  derived  from  an  impedance-mismatch  gage  -  of 
shock  pressure  was  possible.  Rather,  to  derive  shock  pressures  from  measurements  it  was  necessary 
to  make  the  assumption  that  the  medium  Hugoniot  was  known.  The  desirability  of  a  gage  which  gives 
direct  measurements  of  shock  pressure  is  clear,  for  with  this  type  of  gage  the  reliance  on  possibly  erro¬ 
neous  assumptions  regarding  medium  Hugoniots  is  avoided  and,  to  a  certain  extent,  the  requirement  of 
extensive  laboratory  determinations  of  medium  Hugoniots  is  eliminated. 

To  date,  23  attempts  have  been  made  with  simplified  gages  to  obtain  peak  shock  pressures  in  the 
hydrodynamic  region  about  nuclear  explosions.  Except  for  one  in  granite  rock,  all  measurements  were 
made  in  desert  alluvium,  a  porous,  weakly  consolidated  soil.  Of  the  23  measurements,  the  one  in  granite 
and  14  in  desert  alluvium  are  considered  reliable  and  the  remainder  are  felt  to  be  of  doubtful  or  little  value. 

In  desert  alluvium,  shock  pressures  were  measured  using  Plexiglas  gages  placed  in  vertical  drill 
holes  to  the  same  depth  as  that  of  the  explosion.  The  choice  of  Plexiglas  for  gage  material  was  dictated 
by  the  fact  that  it  is  not  possible  to  machine  and  prepare  a  sample  of  poorly  consolidated  alluvial  soil  for 
use  in  8  pressure  gsge  and  by  the  fact  that  the  Hugoniot  of  Plexiglas  matches  t'easonably  well  that  of  desert 
alluvium  over  a  workable  range  of  pressures. 

The  difficulty  of  not  being  able  to  machine  and  work  soil  cores  will  probably  persist,  and  any  shock 
pressure  gage  used  in  soli  measurements  will  probabiy  require  a  subatitute  material  which  can  be  worked 
and  whose  Hugoniot  la  as  almUar  as  possible  to  that  of  the  soil.  In  view  of  this,  it  will  be  necessai^  to 
make  laboratory  determinations  of  soil  Hugoniots  so  that  suitable  aubstitute  ntaterials  can  be  selected. 
Materials  which  can  serve  as  substitutes  for  soil  mtty  be  obtained  by  reviewing  all  available  data  on 
Hugoniots  In  the  hope  of  finding  one  which  is  nearly  the  sante  as  that  of  the  soil  medium,  or  by  synthe- 
aising*^  materials  so  that  their  Hugoniots  match  the  noil  Hugooirt  at  least  over  a  range  of  pressures 
where  measurements  sre  contemplated. 

The  reaaooahly  eloae  match  between  Hufooiote  of  Plexiglas  and  deaert  alluvium  ia  shosm  In 
Figure  I.IS.  From  meaauremeota  of  ahock  velocity  In  PlexiflBa.  shock  preaaure  In  Plantglas  la  obtained 
(aolutlon  of  Squationa  1.1  and  1.6).  Thla  preaaure  Uwn  Is  taken  to  be  the  abook  preaaure  la  deaert 
alluvium  at  tha  posUton  of  the  Plaidglaa  gage. 
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Two  types  of  gages  have  operated  successfully  in  desert  alluvium,  one  a  rectangiilar  slab  gage  and 
the  other  a  ring-type  gage.  The  Plexiglas  slab  gage  is  shown  in  Figure  1.14.  One  or  more  pairs  of  PZT 
shock  detectors  identical  to  those  of  Figure  1.4  were  employed  in  each  Plexiglas  gage  of  Figure  1.14.  A 
pair  of  detectors  were  accurately  spaced  about  2  inches  apart  along  the  directim  of  shock  propagation 
and  spaced  about  1  inch  apart  along  a  direction  transverse  to  that  of  shock  wave  propagation. 


Figure  1.14  Front  and  Side  Views  of  Ftexiglas  Sab  Cage  Used  in  Desert  Alluvium 


Attacbed  to  and  above  the  slab  gage  are  a  gyro«eo|>e  and  a  motor  to  relate  and  position  the  gage  near 


the  bottom  of  its  verileat  heia.  With  gjrroacepe  drift  taken  Into  aeeoimt  the  gage  may  be  aUgned  aalmulh- 
aUy  to  wfthlA  I  degree.  After  the  slab  gage  la  poshioned  the  AriU  hole  ia  ftUed  with  native  aoU  from  the 
same  area  ia  arhldi  the  gage  hole  *as  drilled. 
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To  avoid  the  great  inconvenience  of  orienting  the  slab  gage,  a  ring  gage  has  beat  develt^ed  which 
requirea  no  azimuthal  alignment.  The  manner  of  ring  gage  aeration  ia  illuatrated  in  the  diagram  of 
Figure  1.15.  Two  cylindrical  ringa  PZT  are  the  ahock  detecting  elementa.  Electrical  leada  are 
attached  to  the  inner  and  outer  radii  of  each  PZT  ring  and  taken  cut  to  the  aurface  through  the  axially 
located  ttd>e.  Detaila  of  conatruction  and  c^ration  may  be  found  in  Reference  18.  In  Figure  1.16  ia  ahown 
a  raater  oacillogram  of  ring  gage  aignala  obtained  in  a  field  experiment  with  a  buried  exploaion  in  deaert 
alluvium.  The  amplitude  of  the  pulaea  derived  from  the  PZT  rlnga  ia  about  300  volta.  Theae  pulaea  were 
obtained  at  a  position  where  the  ahock  wave  preaaure  waa  150  kilobara.  Timing  pipa  on  the  record  occur 
at  1  paec  intervale. 


Figure  1.15  Diagram  of  Preaaure  Cage  Whh  PZT  ftinga 


Shock  pressure  data  obtained  to  date  from  the  two  types  of  gages  used  in  desert  alluvium  are  plotted 
in  Figure  1.17  as  a  function  of  scaled  distance  from  explosion  center.  Each  open  circle  point  represents 
a  shock  pressure  measurement  from  one  explosion.  The  triangles  and  squares  occurring  in  pairs  signify 
two  pressure  measurements  from  a  single  explosion.  The  solid  curve  of  Figure  1.17  is  obtained  by  calcu¬ 
lation  using  methods  described  in  Reference  8.  Data  of  Figure  1.17  have  also  been  plotted  in  Figure  4.11 
of  this  report  for  comparison  with  Bishop's  calculations. 

The  reliability  of  shock  pressure  measurements  in  desert  alluvium  may  be  ascertained  by  compari¬ 
son  of  data  with  the  dotted  curve  of  Figure  1.17.  This  dotted  curve  is  an  indirect  and  less  accurate 
indication  of  shock  pressure  versus  distance  obtained  from  independent  measurements  of  shock  wave 
position  as  a  function  of  time. 

Shock  time-of-arrival  data  obtained  from  buried  explosions  in  desert  alluvium  are  shown  scaled  in 
Figure  1.18.  The  solid  curve  of  Figure  1.18  is  calculated  by  methods  given  in  Reference  8  and  the  dotted 
curve  is  a  best  fit  to  all  the  data.  An  analytical  expression  was  obtained  for  the  dotted  curve  from  which, 
by  differentiation,  shock  velocity  as  a  function  of  distance  from  explosion  center  was  estimated.  Since 
shock  velocity  is  now  known  at  any  given  position,  the  shock  pressure  at  any  position  can  be  determined 
from  the  Hugoniot  curve  (Figure  1.13  ,  1.20,  or  1.21)  of  the  medium.  Shock  pressure  as  a  function  of 
distance  from  the  explosion  estimated  in  this  manner  from  shock  time  of  arrival  data  is  the  dotted  curve 
shown  in  Figure  1.17.  Since  direct  shock  pressure  measurements  are  consistent  with  the  independent 
shock  wave  time-of-arrival  measurements,  it  is  believed  that  the  shock  pressure  data  are  reliable  and 
descriptive  of  spherical  shock  wave  propagation  in  desert  alluvium. 

Data  of  Figures  1.17  and  1.18  were  obtained  from  various  explosions  among  which  the  maximum 
ratio  of  two  explosion  energies  was  in  excess  of  50.  Since  data  of  Figures  1.17  and  1.18  scatter  nicely 
about  a  single  curve  in  each  figure,  the  measurements  verify  that  spherical  shock  phenomena  scale  as 
the  cube-root  of  energy  release,  W,  expressed  in  kilotons. 

Scatter  in  the  data  of  Figures  1.17  and  1.18  is  attributed  to  lack  of  precise  values  of  explosion 
energy  release,  to  experimental  errors  in  measurement,  to  the  possibility  that  all  explosions  were  not 
perfectly  spherical,  and  also  to  the  fact  that  the  medium  Hugoniot  was  very  likely  not  the  same  for  all 
explosions. 
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SHOCK  PRESSURE  (KILOBARS) 


Fifur*  1.17  P»ak  Shock  PreMure  Vertiw  Scaled  Dt«t«iic« 
Prom  Exploaiooa  In  Doaort  Ailuvium 
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Pf«Mr«  l.lt  Sr«lMl  V*r«Mi  Sr«l«4  Arrtvmi  TU***  for  Spt»«rical  Shock  Wave*  in  Deaert  Alluvium 


There  is  good  evidence  (see  discussion  below  on  Hugoniots  of  porous  soils)  to  indicate  that  relatively 
small  changes  in  the  porosity  or  water  content  of  desert  alluvium  or  any  other  porous  medium  make  sig¬ 
nificant  changes  in  the  Hugoniot.  Shock  pressure  gages  used  to  obtain  the  data  of  Figure  1.17  were  made 
of  Plexiglas  whose  Hugoniot  is  similar  to  that  of  completely  dry  desert  alluvium.  It  is  possible  that  the 
Hugoniot  of  Plexiglas  was  no  longer  a  good  match  for  the  Hugoniots  of  desert  alluvium  in  those  experiments 
where  the  porosity  and  moisture  content  of  the  medium  were  appreciably  different  from  those  assumed 
(Figui'e  1.13  or  1.20).  If  this  is  true,  then  the  data  of  Figure  1.17  will  require  some  correction  after  the 
influence  of  porosity  eind  moisture  content  on  the  Hugoniot  of  desert  alluvium  is  known.  Another  observa¬ 
tion  that  water  content  may  significantly  alter  the  medium  Hugoniot  has  been  made^^  by  correlating  time-of- 
arrival  data  from  various  experiments  with  some  gross  measure  of  water  content  in  the  medium  obtained 
from  drilling  logs.  If  experiments  are  grouped  as  "dry,*"  “damp,''  and  "average"  according  to  moisture 
condition  indicated  for  the  medium,  the  three  separate  distance-time  curves,  similar  to  the  one  shown  in 
Figure  1.18,  are  obtained  from  the  data  for  desert  alluvium.  A  best  fit  made  to  each  of  these  three  groups 
of  data  reveals  a  smaller  standard  deviation  than  that  for  the  average  curve  given  in  Figure  1.18  when  all 
data  are  taken  together.  The  interpretation  of  this  result  is  that  moisture  content  in  the  medium  does 
noticeably  alter  the  medium  Hugoniot. 

A  somewhat  different  modification  of  the  impedance-mismatch  gage  was  used  to  measure  peak  pres¬ 
sure  of  a  shock  wave  generated  by  a  nominal  5-kiloton  nuclear  explosion  (Hardhat)  in  granite.  For  con¬ 
struction  of  the  gage  a  core  of  granite  was  obtained  from  near  the  point  of  measurement  and  machined  into 
a  circular  disk  7  inches  in  diameter  and  2  inches  thick.  The  granite  disk  was  placed  in  direct  contact  with 
a  7-inch-diameter,  2-1 /2-inch-thick  disk  of  (2024)  aluminum  as  shown  in  Figure  1.19.  Shock  transit 
times  through  each  disk  were  recorded  with  PZT  crystals  placed  at  the  granite  front  surface,  where  con¬ 
tact  with  the  shock  was  first  made,  at  the  granite-aluminum  interface,  and  at  the  aluminum  back  surface. 
Three  sets  of  three  crystals  positioned  in  this  manner  and  on  axes  120  degrees  apart  were  installed  as 
insurance  against  the  possible  failure  of  one  or  more  crystals.  If  at  least  two  of  the  three  sets  of 
crystals  functioned  without  failure,  then,  in  addition  to  shock  velocity  in  each  disk,  the  tilt  of  the  shock 
wave  with  respect  to  the  longitudinal  axis  of  the  circular  disks  would  be  obtained  from  Equation  1.8.  One 
cable  for  each  of  the  three  sets  of  crystals  was  used  to  bring  out  the  signals. 

The  gage  assembly  was  inserted  at  the  bottom  of  a  drill  hole  whose  axes  made  an  angle  of  36  degrees 
with  a  line  joining  the  gage  position  and  the  explosion  center.  The  hole  was  then  filled  with  a  grout  whose 
sonic  impedance  matched  that  of  granite.  Shock  impedance  of  the  grout  mixture  was  not  known. 
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Insufficient  time  was  available  prior  to  the  e:q>eriment  to  prepare  a  grout  whose  Hugoniot  matched  that  of 
granite  at  the  pressure  level  e]q>ected. 

Because  the  instrument  drill  hole  was  not  truly  radial  it  was  desirable  to  obtain  a  measurement  of 
shock  tilt.  R  was  for  this  purpose  that  the  three  sets  of  three  shock  detectors  were  employed. 

Results  of  the  granite  shock  pres;Aire  experiment  are  shown  in  Table  1.1.  From  measured  shock 
velocities  pressures  were  determined  using  Equations  1.1  and  1.6. 

TABLE  1.1 

Shock  Pressures  From  a  Nuclear  Explosion  in  Granite 


E>isk 

Density 

(gm/cm^) 

Shock  path 
(cm) 

Transit  time 
(fisec) 

Shock  velocity 
(mm /p  sec) 

Pressure  in  disk 
(kb) 

Granite 

2.68 

4.07 

6.27 

6.48 

450  ±  50 

Aluminum 

2.79 

5.08 

6.10 

8.33 

514  ±  100 

The  shock  wave  pressure  measured  in  granite  during  the  Hardhat  explosion  is  consistent  with  an 
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independent  measurement  made  by  Lombard  at  a  nearby  position. 

It  must  be  noted  again  that  the  pressures  listed  in  Table  1.1  are  not  direct  measurements  but  were 
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deduced  from  shock  velocities  by  assuming  that  the  Hugoniots  for  granite  and  aluminum  were  known. 

The  purpose  of  placing  ari  aluminum  disk  behind  the  granite  disk  is  now  seen  to  be  an  attempt  to 
check  the  validity  of  the  assumption  regarding  the  granite  Hugoniot.  If  the  Hugoniot  assumed  for  granite 
were  tne  same  as  that  for  the  granite  disk  in  the  gage,  then  the  reflection  Hugoniot  for  granite  (Equation 
1.7),  passing  through  the  measured  granite  pressure  point,  would  intersect,  in  a  P-u  plot  such  as 
Figure  1.2,  the  accurately  known  aluminum  Hugoniot  (Equation  1.6)  at  the  pressure  measured  in  the 
aluminum  disk.  Apparently  the  Hugoniot  assumed  for  granite  was  correct  since  the  514-kilobar  pressure 
measured  in  aluminum  was  within  experimental  error  of  the  540  kilobars  indicated  by  the  solutions  of 
Equations  1.5  and  1.6. 

For  a  more  comprehensive  description  of  the  hydrodynamic  measurement  performed  on  the  Hardhat 
experiment.  Reference  8  should  be  consulted. 

Aside  from  experimental  errors  it  is  seen  from  the  measurements  of  shock  pressure  performed 
thus  far  that  the  greatest  uncertainty  associated  with  the  data  resides  in  the  assumption  that  the  medium 
Hugoniot  is  well  known.  Accuracy  of  the  shock  pressure  data  from  experiments  in  desert  alluvium 


d^ends  on  the  validity  of  the  assumpticm  that  the  Hugoniot  of  u  situ  desert  alluvium  is  the  same  as  the 
laboratory  determined  Hugmuot  of  dry  reconstituted  alluvium  samples  and  nearly  the  same  as  the  Hugoniot 
of  the  Plexiglas  used  in  gages.  As  pointed  out.  there  are  strong  indications  that  these  assumptions  for 
desert  alluvium  are  not  absolutely  correct  because  the  moisture  content  of  the  medium  was  not  zero  for 
all  the  field  experiments  but  apparently  varied.  This  variatiwi  from  zero  moisture  content  results  in  a 
Hugoniot  different  from  that  for  a  perfectly  dry  material.  Accordingly,  the  desert  alluvium  pressure  data 
have  inherited  a  certain  unknown  inaccuracy  resulting  from  the  Hugoniot  assumption. 

The  granite  pressure  of  Table  1. 1  is  consistent  with  the  pressure  measured  in  aluminum  and  indi¬ 
cates  that  within  experimental  error  the  Hugoniot  assumed  (least  square  fit  of  available  data)  for  the  m 
situ  medium  was  correct.  However,  had  the  granite  and  aluminum  pressures  been  widely  inconsistent, 
no  pressure  data  would  have  been  obtained  and  the  only  information  provided  by  experimental  results 
would  be  that  the  assumed  Hugoniot  was  probably  incorrect. 

If  there  is  a  high  degree  of  confidence  that  a  medium  Hugoniot  is  well  known  then  there  is  absolutely 
no  need  for  attempting  elaborate  and  difficult  pressure  measurements.  All  that  is  required  is  a  measure¬ 
ment  of  shock  wave  arrival  time  at  various  positions  from  an  explosion  center  to  obtain  a  distance-time 
curve  such  as  that  shown  in  Figure  1.18.  Since  the  medium  Hugoniot  is  considered  well  known  and  time- 
of-arrival  measurements  can  be  made  simply  and  accurately,  shock  pressure  may  be  deduced  at  various 
positions  from  the  explosion  as  was  done  to  obtain  the  dashed  curve  of  Figure  1.7,  However,  it  seems 
unlikely  that  we  will  ever  have  this  much  confidence  in  a  Hugoniot  assumed  for  any  geologic  material 
in  situ. 

To  make  our  shock  pressure  measurements  more  meaningful,  it  is  important  to  divorce  the 
measurement  technique  from  any  assumptions  regarding  the  medium  Hugoniot  and  to  attempt  some 
relatively  direct  measurements  of  shock  pressure.  Hopefully,  the  Impedance  mismatch  pressure  gage  can 
accomplish  this  task,  at  least  In  rock  type  geologic  media. 

Hugonlots  of  Porous  Earth  Materials 

In  conjunction  with  development  of  shock  pressure  gages  for  use  In  soil  media,  an  examination  of 
soil  Hugonlots  has  been  necessary.  Since,  as  mentioned  previously,  cores  of  poorly  consolidated  soils 
cannot  easily  be  worked  and  prepared  for  fabrication  of  shock  pressure  gages,  a  particular  material  whose 
Hugoniot  is  very  nearly  the  same  as  that  of  the  soil  in  question  is  substituted  in  the  gage.  After  laboratory 
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determination  of  the  soil  Hugoniot,  the  properties  required  by  the  substitute  material  Hugoniot  are  made 
known.  Hugoniot  data  on  soils  are  also  necessary  for  performing  calculations  of  shock  wave  propagation 
in  these  media. 

Data  have  been  obtained  on  the  Hugoniots  of  desert  stUuvium,  volcanic  tuff,  and  a  medium  called  playa. 
Volcanic  tuff  is  a  porous  rock,  well  consolidated  (relatively  well  cemented  and  firm).  Desert  alluvium  is  a 
very  poorly  consolidated  (friable)  soil,  and  playa  is  an  alluvial  soil  relatively  well  consolidated.  These 
materials  are  representative  of  the  media  in  which  shock-pressure  measurements  are  currently  being 
made.  Core  samples  of  desert  alluvium  and  of  txiff  which  underlies  this  alluvium  were  obtained  at  various 
depths  from  Area  3  of  the  Nevada  Test  Site  (NTS).  Playa  samples  were  obtained  from  Frenchman  Flat, 
Area  5  of  NTS. 

Sample  preparation  of  these  materials  and  the  techniques  employed  to  obtain  Hugoniot  data  are  given 
20 

in  another  report.  Majority  of  the  data  are  from  samples  that  were  completely  dry,  i.e. ,  any  moisture 
contained  in  the  samples  obtained  from  cores  was  removed  by  baking  before  Hugoniot  experiments  were 
conducted.  The  remainder  of  the  data  are  results  of  investigations  into  the  effects  of  porosity  and  water 
content  on  the  Hugoniot  of  the  dry  material.  Data  are  summarized  in  Table  1.2  at  the  end  of  this  section. 

In  Figure  1.20  are  plotted  data  for  dry  desert  alluvium.  Porosity  of  this  material  is  about  30  per- 
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cent,  and  the  initial  density  of  samples  ranges  between  the  values  of  1.38  and  1.80  gm/cm  .  The  squares 
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and  triangles  of  Figure  1.20  represent  data  from  the  earlier  work  of  McQueen  and  Marsh.  Shock  wave 

velocity  and  particle  velocity  behind  the  wave,  plotted  in  the  figure,  indicate  raw  data  from  which  shock 

pressure  and  density  may  be  determined  using  Equations  1.1  and  1.2,  It  Is  pointed  out  that  the  squares 

and  triangles  of  the  figure  are  data  obtained  from  soil  samples  taken  from  one  location  in  the  medium. 

Circles  represent  data  from  samples  taken  at  different  locations  in  the  same  medium,  which  probably 

account'  for  the  scatter  in  the  data,  since  the  alluvium  is  highly  inhomogeneous.  Also  shown  in 
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Figure  1.20  for  comparison  is  a  line  through  Hugoniot  data  for  aluminum  and  a  line  estimated  for 
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nonporous  desert  alluvium.  A  synthesized  Hugoniot  for  dry  desert  alluvium  of  40  percent  porosity 
is  also  plotted  on  the  figure. 

Curves  may  be  drawn  through  the  data  of  Figure  1,20  in  a  number  of  ways.  We  have  drawn  two 

curves  to  designate  the  range  of  densities  likely  to  be  encountered  in  a  desert  alluvium  me  dium  and  for 

which  there  are  some  data.  It  is  seen  that  most  of  the  data  fall  within  that  region  of  the  figure  bordered 

by  the  two  curves,  Hugoniot  data  and  the  curves  of  Figure  1.20  have  been  r^lotted  in  Figure  1,21  to 

show  shock  pressure  as  a  functicm  of  relative  specific  volume,  v/v^  ,  where  v  is  shocked  volume  of  the 
medium  and  v^  is  initial  specific  volume.  For  the  dry  samples  of  alluvium  whose  initial  density,  , 
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was  1.8,  the  derived  grain  doislty  was  2.46  gm/cm  ,  and  for  the  samples  of  density  ■  1.54,  grain 

3 

density  is  givai  as  2.24  gm/cm  . 

By  conq>arison  with  Hugoniots  of  most  other  (nonporous)  materials  the  Hugoniots  of  desert  alluvium 

are  most  unusual.  It  is  seen  that  the  alluvium  Hugoniots  are  double  valued  functions  of  density  and  that 

the  Hugoniot  has  a  region  of  pressures  over  which  the  slope.(8P/  8p),  of  the  curve  is  negative.  Hugoniots 
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of  porous  metals  have  been  found  *  which  demonstrate  similar  characteristics. 

A  qualitative  explanation  of  the  shape  of  alluvium  Hugoniot  curves  may  be  as  follows.  Ih  porous 
soil  individual  particles  of  density  p^  are  separated  by  void  spaces  so  that  the  average  density  is  p^.  A 
shock  wave  passing  through  the  material  does  work  by  collapsing  the  void  spaces,  by  overcoming  the 
internal  friction  between  particles  and  by  deforming  and  compressing  the  particles.  The  increase  in 
internal  oiergy  (Equatim  1.3)  as  a  result  of  shock  compression  is  greater  for  a  porous  medium  than  for 
a  similar  nonporous  medium  when  shock  pressures  in  each  medium  are  equal.  This  may  be  seen  from  the 
diagram  of  Figure  1.22. 
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Figure  1.22  ComiNiriscsi  of  Energy  Deposition  by  Shock  Waves  In  Porous  and  Nonporous  Solids 

The  increase  In  apectflc  Internal  aneigy  given  by  Equation  t.S  is  represented  in  Figure  1.22  by  the 
triangular  area,  «  tor  the  nonporous  material  and  is  less  than  the  area  tor  the  porous 


material.  Since  the  increase  in  internal  energy  behind  the  shock  wave  is  equal  to  the  increase  in  material 
kinetic  energy  (from  Equations  1.1  and  1.2)  it  is  also  seen  that  the  kinetic  energy  Imparted  to  a  porous 
material  will  be  greater  than  that  imparted  to  a  similar  nonporous  material  when  shock  pressures  in  both 
materials  are  equal.  It  is  for  this  reason  that  porous  Hugonlot  curves  such  as  shown  in  Figure  1.20  will 
be  below  the  Hugoniot  of  the  nonporous  material.  Consequently,  in  porous  materials  similar  to  desert 
alluvium,  internal  energy  and  particle  velocity  will  be  greater  and  shock  wave  velocity  will  be  less  than 
that  produced  by  a  shock  wave  of  the  same  pressure  in  the  material  with  no  void  spaces.  Because  internal 
energy  in  the  porous  material  is  greater  than  in  the  noiqmrous  material,  temperatures  behind  the  shock 
waves  in  porous  materials  will  be  greater.  These  higher  temperatures  in  porous  materials  may  have 
effects  such  as  shock*  induced  phase  changes  at  pressures  lower  than  those  at  which  phase  changes  occur 
in  nonporous  materials.  Also  the  higher  temperatures  resulting  in  porous  media  from  shock  compression 
may  account  for  the  region  of  negative  slope.  (tP/gp),  demonstrated  by  the  Hugoniots  in  Figure  1.21. 

The  shape  of  the  alluvium  Hugoniots  in  Figure  1.21  between  0  and  100  kilobars  is  possibly  described 
by  the  collapsing  of  void  spaces,  the  principal  mechanism  of  medium  deformation  in  this  pressure  range. 

A  small  change  in  shock  pressure  results  in  a  large  change  in  specific  volume.  Between  about  100  and 
200  kilobars  phase  changes  are  being  initiated  by  shock  waves  passing  over  particles  of  the  medium.  At 
pressures  in  excess  of  200  kilobars  the  effect  of  shock  heating  is  sufficiently  intense  in  the  porous  material 
that  increases  in  shock  pressure  are  acconq>anied  by  decreases  in  density.  Also  In  this  pressure  region 
the  two  alluvium  Hugoniot  curves  appear  to  be  coaloacing.  i.  e. .  aa  shock  pressure  increases  above  about 
300  kilobars  fewer  and  fewer  void  spaces  remain,  and  a  given  ahock  preaaure  coropressea  both  materials, 
of  different  initisl  density,  to  very  nearly  the  same  density.  Eventually,  at  preaaures  of  the  order  of  half 
a  megabar,  both  poroua  and  nonporous  materials  are  compressed  to  the  same  density  by  a  shock  wave  of 
a  given  preaaure.  It  may  be  that  at  tome  preaaure  greater  than  about  0.S  megsbar  the  Hugoniot  curves  tit 
porous  alluvium  are  tdentlcsl  to  a  Hugoniot  of  a  nonporous  alluvium  medium,  h  must  be  remarked  that, 
in  all  the  alluvium  experlmenta  conducted,  apace  in  the  samples  was  occupied  hy  air  at  atmoapherte 
pressure.  What  taflueace  this  amail  mass  of  gaa  may  have  on  the  Hugoniot  la  not  actually  known. 

While  these  cpialttative  apeculetions  may  net  be  psftlcularly  dliuninatbig  they  serve  to  illustrate 
the  complexity  of  phenomena  involvad  and  the  (Nffleultlea  likely  to  be  encountered  in  fermulatlng  theoretical 
deacHpttona  of  poroua  soil  Ibiganiels. 

The  (kSmiegabar  Hugoniot  point  for  desert  altuvlum  represents  tha  upper  Umtt  of  ahock  preaauraa 
attalnabla  to  laboratory  axparlmanta.  M  is  practically  impoaaihla  to  ganorata  praaaurea  much  hii^r  than 
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thla  conventional  methods.  As  a  consequence  ot  this  limitation  additional  Hugoniot  information  in  the 

pressure  region  above  0.5  megabar  must  cone  from  other  sources  such  as  e:q>eriments  about  nuclear 

explosions  or  theoretical  investigations.  While  experiments  with  nuclear  explosions  will  be  considerably 

more  difficult  than  laboratory  experiments,  shock  pressures  as  high  as  10  or  even  100  megabars  may  be 
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achieved  and  valuable  Hugcmiot  and  thermo<^amlc  Information  can  be  obtained.  Attempts  ‘  ‘  • 

have  been  made  to  calculate  Hugoniots  of  materials  for  shock  pressures  between  0. 1  and  1  megabar:  how> 
ever,  the  successes  of  these  methods  have  not  been  striking  and  great  reliance  cannot  be  placed  on  the 
calculated  results.  Much  more  work  is  required  in  the  theoretical  calculation  of  Hugoniots  for  geologic 
media,  particularly  those  media  with  high  void  content. 

An  implication  of  the  alluvium  Hugoniot  data  shown  in  Figure  1.31  is  that  Gruneisen's  parameter. 

r  •  v(gP/tfE)^  .  for  desert  alluvium  is  not  a  constant  but  likely  a  strongly  varying  function  of  specific 

volume.  V  >  ^  .  This  implication  is  significant  to  theoretical  calculations  of  spherical  shock  wave  props* 
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gation  in  desert  alluvium  or  other  porous  geologic  media.  In  these  theories  *  the  Mie-Grunelsen 
equation  of  stale,  which  contains  the  parameter  f.  is  commonly  invoked  and  Gruneisen's  T  is  usually 

taken  to  be  constant.  While  for  many  materials^*  T  is  nearly  constant  or  only  a  slowly  varying  function 

of  specific  volume,  the  data  for  alluvium  indicate  that  f  is  not  constant  uut  varies  appreciably  with 

specific  volume. 

Gruneisen's  parameter  for  solids  is  analogous  to  the  ratio  of  specific  heats  for  gases  and  determines 
how  rapidly  pressure  falls  off  behind  the  shock  front  and  the  amount  of  cooling  eiqierienced  by  a  medium 
particle  in  adiabatic  expansion.  These  processes  are  important  in  determining  how  rapidly  peak  shock 
pressure  falls  off  with  distance  and  the  rale  at  which  a  shock  wave  is  propagated  from  the  explosion 
source. 

From  the  alluvium  data  of  Figure  1.31  estimates  of  Gruneisen's  f  and  Its  variation  with  volume  may 
be  made  in  a  number  of  ways.  For  example  from  the  measured  Hugoniot  and  the  Dugdale*  MacDonald  or 
Sister  relations^  for  f.  values  of  f  for  desert  alluvium  may  be  computed.  Also  estimates  of  f  may  be 
obtained  using  the  method  of  Altshuler**  for  porous  iron,  evaluations  of  Crunetsen's  parameter  for 
desert  aUuvitnn  by  these  various  methods  results  in  a  confusing  array  of  values. 

■  *4  «a 

Thsoretlcsl  cslcwlstlens  *  *  of  sphsrlcsl  shoek  wave  propsgstioos  In  desert  sUuvIwn  with 

constant  value  of  f  (both  0.5  and  I)  have  yisUM  predictions  of  shock  position  versus  time  trhich  are 
generally  la  exceUsnt  sgrssmsnt  with  axpartmanlal  data  (a.g. .  ass  Figure  l.lt  and  also  Figure  t  of 
Itefermwe  15).  Stuee  the  deoert  aUmrium  Hugoniot  data  hidteata  that  r  ia  daflaltaly  aot  cenManl.  one 


wonders  how  the  spherical  shock  calculations  can  give  correct  answers  with  constant  values  of  r .  It  is 
not  known  whether  agreement  between  calculations  and  data  is  fortuitous,  whether  shock  position  versus 
time  is  not  sensi  ive  to  large  changes  in  F  demonstrated  alluvium,  or  whether  the  Mle-Gruneisen 
equation  of  state  is  not  applicable  to  a  porous  soil  medium.  However,  in  order  to  have  confidosce  in  the  cal¬ 
culations  the  r  variation  indicated  by  desert  alluvium  must  be  understood  nd  its  implicatitms  for  theoret¬ 
ical  results  must  be  evaluated  more  fully. 

Figure  1.23  shows  data  obtained  for  a  volcanic  tuff  from  Areas  3  and  16  of  the  Nevada  test  site. 

3 

The  dry  tuff  has  a  density  of  1.46  gm/cm  and  a  porosity  of  37  percent.  These  Hugonlot  ctota  are  different 

2 

from  those  reported  by  Lombard  for  tuff  from  Area  12  at  NTS.  Lombard's  tuff  had  a  density  of  1.7  and 
a  porosity  of  about  34  percent.  His  data  points  for  dry  tuff  lie  above  those  shown  for  the  dry  tuN  in 
Figure  1.23.  Data  for  saturated  tuff  are  also  shown  in  Figure  1.33. 

Figure  1.34  presents  Hugonlot  data  for  the  playa  medium  of  Area  S  at  the  Nevada  test  site.  ‘Fhe 
density  of  samples  received  from  the  field  was  1.47  Completely  dry.  these  samples  had  a 

density  of  1.41  gm/cm\  The  porosity  of  playa  is  estimated  to  be  about  30  percent. 

From  Figures  1.30.  1.33.  and  1.34,  it  is  observed  that  between  30  and  300  kllobars  the  Hugonlots  of 
dry  desert  alluvium,  dry  tuff,  and  plays  are  essentially  identical. 

bt  Figure  1.33  are  shown  the  results  of  a  few  experiments  with  a  fine,  pure  silica  sand.  These 
experiments  were  conducted  to  obtain  information  on  the  effects  of  porosity  snd  water  content  on  Hug^niota. 

1*he  silica  sand  used  is  compnsed  of  particles  10  percent  of  which  have  diameters  less  than  73 

microns.  Maximum  particle  sUe  is  about  130  microns.  Crain  density  of  the  sand  is  1.63  gm/cm  .  {tie 

same  as  that  of  crystaltine  quarts.  For  expeHments.  the  sand  was  {Kicked^**  to  dry  densities  of  1,38  and 

3.03  gm/cm  .  Corresponding  porosities  are  41  percent  and  33  percent.  It  is  seen  that  the  higher  density 

dry  sand  data  tie  above  those  of  the  lower  density  dry  ssnd  and  in  a  direction  toward  the  Itugoniot  curve 
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tor  crystalline  quarts  *  ahown  ai  the  top  of  Figure  1.33.  Aiae  shown  in  the  figure  are  data  far  water- 
saturated  sand  with  a  denaily  of  I.M.  Since  the  porosity  of  the  water-sMurated  aand  is  the  same  as  that 
of  the  loweM  curve  in  the  figure,  it  is  seen  thst  adtBtinn  of  water  to  fill  the  void  space  fat  this  s^  shifts 
the  Hiigooiot  upward  by  a  consldershle  amotMt.  making  the  saturated  sand  much  less  compreasible.  Thus 
a  decrease  In  porosity  or  sn  bicrense  In  wsler  content  of  s  pnrous  mMerial  results  in  a  higher  shock 
pressure  tf  the  same  rolativt  donsity  Is  scMsved  fsse  Figure  I.Hk 
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In  the  close-in  region  about  buried  explosions  in  porous  media,  this  means  that  stronger  than 
expected  shock  waves  will  be  propagated  if  the  water  content  of  the  medium  is  greater  than  expected  or  if 
the  porosity  is  less  than  expected.  For  experiments  with  explosions  in  porous  media  it  will  be  essential 
to  obtain  information  on  the  to  situ  porosity  and  water  content. 

Summary 

An  impedance-mismatch  pressure  gage  has  been  developed  which  is  capable  of  measuring  peak  shock 

wave  pressures  of  amplitude  0.1  to  1  megabar  in  the  hydrodynamic  region  about  buried  explosions.  In 

■i'  '  ■  r  . 

addition  the  gage  provides  information  on  the  medium  Hugo'niot  and  a  means  of  measuring  degree  of  tilt 
of  the  shock  wave. 

Techniques  employed  in  measuring  shock  wave  pressures  have  been  tested,  and  shock  pressures 
have  been  measured  in  experiments  with  buried  nuclear  explosions.  By  comparison  with  independent 
hydrodynamic  measurements  the  shock  pressure  data  have  been  shown  to  be  reliable.  These  data  are 
currently  being  employed  to  evaluate  and  refine  theoretical  calculations  of  spherical  shock  wave  propaga¬ 
tion  to  solid  media. 

Limited  investigations  have  been  made  of  Hugoniots  of  porous  geologic  media.  Some  unusual  features 
of  Hugoniots  for  porous  media  have  been  revealed  which  were  not  previously  observed  in  Hugoniot 
studies  of  other  (nonporous)  materials.  Moisture  content  and  porosity  are  found  to  be  significant  param¬ 
eters  influencing  the  Hugoniot  of  porous  materials.  Possible  applications  of  the  Hugoniot  data  for  porous 
media  to  theoretical  calculations  of  shock  propagation  to  these  media  have  been  pointed  out. 
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TABUS  l.a 


Hufoniot  Otte  tor  Pbrotw  Carlli  Matertola 


'b 

V 

(mm/asec) 

u 

<mm/|isec) 

P 

fkb) 

p  Ip 

Source 

Dry  Oeaert  AUuvlun 

-  Puroaity.  35  to  30% 

1.44 

a.66 

1.00 

38 

0.634 

NTS,  Area  3 

i.sa 

a.»o 

0.97 

44 

0.666 

NTS.  Area  3 

1.S4 

3.41 

1.58 

83 

0.537 

NTS,  Area  3 

1.50 

3.65 

1.57 

86 

0.570 

NTS,  Area  3 

1.70 

3.70 

1.53 

96 

0.589 

NTS,  Area  3 

1.4C 

4.35 

3.45 

156 

0.437 

NTS.  Area  3 

l.M 

4.36 

3.39 

171 

0.453 

NTS,  Area  3 

t.3« 

5.35 

3.37 

337 

0.377 

NTS.  Area  3 

1.77 

5.19 

3.37 

351 

0.438 

NTS,  Area  3 

Dry  VckleaiUc  Tuff  •  Beraaiiy,  37% 

1.4« 

1.34 

0.95 

31 

CL  576 

NTS.  Area  3  (1500>  det^h) 

1.46 

3.70 

1.58 

85 

0.573 

NTS,  Area  3  (|5M‘  depth) 

1.46 

4.38 

3.38 

143 

0.467 

NTS.  Area  3  (|500>  depth) 

l.4i 

4,80 

3.50 

153 

0.405 

NTS.  Area  3  (1500*  depth) 

1.46 

4.76 

3.90 

303 

0.393 

NTS.  Area  3  (1500*  depth) 

1.46 

3.68 

1.00 

39 

0.637 

NTS.  Area  16 

1.46 

3.56 

1.57 

83 

0.566 

NTS.  Area  16 

1.46 

4.03 

3.>0 

147 

0.380 

NTS.  Area  )6 

1.46 

4.34 

3.48 

153 

0.430 

NTS.  Area  )6 

Waiar  8aiurata4  Votranlr  Tuff 

l.T» 

3.43 

0.80 

53 

0.737 

NTS.  Area  )6 

1.74 

4.U 

1.45 

108 

0.660 

NTS.  Area  16 

1.15 

5,48 

3.81 

370 

%488 

NTS.  Area  S 

IS 


TABLE  1.2  (Cent) 


"o  3 

(gea/ctn) 

V  u 

(nun/fioec)  (mm/MSOc) 

P 

(kb) 

Source 

Plnya  •  Poroolly. 

*>'50%:  Water  Content, 

"•6% 

1.41« 

2.70 

1.04 

40 

0.615 

NTS.  Area  5 

1.41* 

4.40 

2.48 

148 

0.436 

NTS.  Area  5 

1.47 

3.00 

1.08 

48 

0.640 

NTS.  Area  5 

1.47 

2.58 

1.04 

39 

0.597 

NTS,  Area  5 

1.47 

3.69 

1.60 

87 

0.566 

NTS.  Area  5 

1.47 

4.47 

2.52 

166 

0.436 

NTS.  Area  5 

1.47 

4.36 

2.50 

160 

0.427 

NTS.  Area  5 

1.47 

5.07 

3.54 

264 

0.302 

NTS.  Area  5 

1.47 

5.24 

3.52 

271 

0.328 

NTS,  Area  5 

tbry  Silica  Sand  -  Poreotty.  40% 

l.&a 

3.13 

1.17 

58 

0.626 

l.M 

3.23 

1.16 

59 

0.641 

1.59 

3.42 

1.61 

88 

0.529 

i.sa 

3.47 

1.70 

93 

0.5)0 

1.56 

4.26 

2.25 

150 

0.472 

1,62 

4.24 

2.23 

153 

0.474 

Ory  Silica  Sand  -  Poreolly.  22% 

2.02 

3.45 

1.07 

75 

0.690 

2.14 

3.70 

1.46 

1)6 

0.605 

2.03 

4.78 

2.03 

197 

0.575 

Water  Saturated  Silica  Sand. 

PartMity. 

40% 

2.02 

4.53 

0.98 

90 

0.784 

1.97 

5.00 

1.45 

143 

0.710 

1.95 

5,63 

1.94 

2)3 

0.655 

2.00 

5.59 

1.93 

218 

0.655 

obtained  frt^  IMhfoe*  de|Mli  wore  rownMely  dry. 
Kenwteder  of  ntprrlmenu  wore  witli  oamidofl  from  abetd 
ll'foot  dofNb  »«  received  rroni  flold.  o^-'oled  in  wm  conulnero- 
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DYNAMIC  RESPONSE  OF  GEOLOGIC  SOLIDS  TO  LARGE  AMPLITUDE  STRESS  WAVES 

R.  C.  Bass,  H.  L.  Hawk,  and  A.  J.  Chabai 

Introduction 

Propagation  of  large  amplitude  stress  waves  in  solids  is  reasonably  well  understood  when  com¬ 
pressive  stresses  are  in  excess  of  about  100  kilobars.  At  these  stress  levels  material  rigidity  is 
negligible,  the  stress  wave  is  a  steep-fronted  shock,  and  a  fluid  model  is  used  to  describe  wave  propa¬ 
gation  in  the  solid.  For  small  amplitude  stresses,  say  less  than  about  1  kilobar,  an  elastic  model  is 
often  sufficient  to  describe  the  response  of  the  solid.  Between  these  limits,  when  amplitudes  of  com¬ 
pressive  stresses  are  of  the  order  of  1  to  100  kilobars,  there  is  no  general  model  which  effectively 
describes  the  dynamic  behavior  of  solids,  particularly  geologic  solids. 

In  this  range  of  nonfluid-nonelastic  stresses  phenomena  such  as  yielding,  fracture,  strain  rate 
and  relaxation  effects,  anomalous  compression  and  viscosity  are  expected  to  occur  and  may  play 
dominant  roles  in  determining  the  detailed  structure  of  the  compressive  wave  profile.  Principally 
because  an  adequate  model  does  not  exist  for  the  description  of  nonfluid-nonelastic  behavior  of  geologic 
solids,  there  is  at  present  no  theory  of  spherical  wave  propagation  from  buried  explosions  which  is 
capable  of  reproducing  wave  profiles  that  have  been  observed  in  experiments.  It  is  believed  that,  once 
a  reasonably  realistic  model  were  formulated,  theoretical  calculations  of  spherical  wave  disturbances 
could  be  made  which  would  yield  wave  profiles  similar  to  those  observed  in  experiments. 

If  a  model  can  be  obtained  allowing  accurate  calculations  of  nonfluid-nonelastic  stress  wave  propa¬ 
gation,  then  elastic  disturbances  resulting  fi’om  buried  explosions  In  geologic  media  can  be  predicted 
and  seismic  phenomena  may  thereby  be  better  understood. 

To  study  nonfluid-nonelastic  behavior  of  solids,  use  is  made  of  plane  compressive  waves  of  large 
amplitude,  o  ,  produced  by  explosives  or  by  high  velocity  projectile  Impacts.  In  these  experiments,  a 
one- dimensional  state  of  strain,  ^  =  (1-p^/p),  Is  achieved.  Often  the  single,  step-function  type  of 
wave  pulse  introduced  into  the  test  specimen  by  these  means  Is  unstable  because  the  condition  for 
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stability,  8  a/de  >  0,  has  been  violated.  As  a  result  two  waves  are  propagated  through  the  specimen. 
Simple  examples  of  some  stress-strain  curves  leading  to  double  wave  profiles  are  shown  in  Figure  2.1. 


IDEAL 

FLUID  MEDIUM 


a 


(a) 


(b) 


MEDIUM  WITH 
MIXED  PHASE 
REGION 


MEDIUM  WITH 

ANOMALOUS 

COMPRESSION 


Figure  2. 1  Some  Stress-Strain  Curves  Leading  to  Two-Wave  Profiles 


Figure  2.1(b)  is  an  example  of  a  stress-strain  relation  or  *Hugoniot"  (Equation  1,5  of  Section  1) 


for  an  ideal  solid  in  which  waves  of  any  amplitude  are  stable.  In  Figure  2.1(b)  is  illustrated  a 


stress-strain  curve  exhibiting  a  finite  yield  strength  or  "Hugoniot  elastic  limit.*  This  curve  gives  rise 
to  a  double  wave  structure  consisting  of  an  elastic  precursor  which  propagates  with  the  sonic  velocity 

of  the  specimen  and  a  second  higher  amplitude  wave  of  slower  velocity.  Double  wave  profiles  resulting 
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from  yielding  have  been  observed  in  metals,  '  '  in  quartz,  ’  and  also  in  rock  materials, 

A  double  wave  structure  as  shown  in  Figure  2.1<c)  may  also  result  from  polymorphic  phase 

transitions  induced  in  the  test  specimen  by  the  incident  wave.  The  phase  transition  mechanism  has  been 
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invoked  to  explain  double  waves  observed  in  iron,  quartz,  '  and  rocks,  '  ’  Phase  transitions 
induced  in  solids  by  shock  waves  seem  to  occur  generally  at  stresses  in  excess  of  100  kilobars.  How- 

7 

ever,  Gregson  and  Grine  suggest  that  in  carbonate  rocks  phase  transitions  may  occur  at  stresses  as 
low  as  30  kilobars. 

Figure  2.1(d)  illustrates  a  medium  with  a  stress-strain  curve  having  a  region  of  anomalous 

compressibility.  Over  this  range  of  stresses  a  steep-fronted  wave  cannot  be  propagated  and  becomes 

0 

elongated  as  it  travels  through  the  specimen.  Plane  wave  experiments  have  shown  that  fused  silica 
has  a  stress-strain  curve  similar  to  that  of  Figure  2.1(d). 

If  a  medium  is  characterized  by  viscous  or  dissipative  effects  the  rise  time  of  a  wave  front  such 
as  shown  in  Figure  2.1(a)  will  not  be  steep  but  will  appear  more  like  that  illustrated  in  Figure  2, 2(a). 

That  viscous  effects  may  be  significant  has  been  demonstrated  by  noticeable  variation  in  rise  times 
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observed  in  plane  wave  experiments  with  metals.  ' 

If  relaxation  phenomena  occur  in  conjunction  with  yielding,  wave  profiles  such  as  shown  in 
Figure  2.1(b)  may  appear  more  like  those  of  Figure  2.3(b).  Very  definite  relaxation  effects  have  been 
witnessed  in  Iron,^'*^  steels,^  and  quartz.^*® 

It  is  apparent  that  a  realistic  model  describing  nonfluid-nonelastic  behavior  of  geologic  solids 

may  be  quite  complex  since  a  natural  material  may,  for  example,  exhibit  several  of  the  phenomena 

illustrated  in  Figures  2,1  and  2.2,  Mechanisms  of  these  various  phenomena  are  largely  unknown,  oven 

for  metals.  For  comprehensive  discussions  of  wave  instability,  viscosity,  anomalous  compression, 

relaxation,  and  yielding,  the  reader  is  referred  to  the  literature,  particularly  two  excellent  review 
12  13 

articles  by  Duvall  *  and  the  references  contained  therein. 

The  extent  to  which  geologic  solids  demonstrate  nonfluid-nonelastic  effects  such  as  considered 
above  has  not  been  sufficiently  investigated.  The  purpose  of  the  work  described  here  has  been  to  Initiate 
experiments  by  which  these  various  effects  may  be  observed  and  to  perform  a  survey  of  a  few  rock 
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materials  with  large  amplitude  <1  to  100  kllobars)  stress  waves  of  high  strain  rates  (10  to  10  sec  ) 
under  conditicms  similar  to  those  created  by  buried  explosions.  These  experiments  were  designed  to 
yield  sufficient  data  to  at  least  characterize  the  responses  of  various  rock  materials  and  to  provide 
guiding  information  which  might  assist  in  formulations  of  reasonable  models.  More  specifically, 
objectives  were  to  investigate  in  rock  materials  of  current  interest  the  magnitude  of  the  dynamic  yield 
point  if  one  exists,  the  extent  of  viscous  properties  which  might  Influence  the  rise  times  of  propagated 
waves,  the  possibility  of  occurrence  of  anomalous  compression  in  rocks,  and  whether  or  not  rock 
materials  demonstrate  relaxation  effects. 


(a)  Wave  in  Viscous  Medium 


(b)  Wave  in  Medium  With  Nonzero  Elastic  Limit  Demonstrating 
Stress  Relaxation 


Figure  2.3  Wave  Profiles  for  a  Viscous  Medium  and  Wave  Profile  With  Stress 
Relaxation  After  Yielding 

b)  order  to  realize  these  aims  it  has  been  necessaiy  to  develop  an  instrument  capable  of  measuring 
large  amplitude  stresses  with  high  time  resolution. 
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Description  of  Quarts  Stress  Gage 
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The  new  X-cut*  quartz  crystal  technique  perfected  by  Neils<m  '  et  aL  for  examination  of 
stresses  in  metals  has  been  adqited  for  the  observation  of  stresses  in  geologic  solids.  Since  quartz 
has  a  higher  cfynamic  elastic  limit  than  that  expected  of  most  geologic  materials,  it  is  an  ideal  trans¬ 
ducer  for  investigaticHi  of  stress-time  wave  profiles.  The  (iynamic  elastic  limit  of  X-cut  quartz  has 
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been  rqx>rted  '  as  about  SO  kilobars.  X-cut  quartz  has  also  been  shown  to  have  a  linear  piezoelec¬ 
tric  response  to  at  least  25  kilobars. 

Quartz  transducers  used  to  obtain  stress-time  profiles  in  solids  have  been  found  to  compare 
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favorably  with  other  techniques  such  as  the  <9tical.  *  capacitor,  *  slanted  resistance  wire,  and 
2  3 

pin  '  methods,  all  of  which  measure  free  surface  velocities  rather  than  stress  directly.  The  great 
advantages  of  quartz  transducers  are  their  ability  to  measure  stresses  directly,  their  high  time  resolu¬ 
tion.  and  the  relatively  little  effort  required  for  their  use. 

Illustrated  schematically  in  Figure  2.3  are  two  experimental  arrangements  for  observing  stress- 
time  profiles  in  rock  specimens.  An  explosively  generated  plane  shock  wave  is  transmitted  into  the 
specimen  and  then  into  the  positively  oriented^  X-cut  quarts  crystal. 


Figure  2.3  Schematic  of  Wave  Profile  Experiments 


*  Crystals  oriented  so  that  the  stress  wave  is  made  to  propagate  along  the  crystallogrtgihlc  X-direction. 

^  The  crystal  is  positively  oriented  when  a  con^reasive  stress  produces  on  the  face  farthest  from  the 
specimen  a  voltage  wliich  is  positive  with  respect  to  that  on  the  other  face. 


For  small  values  of  the  resistor,  R.  the  circuit  current,  i,  is  given^’  by 


/•  B  a<T  {x,t) 

w  4c„  ax  ~ 

F 
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where  A  is  the  area  and  w  =  (Xg  -  x^,),  the  thickness  of  the  quartz  crystal.  The  sonic  velocity  in 
quartz  is  C^;  d  =  ^jj/c^^  is  the  one- dimensional  strain  piezoelectric  coefficient;  CTQ(Xp,t)  is  the 
time  variation  of  the  stress  wave  in  the  crystal  at  the  front  surface,  x„,  next  to  the  shocked  specimen; 
and  OQfXg,  t)  is  the  time  variation  of  the  wave  leaving  the  back  face,  Xg  .  Note  that  for  t  <w/Cq  , 
before  the  stress  wave  front  reaches  the  back  surface,  a_(x„,  t)  =  0,  so  that  time  variations  in 
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current  give  a  direct  measure  of  the  pressure- time  profile  in  quartz  at  the  quartz- specimen  interface. 

For  times  t  >  w/C„  the  current  is  proportional  to  the  stress  difference  across  the  crystal. 

W 

To  determine  the  amplitude  of  the  stress  history  in  the  shocked  specimen  from  the  quartz  stress 
history  at  the  interface,  the  degree  of  mismatch  at  the  interface  must  be  examined. 


Consider  the  one-dimensional  wave  system  at  a  time  just  after  transmission  of  the  wave  into 


quartz. 


CT  =<  0 

oQ 


u  _  =  0 
oQ 


Interface 


By  conservation  of  momentum  across  the  transmitted  wave  in  quartz,  we  have 
Q  oQ  Q  Q  Q  Q 

where  and  u,_  are  stress  and  particle  velocity,  respectively,  behind  the  transmitted  wave, 

W  vj 

is  the  quartz  density  ahead  of  the  wave,  is  the  wave  velocity  in  quartz,  and  Zq  is  the  quartz  im¬ 
pedance.  Similarly,  across  the  incident  wave  in  the  sample  we  had  before  the  wave  reached  the  interface 
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Aci’oss  the  reflected  wave,  momentum  consexwatlon  gives 


'’2  ■  "l  °  "  "2>  ' 

where  is  the  reflected  wave  impedance  of  the  sample  relative  to  the  flow  ahead  of  the  wave. 
Dividing  Equation  2.4  by  Equation  2.3  and  making  use  of  the  interface  boundary  conditions,  Ug  =  Uq  , 


(2.4) 


Og  =  <Tq  results  in 
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(2.5) 
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Since  —  =  — ^  ,  Equation  2.5  reduces  to 

"l  *^1 


u. 


1  +  r 


Ir 


Q  . 


1  + 


7 

"Q 


Thus  the  stress  <jj(t)  behind  the  incident  wave  in  the  specimen  is  given  in  terms  of  the  quartz 
interface  stress,  a,~(t),  and  an  impedance  factor  as 
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aj(t) 
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(2.6) 


provided  the  transmitted  wave  in  quartz  has  not  reached  the  back  face,  Xg 

Having  limited  ourselves  to  stresses  in  quat^z  less  than  the  elastic  limit,  the  crystal  impedance 
is  simply  Z^  =  p  ,-,C_  ,  where  C_,  is  the  appropriate  elastic  wave  velocity. 

To  completely  determine  the  impedance  factor  of  Equation  2.6,  it  is  necessary  to  make  the  approxi¬ 
mation  that  Z.  =  Z  =  p  ,C,  ,  where  C.  is  the  elastic  wave  velocity  in  the  specimen.  The  close 
Ir  I  01  1  X 

agreement  of  stress  time  profiles  obtained  by  quartz  transducers  with  those  obtained  by  other  means 
mentioned  above  indicates  that  the  impedance  approximation  is  sufficiently  accurate. 

The  use  of  quartz  crystals  for  stress  transducers  functioning  as  described  and  according  to 
Equation  2.1  is  novel  and  different  from  the  usual  application  of  quartz  crystals  to  measure  stresses. 

In  the  common  application,  dimensions  of  the  crystal  transducer  are  small  compared  to  a  stress  wave¬ 
length;  the  wave  transit  time  in  quartz  is  much  less  than  the  time  variations  in  the  stress  wave,  and  the 
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total  charge  released  by  the  crystal  is  proportional  to  the  average  stress  imposed  on  it.  In  the  Neilson 
application,  wave  transit  time  in  quartz  is  much  greater  than  time  variations  in  the  applied  stress,  waves 
of  length  less  than  the  transducer  dimensions  may  be  measured,  and  the  specimen-quartz-interface  stress 
is  projjortional  to  the  quartz  current  output. 

In  all  our  experiments  X-cut  crystals  of  synthetic  alpha  quartz*  were  used.  The  crystals  were 
positively  oriented^  and  generally  were  1/2  inch  in  diameter  and  1/8  inch  thick.  Crystals  of  this  size 
were  calibrated  by  R.  A.  Grah»,  ,n  of  Sandia  Corporation  with  step-function  type  wave  profiles  at  two  pres¬ 
sures,  14.5  and  23.7  kilobars,  resulting  from  high  velocity  projectile  impacts.  At  these  pressures  the 
piezoelectric  constant,  d  (average  value  =  2.05  x  10  coulombs/cm  -kilobars),  of  Equation  2.1  was 
found  to  vary  less  than  5  percent  for  the  initial  step-function  rise  of  the  input  wave.  At  a  time  when  the 
wave  has  progressed  halfway  through  the  crystal  the  value  of  d  has  increased  by  about  10  percent.  When 
the  wave  has  reached  the  back  surface  of  the  quartz  crystal  the  value  of  d  is  on  the  average  40  percent 
higher  than  it  was  initially,  but  at  this  time  the  stress-profile  measurement  is  completed.  That  d  is  not 
strictly  constant  is  a  consequence  of  both  mechanical  and  electrical  edge  effects  in  the  crystal.  In  our 
survey  of  stress  profiles  from  rock  materials  we  have  neglected  the  variation  of  d  with  distance  of 
wave  travel  into  the  crystal  and  have  used  a  constant  value,  d  =  2  x  10  coulombs/cm  -kilobar. 

The  primary  purpose  for  using  the  configuration  of  Figure  2.3(b)  was  to  increase  the  separation 
between  the  elastic  precursor  and  the  main  shock  wave  entering  the  transducer  crystal  from  the 
shocked  rock  specimen.  For  some  of  the  rock  materials  studied  it  was  found  that  the  elastic  precursor 
and  main  shock  wave  had  velocities  which  were  not  greatly  different.  As  a  result  the  main  shock  wave 
entered  the  transducer  crystal  immediately  after  entry  of  the  precursor  and  made  examination  of  the 
precursor  amplitude  and  profile  difficult  if  not  impossible. 

The  time-distance  diagram  of  Figure  2.4  illustrates  the  difficulty  encountered.  For  a  given  driving 
pressure  at  the  rock- explosive  interface,  two  waves  are  generated  in  the  rock  specimen,  the  elastic 
precursor  wave  of  nearly  constant  velocity,  Cg,  and  a  shock  wave  of  velocity  Vg,  whose  value  depends 
on  the  driving  pressure.  At  the  rock-quartz  interface,  the  elastic  precursor  first  enters  the  transducer 
at  time  t^  and,  after  adjusting  its  eimplitude  according  to  the  degree  of  impedance  mismatch  (Equation  2.6) 


*ObtaIned  from  the  Valpey  Crystal  Corp. ,  1244  Highland  Street,  Holliston,  Mass. 

+  Negatively  oriented  crystals  may  also  be  used  but  only  over  a  limited  range  of  stresses.^® 
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between  rock  and  quartz,  proceeds  through  quartz  with  the  elastic  velocity  of  quartz,  C,_,  emerging  at 
time  tg.  However,  before  time  tg,  the  rock  shock  wave  enters  quartz  at  time  tg,  transmitting  a  shock 
wave  into  the  transducer. 


TIME 


Figure  2.4  Time- Distance  Diagram  of  Double-Wave  System  in  Rock  and  Quartz  Transducer 

If  the  transmitted  wave  of  velocity  V_  has  a  sufficiently  large  amplitude,  impurity  atoms  of  .he 
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quartz  lattice  in  the  region  behind  this  wave  are  ionized,  as  suggested  by  Neilson  et  al.  As  a  result 

of  ionization  the  electric  field  breaks  down  and  conduction  of  electrons  occurs  behind  the  shock  wave. 
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Thus,  after  time  tg,  the  quartz  no  longer  operates  according  to  Equation  2.1  but  by  other  means. 

In  order  to  obtain  information  on  the  elastic  precursor,  we  must  try  to  satisfy  (tg  -  tj)  >  (tg  -  t^), 
which  gives  us  a  condition  on  the  ratio  of  specimen-to-transducer  thicknesses,  S/w  >  ^s^s^^Q^^S  ' 
Since  is  fixed  and  the  quarts  thickness,  w,  cannot  be  greatly  varied,  specimen  thickness,  S,  is 
governed  by  the  magnitudes  of  precursor  velocity,  Cg,  and  shock  velocity,  Vg,  in  the  specimen. 
Another  condition  of  the  thickness,  S,  is  that  the  thickness-to-diameter  ratio  of  the  specimen  be  suffi¬ 
ciently  small  so  that  a  true  plane  wave  is  maintained  in  the  specimen  as  it  progresses  into  the  quartz 
transducer. 

At  some  shock  pressure  (overdriving  pressure)  in  the  specimen,  which  is  determined  by  the  shape 
of  the  specimen  Hugoniot,  the  clastic-plastic  two-wave  structure  is  no  longer  generated,  the  elastic 
wave  is  said  to  be  overdriven  by  the  larger  amplitude  shock  or  plastic  wave,  and  at  this  pressure  Vg  is 
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equal  to  Cg.  At  shock  pressures  somewhat  below  this  level.  Cg  >  Vg  but  the  difference  (Cg  -  Vg)  is 
small  and  the  condition  on  S/w  cannot  be  met.  In  the  range  of  pressure  from  the  Hugoniot  elastic  limit 
pressure  to  the  overdriving  pressure.  Hugoniots  of  rock  materials  under  test  are  either  not  well  defined 
or  not  known  at  all.  As  a  result,  it  is  not  possible,  prior  to  making  a  wave  profile  measurement,  to 
accurately  estimate  magnitudes  of  Vg  which  may  be  achieved  by  different  driving  explosives  so  that  the 
condition  on  £>/w  can  be  realized. 

By  making  the  quartz  crystal  shown  in  Figure  2.4  a  stress  transmitter  rather  than  the  transducer 
and  placing  another  quartz  crystal  as  transducer  behind  it.  the  difficulty  of  an  early  enti^  of  the  shock  or  ,) 
plastic  wave  is  alleviated.  In  this  configuration,  the  first  and  second  wave  arrivals  are  separated  in 
time  by  an  amount  (t  -  t  )  which  is  greater  than  the  separation  time  (t_  -  t, )  when  no  stress  transmitter 
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is  used. 

In  Figure  2.5  are  shown  oscUlograms  of  the  elastic  precursor  in  Armco  iron,  which  was  used  to  test 
the  two  transducer  assemblies  of  Figure  2.3.  The  amplitude  of  the  elastic  wave  just  after  the  initial  rise 
is  about  6  kUobars.  Sweep  speed  for  all  traces  in  Figure  2.7  was  0.2  psec/cm,  except  for  the  top  trace 
in  (a),  where  it  was  0.1  psec/cm.  All  the  quartz  crystals  used  to  obtain  the  records  shown  had  diameters 
of  1/2  inch  and  thicknesses  of  1/8  inch.  TNT  was  the  driving  explosive  in  each  experiment. 

The  record  traces  of  Figure  2.5(a)  were  obtained  from  quartz  crystals  in  direct  contact  with  a 

1/2-inch-thick  piece  of  Armco  iron,  which  served  as  one  electrode  for  the  transducer.  The  symbol  A-Q 

is  used  to  note  this  type  of  assembly.  In  this  experiment  the  plastic  wave  from  Armco  entered  the  front 

surface  of  the  quartz  transducer  after  entry  of  the  precursor  but  before  the  precursor  arrived  at  the 

transducer  back  surface.  Thus  one  sees  the  6-kilobar  precursor  followed  in  time  by  the  large- amplitude 

(>  lOO-kUobar)  plastic  wave  which  takes  the  trace  off  scale.  These  records  are  similar  to  those 
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observed  by  others  '  '  in  Armco  iron. 

Figure  2.5(b)  compares  the  signals  from  one  experiment  with  two  quartz  transducers  on  a  1-inch- 
thick  piece  of  Armco.  The  top  trace  of  (b)  is  from  a  transducer  in  direct  contact  with  Armco  (A-Q). 
and  the  lower  trace  is  from  a  transducer  with  a  0.2-mil  aluminum  foil  placed  between  the  Armco  and 
quartz  (A-Al-Q).  as  shown  in  Figure  2.3(a).  The  traces  are  essentially  identical,  indicating  that  the 
0.2-mll  aluminum  foil  is  not  a  serious  perturbation  to  the  wave  from  Armco.  In  this  experiment  the 
plastic  wave  from  Armco  docs  not  arrive  at  the  front  surface  of  the  quartz  until  after  the  precursor 
transmitted  into  quartz  has  traversed  the  thickness  of  quartz.  The  time  at  which  the  precursor  wave  in 
quartz  reaches  the  quartz  back  surface  is  noted  by  an  abrupt  swing  of  the  trace  downward,  owing  to  the 


GC 


ELECTRODE 

ASSEMBLY 


RELATIVE 

GAIN 


SWEEP  ARMCO 

THICKNESS 


0.1 


0.2 


1/2 


0.2 


I 


0.2 


0.2 

0.2 


I 


Figure  2.5  Quartz  Transducer  Oscillograms  from  Armco  Iron 

With  Three  Types  of  Electrode  Assemblies: 

(a)  quartz  directly  on  Armco  iron, 

(b)  quartz  with  aluminum  foil  on  Armco, 

(c)  quartz  on  quartz  on  Armco  with  gold  electrode. 


term  of  Equation  2.1.  The  time  difference  between  the  start  of  the  trace  and  the  time  at  which 

it  turns  downward  is,  within  experimental  error,  the  calculated  elastic  wave  transit  time, in  quartz. 

Records  of  Figure  2.5(c)  are  from  an  experiment  identical  to  that  described  for  Figure  2, 5(b) 
except  that  the  lower  trace  of  (c)  is  from  a  transducer  in  contact  with  a  quartz  blank  (A-Q-Q),  as  shown 
in  Figure  2.3(b).  This  transducer  arrangement  is  also  shown  In  Figures  2.6  and  2.7.  The  traces  of 
oscillograms  (b)  and  (c)  indicate  that  the  Armco  elastic  precursor  is  very  well  reproduced  by  quartz 
transducers  in  either  assembly  configuration  of  Figure  2.3.  From  these  experiments  with  Armco^iron 
it  was  concluded  that  the  quartz  transducer  technique  would  be  adequate  to  examine  strees-clme  profiliS, 
in  rock  materials. 
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Figure  2.6  Assembly  Drawing  of  Quartz  Crystal  Transducer 
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Figure  2. 7  Photo  of  Quartz  Stress  Transducer 


Experimental  Results 

Rock  samples  examined  with  the  quartz  crystal  transducer  were  obtained  from  drill  cores,  cut  to 
various  thicknesses  and  ground  smooth  and  flat.  The  diameter-to-thickness  ratio  of  samples  was 
always  in  excess  of  2.5,  Aluminiun  (2024)  or  steel  (4340)  was  the  driver  plate  placed  between  sample 
and  explosive  as  shown  in  Figure  2.3(b).  Explosives  used  were  TNT,  Baritol,  and  Nitroguanidine. 

The  quartz  transducer  employed  for  all  experiments  described  here  was  that  of  Figures  2.3(b),  2.6,  or 
2.7.  Current  output  of  transducers  was  measured  on  Tektronix  555  and  581  oscilloscopes  as  the  voltage 
drop  across  a  50-ohm  resistor.  Timing  markers  and  voltage  calibration  signals  were  also  recorded 
along  with  signal  traces.  In  many,  but  not  all,  experiments  velocity  of  the  first  wave  was  measured 
by  recording  wave  transit  time  through  the  specimen. 

Granite 

Granite  rock  samples  utilized  in  experiments  reported  here  were  obtained  from  a  location  near 
Fallon,  Nevada  and  are  referred  to  as  Shoeil  granite.  Typical  wave  profiles  recorded  by  the  quartz 
transducer  on  two  thicknesses  of  Shoal  granite  are  shown  in  Figure  2.8.  Rise  times  of  waves  to  the 
first  peak  are  noted  by  arrows  labeled  t^.  Time  of  arrival  of  the  second  wave  is  noted  by  arrows  labeled 
t  =  (t .  -  t  ).  The  scale  to  the  right  of  the  figure  indicates  the  magnitude  of  stresses  in  the  granite  wave 

S  4  o 
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and  is  derived  from  the  voltage  scale  to  the  left  using  d  =  2  x  10  coulombs/cm  -kilobar  in  Equation  2.1 

and  Z  =  1  in  Equation  2.6.  Between  the  times  t^  and  there  appears  to  be  a  relaxation  of  peak  stress 

in  the  first  wave.  The  stress  profiles  of  Figure  2.8  prior  to  the  time  t  are  believed  to  be  elastic 

s 

precursor  waves  in  granite. 
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TIME  ( M  SEC) 

Figure  2.8  Stress  Wave  Profiles  Observed  in  Two  Thicknesses  of  Granite 
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Although  the  dynamic  yield  stress  of  granite  is  expected  to  be  quite  high  (Grine  has  reported 

86  kilobars  for  another  type  of  granite),  it  cannot  be  concluder  from  the  records  of  Figure  2.8  that 

Shoal  granite  has  an  elastic  limit  of  60  kilobars.  There  are  two  reasons  for  this.  First,  the  dynamic 

piezoelectric  coefficient,  d,  for  quartz  is  not  veil  known  at  stress  levels  near  or  above  50  kilobars, 

5  6 

and  second,  at  about  30  kilobars  quartz  itself  begins  to  yield  '  and  is  no  lunger  an  adequate  transducer 


Calibration  of  quartz  crystals  with  diameter-to-thickness  ratio  of  5  by  projectile  impact  experi- 

21  "“8 
ments  has  shown  that  d  increases  fairly  linearly  with  increasing  stress,  being  2.05  x  10  and 

*8  2 

2.48  X  10  coulombs/cm  -kilobar  at  25  and  45  kilobars.  respectively. 

From  this  information  we  may  surmise  that  the  peak  amplitudes  of  waves  in  Figure  2.8  are  some¬ 
what  more  than  40  kilobars  and  probably  less  than  SO  kilobars.  Because  the  precursor  amplitude  in 
granite  is  higher  than  was  expected,  the  wave  profiles  observed  were  obtained  with  the  quartz  transducer 
operating  in  a  region  where  its  piezoelectric  prt^erties  are  not  well  known.  Consequently,  the  wave  pro¬ 
files  obtained  for  granite  (except  for  experiment  number  188)  must  be  regarded  as  being  of  a  qualitative 
nature. 

Table  2.1  is  a  list  of  experiments  performed  with  various  thicknesses  of  Shoal  granite  showing  the 
explosives  used  and  the  rise  times  and  second  wave  arrivals  observed. 

TABLE  2.1 


Shoal  Granite  Experiments 
(1  cm  thick  Aluminum  Driver) 


Sample 

Rise  time 

Second  wave 

Experiment 

Explosive 

thickness 

t 

arrival  time 

No. 

type 

Sglcm) 

r 

(fisec) 

tg  (psec) 

215-1* 

TNT 

0.340 

~0.1** 

~0.24 

215-2* 

TNT 

0.986 

0.14 

0.40 

189 

TNT 

2.480 

0.34 

0.47 

197 

TNT 

4.032 

0.71 

0.80 

216 

Baritol 

0.340 

0.06 

0.22 

208-1+ 

Baritol 

0.582 

0.08 

0.27 

204 

Baritol 

0.824 

0.16 

0.28 

203 

Baritol 

1.062 

0.16 

0.38 

208-2+ 

Baritol 

1.320 

0.28 

0.42 

192 

Baritol 

2.415 

-0.55 

-- 

188 

Nitroguanidine 

2.615 

"O.O 

-- 

*  Both  granite  samples  on  one  TNT  experiment 
"'''^Doubtful  result 

t  Both  granite  samples  on  one  Baritol  experiment 


hi  Figure  2,9  are  plotted  granite  wave  profiles  obtained  from  experiments  with  TNT.  Since  from 
Table  2.1  rise  time,  t^,  appears  to  increase  linearly  with  increasing  sample  thickness  we  have  plotted 
as  abscissa  the  dimensionless  time,  r  =  tC^/S^  ,  using  the  average  value  of  first  wave  velocity  in 
granite,  C^  =  0.57  ±  0.03  cm/jusec,  obtained  from  independent  measurements.  Small  arrows  in  the 

figure  indicate  the  times,  t  ,  considered  to  be  second  wave  arrivals.  The  pressure  scale  on  the  right  of 
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the  figure  is  provided  to  show  only  the  magnitude  of  stresses  involved  and,  just  as  in  Figure  2.8,  over¬ 
estimates  stresses  particularly  at  higher  stresses  and  later  times.  From  Figure  2.9  it  is  seen  that  the 
reduced  rise  time  is  roughly  constant  for  all  thicknesses  of  granite  except  for  the  thinnest  sample, 

Sq  =  0.34  cm.  The  wave  profile  for  this  sample  seems  to  be  anomalous.  The  reason  for  the  anomaly  is 
unknown;  it  may  be  that  the  wave  in  the  sample  was  not  plane  or  was  tilted  or  possibly  that  the  sample 
was  not  uniform  and  homogeneous.  That  the  reduced  rise  time  is  nearly  constant  for  all  the  other  sample 
thicknesses  suggests  that  granite  exhibits  dissipative  or  viscous  effects. 

Granite  wave  profiles  produced  by  Baritol  explosives  are  shown  plotted  in  Figure  2.10.  Here  again 
the  reduced  rise  time  for  all  thicknesses  of  granite  is  roughly  constant.  Note  that  the  reduced  rise  time 
for  the  sample  of  0.34  cm  thickness  is  about  the  same  as  for  the  other  thicker  samples  and  does  not 
appear  to  be  anomalous  as  was  the  case  for  the  0.34  cm  thick  sample  driven  by  TNT. 


Figure  2.9  Granite  Stress  Wave  Profiles  as  a  Function  of  Reduced  Time  for  TNT  Kxpcriments 


Figure  2,10  Granite  Stress  Wave  Profiles  as  a  Function  of  Reduced  Time  for  Baritol  Experiments 
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The  influence  of  driving  pressure  on  the  rise  time  of  waves  in  granite  is  seen  in  Figure  2.11  where 
profiles  initiated  by  TNT,  Baritol,  and  Nitroguanidine  in  samples  about  2,5  cm  thick  are  shown.  Driving 
pressures  produced  in  granite  by  TNT,  Baritol  and  Nitroguanidine  are  about  200,  130  and  40  kilobars, 
respectively.  Evidently  the  rise  time  decreases  as  driving  pressure  increases.  The  wave  profile  from 
the  Nitroguanidine  experiment  is  particularly  interesting.  Peak  amplitude  of  the  wave  is  about  35  kilobars, 
and  this  value  may  be  considered  as  having  some  reliability  since  quartz  transducers  function  properly  at 
these  stress  levels.  There  is  no  evidence  of  a  second  wave  arrivt.^  before  one  transit  time  of  the  wave  in 
quartz,  from  which  we  may  conclude  that  the  Hugoniot  elastic  limit  in  Shoal  granite  is  in  excess  of  35  kilo¬ 
bars.  If  the  elastic  limit  were  less  than  about  35  kilobars  a  double  wave  structure  would  certainly  have 
been  observed.  The  Nitroguanidine  plane  wave  generator  used  in  this  experiment  was  fabricated  by  W. 

3 

Benedick  of  Sandia  Corporation.  Density  of  the  explosive  was  0.41  gm/cm  .  Shock  pressures  generated 

22 

in  aluminum  by  Nitrogueinidine  with  this  density  are  found  to  be  between  40  and  50  kilobars. 
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Figure  2.11  Effect  of  Driving  Pressure  on  Kise  Time  of  Waves  in  Granite 

Because  the  first  wave  stresses  in  granite  were  so  high  the  quartz  crystal  transducers  found 
themselves  operating  at  stress  levels  where  their  piezoelectric  properties  are  not  known  and  at  stress 
levels  near  which  quartz  itself  yields.  Consequently  we  do  not  have  great  confidence  in  the  times  denoted 
for  second  wave  arrivals  in  Figures  2.9  and  2,10,  There  are,  however,  other  indications  that  a  double 


wave  structure  is  exhibited  by  Shoal  granite.  First  wave  velocities  appear  to  be  nearly  constant, 

Cq  =  0,57  ±  0,03  cm/psec,  regardless  of  sample  thickness  and  driving  explosive.  If  a  single  shock  wave 

were  stable  in  granite  the  more  energetic  explosives  would  produce  waves  of  higher  velocity;  this  was 

not  observed.  Also,  if  the  first  wave  velocity  is  constant  the  second  or  shock  wave  velocity,  ,  in 

o 

granite  may  be  calculated  from 


and  compared  with  other  Hugoniot  data  for  Shoal  granite.  From,  the  values  of  t  and  S.  listed  in 

S  Cl 

Table  2.1  one  obtains  the  average  values  of  second  wave  velocity,  =  0.500  cm/psec,  from  TNT 
experiments  and  =  0.493  cm/psec  from  Baritol  experiments.  Knowledge  of  the  shock  pressure 
generated  in  the  aluminum  plate  by  TNT  and  Baritol  enables  one  to  determine  the  pressure  and  particle 
velocity  associated  with  the  second  wave  in  granite.  By  the  steuidard  impedance-mismatch  technique 
(e.  g. ,  see  Figure  1  of  Reference  6)  one  obtains  the  following  Hugoniot  data  for  Shoal  granite: 

P„  =  163  Kbar,  u„  =  0.118  cm/psec,  y_  =  0.500  from  TNT 
■  Pg  =  133  Kbar,  Ug  =  0,096  cm/psec,  =  0.493  from  Baritol. 

To  obtain  the  above  figures  it  was  necessary  to  estimate  the  Hugoniot  elastic  limit,  Pj,  of  Shoal 

granite.  This  Was  taken  as  50  kilobars.  Other  assumed  values  of  Pj  between  40  and  60  kilobars  change 

the  values  of  listed  above  by  less  than  10  percent.  The  two  Hugoniot  data  obtained  from  the  quartz 
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crystal  measurements  are  shown  along  with  other  Hugoniot  data  in  Figure  2.12.  That  the  quartz 
crystal  data  are  consistent  with  Hugoniot  data  obtained  by  other  means  suggests  to  us  that  a  double  wave 
structure  is  actually  being  observed  by  the  quartz  transducers. 

Results  of  wave  profile  experiments  with  quartz  crystal  transducers  strongly  point  to  the  existence 
of  a  two-wavt  structure  in  Shoal  granite  at  driving  pressures  below  about  300  kilobars.  The  experi¬ 
mental  records  do  not  demonstrate  this  positively,  however.  If  a  double  wave  structure  exists,  as 
appears  to  be  the  case,  the  granite  precursor  amplitude  or  dynamic  yield  strength  is  felt  to  be  in  excess 
of  35  kilobars  and  probably  between  40  and  SO  kilobars. 

In  the  wave  profiles  measured  there  appears  to  be  a  relaxation  of  stress  after  the  first  peak  or 
yield  stress  in  granite  is  achieved.  As  noted  before,  this  relaxation  may  be  Instead  some  idiosyncrasy 
of  the  quartz  transducer.  Further  calibration  experiments  on  quartz  transducers  are  required  with  50- 
kilobar  waves  and  with  various  wave  profiles  or  strain  rates  before  the  relaxation  can  be  definitely 
assigned  to  the  granite  waves. 
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Figure  2.12  Quartz  Crystal  Data  Plotted  With  Hugoniot  Data  for  Granite 

A  more  definite  result  of  the  granite  experiments  is  that  the  rise  time  of  the  first  (and  possibly  only) 
wave  increases  with  increasing  distance  of  wave  travel  in  granite.  This  observation  is  considered  as 
evidence  for  viscous  behavior  of  granite.  ■ 

It  must  be  emphasized  that  the  interpretations  given  the  records  obtained  from  Shoal  granite  with 
quartz  crystal  transducers  are  largely  qualitative  and  certainly  not  unique.  Before  more  definitive 
information  on  granite  may  be  obtained  further  experiments  are  required. 

Salt 

Samples  of  rock  salt  used  in  wave  profiles  experiments  were  obtained  from  a  location  near 
Winnfleld.  La. ,  and  are  referred  to  as  Cowboy  salt.  All  samples  were  highly  nonuniform  in  the  sense 
that  they  were  permeated  by  numerous  cleavage  and  shatter  cracks.  Density  of  Cowboy  salt  was 

3 

2.15  gm/cm  and  the  average  first  wave  velocity  measured  was  0.40  t  0.02  cm/psec,  slightly  less  than 
24 

the  reported  seismic  velocity.  0.438  cm/psec. 

The  crystal  transducer  of  Figure  2,6  was  used  In  all  experiments.  Because  a  very  low  precursor 

23 

amplitude  In  salt  was  expected,  of  the  order  of  1  or  2  kllobars,  explosives  In  direct  contact  with  salt 
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or  explosives  with  a  driver  plate  such  as  aluminum  could  not  be  used  as  it  was  feared  that  the  high  driving 
pressures  generated  would  overdrive  the  elastic  precursor.  The  driver  plate  chosen  for  the  salt  experi- 

4 

ments  was  4340  steel  tzom  which  waves  of  amplitude  16  to  25  kilobars  are  obtained,  depending  on  the 
Rockwell  C  hardness  of  the  steel.  These  waves  in  4340  are  actually  the  elastic  precursors  of  steel.  TNT 
explosive  was  used  to  initiate  the  waves  in  steel. 

In  Figures  2.13  and  2.14  are  shown  the  results  of  wave  profile  measurements  in  salt.  For  each 
experiment  a  quartz  transducer  was  placed  on  the  back  surface  of  the  4340  driver  plate  as  well  as  on  the 
salt  sample.  By  this  means  the  wave  profile  produced  in  steel  by  TNT  and  transmitted  into  salt  was 
measured.  The  topmost  profiles  of  Figures  2.13  and  2.14  are  the  waves  entering  the  salt  samples;  lower 
profiles  in  the  figures  are  those  from  salt. 

Salt  sample  thicknesses  for  Figures  2.13(a)  and  2.13(b)  were  2.54  and  2.01  cm,  respectively.  In  each 
of  these  experiments  the  4340  driver  plate  was  7.62  cm  thick.  The  quartz  transducer  used  to  obtain  the 
profiles  of  Figure  2.13(a)  was  1/4  inch  in  diameter  by  1/8  inch  thick,  while  that  used  to  measure  profiles 
of  Figure  2.13(b)  was  1/2  inch  in  diameter  by  1/4  inch  thick. 

In  Figures  2.14(a)  and  2.14(b)  the  4340  driver  plate  was  2.54  cm  thick,  and  the  salt  sample  thick¬ 
nesses  were  1.05  and  0.886  cm.  respectively.  Diameter  and  thickness  dimensions  of  quartz  transducers 
were  1/2  inch  by  1/8  Inch  for  both  experiments  of  Figures  2.14(a)  and  2.14(b). 

Impedance  factors,  Z.  of  Equation  2,6  used  to  reduce  records  were  2.02  for  steel  to  quartz  and 
0.784  for  salt  to  quartz. 

That  the  wave  profiles  measured  in  4340  steel  are  not  identical  in  Figures  2.13  and  2.14  is  due  to 
the  fact  that  the  steel  driver  plates  did  not  all  receive  the  same  heat  treatment,  and  consequently  they  are 

4 

expected  to  yield  different  profiles.  The  Rockwell  C  hardness  of  steel  drivers  in  Figure  2.13  is  not 
known.  For  the  steel  drivers  of  Figure  2.14  the  hardness  is  estimated  to  be  between  35  and  50. 

It  is  difficult  to  make  a  comparison  of  the  stress  profiles  from  salt  because  of  the  different  stress 
profiles  transmitted  into  it  and  because  of  its  different  thicknesses.  Nevertheless,  it  can  be  seen  that 
under  roughly  similar  conditions  (see  Figure  2.13)  the  salt  profile  amplitudes  at  a  given  time  may  differ 
by  as  much  as  a  factor  of  two.  Perhaps  these  differences  in  salt  profiles  are  a  result  of  the  nonuniform¬ 
ities  in  samples.  All  of  the  salt  stress  profiles  do.  however,  have  the  same  general  shape,  a  slowly 
rising  wave  form  without  prominent  maxima  or  peaks.  If  Cowboy  salt  had  a  Hugoniot  elastic  limit  at 
about  1  or  2  kilobars  and  a  stress-strain  curve  similar  to  that  of  Figure  2.1(b),  a  double  wave  structure 
should  have  been  observed  from  the  driving  pressure  of  4340  steel.  Since  the  wave  profUea  observed  in 


76 


(«) 


(b) 

Ficure  2.13  Stres*  Wave  ProfUee  in  Salt,  3- Inch  4340  Steal  Driver 


TIMK  !)••»» 


(b) 

Figure  3.14  atresi  Wave  ProfllM  in  Suit.  l>lRch  4340  Steel  Driver 


salt  are  more  nearly  similar  to  those  of  Figure  3.1(d)  it  is  felt  that  the  stress-strain  relation  of  Cowboy 
salt  has  a  region  td  anomalous  compressibility  for  stresses  in  salt  less  than  20  kilobars. 

Basalt 

From  Area  18  (Danny  Boy  site)  of  the  Nevada  test  site  a  sample  of  basalt  rock  was  obtained  for 
wave  profile  experiments.  A  stress  wave  profile  measured  in  basalt  is  shown  in  Figure  2.15.  A  plane 
wave  explosive  generator.  Baritol  explosive,  a  1  cm  thick  aluminum  driver  plate,  and  the  transducer 
assembly  of  Figure  2.6  were  used  to  obtain  the  profile  of  Figure  2.15.  Density  the  basalt  sample  was 

3 

2.67  gm/cm  and  measured  first  wave  velocity  was  0.524  cm/psec. 


Figure  3,15  Stress  Wave  Profile  la  Basalt 

The  peak  stress  occurring  at  0.07  psec  indicates  a  Hugonlot  elastic  limit  of  38  kilobars  for  basalt. 
This  precursor  amplitude  may  be  compared  with  the  50-ki)obar  Hugonlot  ela^ic  limit  observed  by 
Frkman  in  a  basalt  of  density  2.82  gm/em.  Between  the  time  of  peak  stress  and  the  time  of  setond 
wave  arrival  a  relaxation  in  stress  of  about  7  kilobars  occurs. 
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The  second  wave  arrival  time  is  clearly  noted.  From  Equation  2.7,  we  may  estimate  the  second 
wave  velocity,  as  was  done  for  granite,  and  calculate  a  point  on  the  basalt  Hugoniot.  Shock  velocity  is 
found  to  be  0.485  cm/psec.  Using  38  kilobars  as  the  Hugoniot  elastic  limit  and  the  appropriate  aluminum 
reflection  Hugoniot,  shock  pressure  and  particle  velocity  associated  with  the  second  wave  are  1 27  kilobars 
and  0.102  cm/psec.  This  datum  is  consistent  with  other  Hugoniot  data  for  basalt. 

Unlike  the  transducer  used  with  granite,  the  quartz  transducer  which  measured  the  basalt  profile 
of  Figure  2,15  operated  at  stress  levels  below  the  dynamic  yield  stress  of  , quartz  where  the  piezoelectric 
response  of  quartz  is  relatively  well  known.  As  a  result,  considerably  more  reliability  may  be  attached 
to  the  basalt  profile  than  to  the  granite  profile. 

Conclusions 

14 

The  quartz  crystal  technique  developed  by  Neilson  for  the  measurement  of  stress  wave  profiles 
in  metals  has  been  adapted  for  studies  of  stress  waves  in  rock  materials.  -This  new  technique  is  capable 
of  measuring  stresses  directly  with  high  time  resolution  at  stress  levels  below  about  45  kilobars. 

A  cursory  examination  has  been  made  of  three  rock  materials  to  investigate  the  properties  which 
characterize  nonfluid-nonelastic  stress  wave  propagation.  Results  of  experiments  indicate  that  Shoal 
granite  exhibits  a  viscous  nature  whose  effect  is  to  lengthen  the  rise  time  of  a  step-function  type  stress 
wave  as  it  propagates  in  granite.  Rise  times  of  waves  in  granite  appear  to  be  influenced  by  the  driving 
pressure.  The  Hugoniot  elastic  limit  of  Shoal  granite  is  in  excess  of  35  kilobars  and  probably  is 
between  40  and  50  kilobars. 

Cowboy  salt  appears  not  to  have  a  Hugoniot  elastic  limit  during  times  over  which  observations  were 
made;  rather  the  wave  profiles  observed  in  salt  suggest  that  the  stress-strain  relation  for  salt  is  such 
that  compressions  are  anomalous  at  stresses  below  about  20  kilobars. 

A  wave  profile  measured  in  basalt  established  the  Hugoniot  elastic  limit  as  nearly  40  kilobars  and 
also  indicated  that  relaxation  effects  exist  in  basalt. 

While  these  initial  experiments  have  not  provided  reliable  quantitative  information  on  the  rock 
materials  examined,  they  have  pointed  out  some  effects  which  will  require  consideration  in  the  formula¬ 
tion  of  models  of  nonfluid-nonelastic  behavior  of  geologic  solids.  Before  accurate  calculations  can  be 
made  of  nonfluid-nonclastic  stress  wave  propagation  from  buried  explosions,  the  results  of  laboratory 
plane  wave  experiments  must  first  be  understood. 
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Section  3 


LONG -PERIOD  DISPLACEMENT  GAGES 
J.  W.  Wistor,  J.  A.  Beyeler,  andG.  J.  Hansen 

Instrument  development  for  observation  of  long-period  earth  displacement  produced  by  under¬ 
ground  explosions  has  been  directed  toward  relative  displacement  devices  (long-base  strain  gages), 
inertial  displacement  gages,  and  electronic  single  and  double  integrators. 

Strctched-Wire  Gages 

The  spring-loaded  wire  gages  used  by  Sandia  Corporation  on  Plumbbob  Projects  1.5  and 
26.4b  are  typical  examples  of  state-of-the-art  stretched -wire  relative  displacement  gages.  The 
gage  system  comprises  a  spring-loaded  music  wire  stretched  under  continuous  tension  between  a 
fixed  anchor  and  a  spring-loaded  drum  in  the  transducer  unit.  An  auxiliary  helical  extension  spring 
is  used  to  preload  the  main  music  wire  to  ensure  optimum  frequency  response  from  the  system. 

The  music  wire  is  operated  within  a  flexible  metal  tube,  which  pcj-mits  free  movement  of  the  wire 
and  effectively  decouples  it  from  the  surrounding  medium  (Figure  3. 1).  This  instrument  had  a 
usable  displacement  range  of  about  +20  inches  and  a  response  of  20  milliseconds  per  foot  to  a  step 
displacement.  A  modified  form  of  this  gage  was  used  at  several  locations  on  the  Small  B03’  event 
at  NTS, 

The  Small  Boy  instrumentation  consisted  of  three  deep  drill  holes  which  contained  particle 
motion  gages  at  various  depths  between  the  surface  and  400  feet.  The  displacement  gages  used  in 
those  drill  holes  were  constructed  in  the  following  manner: 

• 

A  large  canister  at  the  bottom  end  of  each  instrument  array  housed  the  spring-loaded  dis¬ 
placement  pickoffs  and  the  preload  springs.  The  remote  displacement  zeroing  and  spring  tension¬ 
ing  mechanisms  were  housed  in  satellite  canisters  at  various  locations  in  the  uphole  portion  of  the 
array.  These  satellite  canisters,  which  also  contained  velocity  gages  and  accelerometers,  acted 
as  the  forward  anchor  points  for  the  displacement  gage  arrays  (see  Figures  3.2,  3.3,  3.4,  and 
3.5).  These  gages  could  measure  up  to  70  inches  of  relative  displacement  and  were  operated  over 
spans  of  from  100  feet  to  375  feet. 


Figure  3.3  Installation  of  Small  Boy  375-foot  Instrument  Array 
(Note  main  canister  and  one  satellite  canister.) 
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Figure  3.4  Small  Boy  Main  Canister  Detail 
Showing  Four  Displacement  Pickoffs 


* 


Figure  3.5  Small  Boy  Satellite  Canister  Detail  Showing 
Music  Wire  Takeup  Winch  and  Canister  Orienting  Unit 


In  retrospect  several  recommendations  can  be  made  with  respect  to  long-span  displacement 
gage  assemblies.  First,  any  long-base,  spring-loaded  music  wire  gage  system  should  be  limited 
to  about  a  200-foot  span.  For  longer  spans  frictional  forces  between  the  music  wire  and  flexible 
metal  decoupling  tubing  become  so  large  that,  even  with  additional  wire  preloading,  gage  response 
time  is  seriously  affected.  Second,  for  the  particular  type  of  relative  displacement  gage  under 
discussion  —  i.e. ,  one  in  which  the  displacement  pickoffs  and  music  wire  preloading  springs  are 
located  at  the  bottom  or  remote  end  of  the  instrument  array  —  overall  system  response  time  is  in¬ 
fluenced  by  several  additional  considerations.  In  this  system  the  ground  displacement  signal  ar¬ 
riving  at  the  forward  anchor  at  time  t^  is  not  sensed  by  the  displacement  pickoff  until  some  later 
time  (tg  +  At),  where  the  additional  delay  time  At  is  determined  by  the  transit  time  of  a  sonic  wave 
in  the  associated  length  of  music  wire  and  the  response  time  of  the  displacement  transducer-music 
wire  system.  In  this  system  the  frictional  and  inertial  forces  required  to  accelerate  the  music  wire 
assembly  degrade  the  system  response.  In  the  alternate  system  (displacement  transducer  and  pre¬ 
load  spring  forward),  the  response  time  of  the  displacement  pickoff  is  the  most  significant  part  of 
the  delay  time.  Fortunately  the  portion  of  the  delay  lime  associated  with  the  transit  time  of  the 
sonic  wave  can  be  calibrated  out  of  the  response  equation  with  reasonable  accuracy.  For  the  Small 
Boy  event  the  advantages  associated  with  remotely  locating  the  displacement  pickoff  were  felt  to 
outweigh  the  disadvantages,  particularly  because  of  the  peculiar  geometry  of  the  Small  Boy  instru¬ 
mentation  plan.  Primary  advantages  were  of  course  associated  with  the  additional  shielding  afforded 
by  remote  location.  Thus  instrument  disturbances  due  to  nuclear  and  electromagnetic  radiations 
and  to  ground  shock  were  minimized.  For  performance  details  consult  Reference  1. 

For  additional  information  on  the  Plumbbob  type  relative  displacement  gage  system,  .see  Refer¬ 
ences  2  and  3. 

Rigid- Rod  Gages 

A  second  type  of  relative  displacement  device  consists  of  a  rigid  rod,  one  end  of  which  is  se¬ 
curely  anchored  to  the  test  medium.  The  rod  is  enclosed  within  a  flexible  metal  housing  containing 
self-align,..g  bearings  at  fixed  intervals.  The  metal  housing  decouples  the  rod  from  the  test 
medium,  and  the  bearings  permit  free  rod  movement  and  damp  out  transverse  rod  vibrations,  A 
distance  transducer  is  used  to  measure  the  relative  displacement  between  the  free  rod  end  and  a 
second  fixed  anchor  (Figures  3.6  and  3.7).  Several  lO-foot  gages  in  various  configurations,  set 
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up  in  the  laboratory  to  simulate  the  buried  condition,  were  excited  by  an  electro<fynamic  shaker. 

First  significant  resonances  for  3/8-  and  1/ 2-inch -diameter  rods  occurred  in  the  region  of  300  to 
450  cycles  per  second.  Smaller  diameter  rods  resonate  laterally  at  frequencies  so  low  as  to  pre¬ 
clude  their  use.  Two  10-foot  gages  were  buried  near  each  of  two  5000-pound  HE  bursts  and  per-  > 

formed  satisfactorily  (Figure  3.8a).  Several  10-foot  gages  were  built  with  a  remote  zeroing  capa¬ 
bility  for  Gnome  and  NTS.  A  record  obtained  from  this  gage  at  a  distance  of  950  feet  from  the  Gnome  * 

explosion  is  shown  in  Figure  3.8b,  and  an  NTS  record  in  Figure  3.8c. 

Recently  the  original  10-foot  relative  displacement  gage  has  been  extensively  redesigned. 

These  modifications  have  produced  several  advantages.  Smaller  and  more  sensitive  displacement 
pickoffs  can  now  be  used.  This  reduces  the  active  gage  span  from  10  feet  to  3  feet  with  no  reduc¬ 
tion  in  gage  sensitivity.  Strains  between  the  limits  2  microinches  per  inch  tc  20,000  microinches 
per  inch  can  be  measured  with  this  device.  For  deep  hole  installations  the  remote  rod  unlatching 
and  transducer  zeroing  mechanisms  have  been  redesigned;  this  has  produced  a  shorter  and  smaller 
diameter  motor  and  transducer  housing,  consistent  with  installation  in  smaller  diameter  drill  holes. 

The  length  has  been  reduced  from  a  nominal  16  inches  to  10  inches  and  the  diameter  from  5  inches  to 
2.25  inches.  A  recent  test  of  this  device  near  a  buried  1000-pound  HE  detonation  at  the  Coyote  test 
field  yielded  good  displacement  data. 

Inertial  Displacement  Gages  ' 

In  order  to  measure  ground  motion  directly,  a  reference  point  must  bo  cho.scn  which  remains 
fixed  v/ith  respect  to  the  earth.  Since  such  a  reference  point  must,  in  general,  be  located  at  a  great 
distance  from  the  point  of  motion,  it  is  more  convenient  to  use  an  inertial  typo  of  displacement  gage. 

Furthermore,  when  studying  largo-amplitude  displacements,  a  gage  having  a  sensing-element  motion 
(with  respect  to  the  gage  case)  of  only  a  fraction  of  the  case  movement  is  desirable. 

Two  mechanical-type  inertial  displacement  gages  were  designed  and  constructed.  In  both  of 
these,  the  relative  motion  between  sensing  mass  and  case  is  reduced  by  coupling  the  linearly  moving 
sensing  mass  to  a  flywheel.  This  is  accomplished  with  a  rack-and-pinlon  gear  in  the  first  gage  and 
a  recirculating  ball  nut  in  the  second. 

The  moving  rack  in  the  first  instrument  is  fastened  to  a  mass  carried  on  a  ball  bushing  which 
rolls  on  a  splined  shaft.  The  flywheel-and-pinlon  shaft  is  mounted  on  precision  ball  bearings  (Fig¬ 
ure  3.9).  If  M  is  the  mass  of  the  ball  bushing  assembly,  R  is  the  pinion  gear  radius,  and  I  is  the 
moment  of  inertia  of  the  flywheel  assembly,  then  the  case  motion,  D,  is  given  by: 


Kelative  Displacei 
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i''iguro  3.9  Diaplaconicnt  Measuring  Mechanism  Assembly  -  Model  II 


(3.1) 


where  y  is  the  relative  motion  between  the  case  and  the  ball  bustling. 

The  output  signal  i.s  generated  from  a  rotary  diff»>rential  transformer. 

The  .second  instrument  t  onsists  of  a  helical'Splined  shaft  whicti  supports  a  rei  irculating  ball- 
screw  nut  of  the  type  used  on  .some  automotive  .steering  gears.  The  inductive  pickoff  im  ludes  a 
ferrite  ring  mounted  on  the  ball-.screw  nut  and  four  laminated,  soft  iron  armatures  attached  to  ttie 
rase.  The  armatures  are  tupe.-ed  and  tlieir  coils  wired  in  a  four-arm  bridge  configuration  (Figure 
3.10).  In  this  gage,  the  case  motion,  O,  is  given  by; 


D  -  (l  +  )  y,  (3.2) 

'  mp  • 

where  I  is  the  moment  of  inertia  of  the  ball-screw  nut  assembly,  m  is  the  mass,  p  is  the  pitch  of  the 
screw  (linear  advance  per  revolution),  and  y  is  the  linear  displacement  of  the  nut  with  respect  to  the 
case. 
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Kiguri*  3.  10  Splim*  and  UltH'k  Assembly 


A  .small  sit'd  was  constructi'd  for  test  and  evaluation  of  displai  ement  gages  (see  l■'igul■e  3.  11). 
The  sled  consists  of  a  carriage  approximately  I  foot  high  which  rides  on  two  Thompson  ball  bushiitgs 
guided  by  a  pair  of  parallel  steel  rods.  The  carriage  is  driven  by  a  quick-release  coil  spring  and  re¬ 
bounds  from  a  second  coil  spring  at  the  far  end  of  the  track.  The  sled  displacement  output  as  a  func¬ 
tion  of  time  is  essentially  a  triangular  pulse  having  an  amplitude  which  is  adjustable  from  zero  to 
50  inches.  The  carriage  normally  makes  several  round  trips,  the  number  depending  on  the  strength 
of  the  initial  impulse.  The  largest  drive  spring  can  impart  an  acceleration  of  well  over  100  g  to 
the  carriage  and  gage  assembly  at  full  compression. 
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Figure  3.11  Spring- Driven  Test  Sled 


Displacement-versus-time  curves  were  plotted  from  the  displacement  sled  test  data  of  both  the 
rack -and -pinion  and  helical-splincd  inertial  displacement  gages  (Figures  3. 12  and  3. 13). 

Maximum  deviation  of  the  rack-and-pinion  curve  from  the  true  displacement  was  less  than 
5  percent  of  the  peak  displacement  amplitude  for  the  first  1.5  seconds  on  one  run,  and  8  percent  for 
the  first  2.3  seconds  on  a  second  run.  These  runs  were  made  with  different  initial  driving  forces. 
As  the  sled  slows  to  a  stop,  the  deceleration  forces  decrease  to  a  point  where  they  are  comparable 
to  the  frictional  forces  in  the  gage,  and  the  gage  error  increases.  In  both  of  the  runs  mentioned 
above,  the  peak  displacement  amplitude  was  about  20  inches,  and  the  sled  came  to  rest  in  approxi¬ 
mately  2-1/2  seconds.  The  present  threshold  sensitivity  of  the  gage  is  0.04  g. 
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Figure  3.12  Sled  Test  of  Rack-and-Pinion  Gage  -  Model  II 


Figure  3. 13  Sled  Test  of  Helical -Splined  Gage 

The  internal  friction  of  the  helical-splined  gage  was  high  enough  to  produce  errors  as  large  as 
30  percent  of  full  scale  after  the  first  oscillation  of  the  sled.  For  this  reason,  together  with  the 
facts  that  its  physical  size  becomes  excessive  when  designed  for  large  displacements  and  its  cost 
is  much  greater  than  the  rack -and -pinion  device,  our  efforts  were  directed  toward  improvement  of 
the  rack-and-pinion  gage. 


The  prototype  rack-and-pinion  and  helical-splined  inertial  displacement  gages  were  exposed 
on  the  HE  test  mentioned  previously.  Because  of  the  small  motion  expected,  the  flywheel  was  removed 


Displacement  (inches) 


from  the  rack-and-pinion  gage.  Performance  of  the  helical-splined  gage  was  unsatisfaOtory  because 
of  excessive  friction.  The  rack-and-pinion  gage  indicated  a  peak  displacement  of  1.0  inch  and  a 
residual  permanent  offset  of  0.25  inch  (Figure  3. 14a).  A  precise  preshot  and  postshot  survey  of  a 
gage-mount  bolt  showed  a  permanent  displacement  of  0.28  inch.  Since  the  maximum  relative  error 
of  this  type  of  gage  usually  occurs  near  the  end  of  the  record,  and  the  recorded  peaks  coincided  with 
the  crossover  points  of  a  velocity  gage  located  nearby,  it  is  believed  "that  the  record  was  a  good  rep¬ 
resentation  of  actual  ground  motion.  This  gage  was  later  used  at  NTS  for  surface  motion  measure¬ 
ments. 
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(a)  5000  pounds  HE  surface  burst  -  radial  displaccinent 
50  feet  from  burst  at  depth  of  2  feet. 


(b)  Hardhat  event,  NTS.  radial  displacement  at  shot  depth, 
600-foot  slant  range. 


Figure  3.14  Typical  Records  From  Inertial  Displacement  Gage 


The  rack-and-pinion  gage  has  been  repackaged  and  modified  to  permit  installation  in  deep  drill 
holes  (both  vertical  and  horizontal).  The  new  external  configuration  is  a  6-inch-diameter  sphere, 
and  provision  has  been  made  for  a  remote  leveling  and  calibration  capability  after  the  case  has  been 
grouted  in  position  (see  Figures  3. 15  and  3. 16).  Twelve  were  built  for  use  at  Gnome  and  NTS.  The 
record  obtained  at  the  950-foot  station  at  Gnome  (Figure  3. 17a)  indicates  a  residual  displacement  to¬ 
ward  ground  zero  which  is  probably  not  true.  This  discrepancy  is  explained  by  reference  to  an  ac¬ 
celerometer  record  from  the  same  location  (Figure  3.17b).  It  is  seen  that  accelerations  are  small 
after  about  0. 15  second.  Frictional  forces  in  the  gage  are  therefore  becoming  more  significant  com¬ 
pared  to  acceleration  forces,  and  the  gage  indication  is  not  reliable  after  this  time.  This  character¬ 
istic  was  also  evident  on  the  sled  tests.  A  record  from  Haidhat  at  NTS  is  shown  in  Figure  3. 14b. 


Figure  3. 15  External  Configurution  of  Inertial  Displacement  Gage 
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Figure  3.  16  Disassembled  View  of  Inertial  Displacement  Gage 


A  cceleration  (g)  Displacement  (inches) 


(a)  Ground  displacement  record,  950  feet  from 
Gnome  event  at  shot  depth. 


(b)  Ground  acceleration  record,  960  feet  from 
Gnome  event  at  shot  depth. 

Figure  3,17  Displacement  and  Acceleration  Records  From  a  Nuclear  Explosion  in  Halite 
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Several  ideas  have  evolved  for  gage  designs  which  use  the  same  basic  operating  principle  but 
which  should  exhibit  less  friction  than  present  models. 

Velocity  Gages 

Recently  considerable  effort  has  been  directed  toward  a  redesign  of  the  SRI  type  velocity  gage. 
This  device  is  basically  a  highly  overdamped,  low  natural  frequency,  pendulum  accelerometer. 

This  new  gage,  in  conjunction  with  a  single-stage  electronic  integrator,  has  yielded  good  displace¬ 
ment  data  on  experiments  conducted  with  buried  HE  detonations  and  the  displacement  test  sled. 

This  modified  velocity  gage  has  been  specifically  designed  to  extend  the  range  of  particle  velocity 
measurements  into  the  300  ft/sec  region  and  to  withstand  shock  loads  in  excess  of  1200  g.  To  effect 
this  capability  a  complete  redesign  of  the  pendulum  suspension  system  and  pendulum  and  armature 
assembly  was  necessary.  Several  other  changes  were  made  in  order  to  improve  the  accuracy  and 
repeatability  of  these  gages.  A  system  for  filling  these  gages  with  DC  200  damping  fluid  under 
vacuum  was  designed.  This  eliminates  damping  ratio  changes  (with  consequent  amplitude  and  fre¬ 
quency  response  errors)  caused  by  air  admixture  in  the  damping  fluid.  To  accommodate  tempera¬ 
ture  expansion  of  the  damping  fluid  a  flexible  bellows  assembly  was  incorporated  in  the  base  of  the 
velocity  gage.  Figures  3.18  and  3. 19  show  various  views  of  the  present  design. 


Figure  3.18  Experimental  Velocity  Gage.  Type  DS-B(H) 


Figure  3. 19  Disassembled  Experimental  Velocity  Gage  Showing  Two  Case  Configurations 


During  the  recent  experiments  conducted  on  subsurface  1000-pound  HE  detonations,  these 
gages  successfully  measured  peak  particle  velocities  between  0. 1  ft/ sec  and  235  ft/ sec;  corre¬ 
sponding  peak  accelerations  ranged  from  0.5  g  to  greater  than  1500  g.  The  235  ft/ sec  peak  velocity 
was  measured  at  shot  depth,  7  feet  from  the  center  of  a  1000-pound  HE  sphere,  detonated  in  desert 
alluvium.  A  reproduction  of  this  velocity  gage  record  is  shown  in  Figure  3.  20. 

This  velocity  gage  was  later  recovered.  Extensive  postshot  tests  showed  no  gage  damage; 
however,  a  slight  loosening  of  the  E-coil  hold-down  screws  was  observed.  As  a  result  of  this  test 
the  E-coil  hold-down  assembly  has  been  redesigned.  The  hold-down  fixture  has  been  strengthened, 
and  epoxy  is  used  to  lock  the  E-coll  and  hold-down  screws  in  place.  Transient  gage  response  records 
derived  from  experiments  conducted  on  the  displacement  test  sled  are  shown  in  Figures  3.21  and 
3.22. 

The  experimental  velocity  gages  discussed  above  fall  into  two  categories.  One  type  is  basic¬ 
ally  a  low-range  model  using  a  brass  pendulum,  used  to  measure  particle  velocities  in  the  range 
0. 1  ft/sec  to  25  ft/sec.  Production  models  of  this  gage  are  referred  to  as  type  DS-B(H).  The  "B" 
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Figure  3.21  Experimental  Instrumentation  Responses  to  Impulse  From  Test  Sled 


Figure  3,22  Experimental  Instrumentation  Responses  to  Impulse  From  Test  Sled 


refers  to  the  brass  pendulum  material  and  the  "H"  refers  to  the  horizontal  mode  of  operation.  The 
second  typo  of  velocity  gage  is  similar  to  the  above  in  overall  design;  however,  the  pendulum  is 
constructed  of  aluminum  alloy.  This  model,  referred  to  as  type  DS-A(H),  Is  designed  for  operation 
in  high  g  regions  and  is  normally  used  to  measure  particle  velocities  in  the  range  of  10  ft/sec  to 
300  ft/ see.  The  experimental  velocity  gage  record  of  Figure  3. 20  was  obtained  with  a  gage  of  this 
type.  The  two  experimental  velocity  gages  referred  to  in  Figure  3.21  as  X,  and  U2^  identi¬ 
fied  as  follows:  X  is  an  early  prototype  of  the  DS-B(H)  model;  U2X  early  prototype  of  the 

DS-A(H)  model. 


Electronic  Integrators 

Because  of  the  problems  of  design,  calibration,  and  installation  of  velocity  and  displacement 
instrumentation,  we  have  long  desired  to  be  able  to  obtain  velocity  and  displacement  data  by  integra¬ 
tions  of  accelerometers.  Several  attempts  to  perform  integrations  using  operational  amplifiers 
with  capacitive  feedback  have  been  made  in  the  past.  However,  long-term  carrier  system  drift  and 
operational  amplifier  zero  drift  liave  caused  the  results  to  be  unsatisfactory.  This  section  reports 
progress  made  in  overcoming  these  difficulties  and  describes  the  development  of  long-term  stabilized 
single  and  double  integrators. 

Single  Integrato r s 

The  basic  circuit  for  performing  integration  with  operational  amplifiers  is  shown  in  Figure 
3.  23a.  The  current  through  R  is  i  =  (o-  -  e^)/R.  Since  an  ideal  amplifier  draws  no  input  current, 
thi.s  current  is  also  the  charging  current  for  the  capacitor,  C.  The  voltage  across  C  is 

V  =  e  -e  =  (1+A)c  . 

c  s  o  s 

Therefore  if  A»l,  V  can  be  written 
c. 

V  =  Ac  =  -e  =  i  /*  idt  =  f  (c.  -  e  )dt. 

c  s  o  C  J  lie  J  1  s 
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Since  e  *  -e  /A  and  A  is  very  large,  the  equation  becomes  e  »  -  i  f  e.dt,  which  shows  that 

8  0  O  KC-  w  A 

the  output  voltage  is  directly  proportional  to  the  integral  at  the  input  voltage.  In  general  an  opera- 

4  8 

tional  amplifier  circuit  can  be  shown  as  in  Figure  3.  23b.  If  A  is  very  large  (10  to  10  ),  then 

e„  0.  Therefore,  E  /E  -  Z,/Z. . 
s  O  1  I  1 

Using  Laplace  transform  methods  to  analyse  the  circuit  shown  in  Figure  3.23a  we  obtain 

Zf  «  ^  and  Zj  =  R  where  S  =  Js) 

E  Z,  ,  E, 

_S  =  .  s  .  _I_  or  E  *  -  — ^ 

Ej  Zj  RCS  o  RCS  * 

Taking  the  inverse  transform  we  obtain 


which  is  the  same  answer  obtained  previously.  Therefore  the  ratio  of  feedback  to  input  impedame 
will  be  used  throughout  this  report  to  determine  cireult  performance. 

The  circuit  of  Figure  3.33a  suffers  a  number  of  difficulties  in  its  praetiral  application.  A 
typical  instrumentation  system  consists  of  a  velocity  or  acceleration  gage  whose  response  is  ampli¬ 
fied  by  a  carrier  amplifier;  output  uf  the  carrier  amplifier  is  then  fed  to  the  integration  circuitry. 
The  output  voltage  of  the  carrier  amplifier  varies  slowly  (approximately  c  30  my)  with  lime.  and. 
when  applied  to  the  circuit  of  Figure  3.33a.  eauses  e^  to  vary  several  volts,  whleh  is  sufficient  to 
saturate  the  signal  recorder.  In  addition,  all  actual  operational  amplifiers  require  a  small  iX'  cur¬ 
rent  flowing  through  Z^  to  stabilise  them.  This  curre.tl  varies  from  amplifier  to  amplifier,  and  for 
the  amplifiers  being  used  (the  Astredata  HA -1 10)  is  about  33  to  30  micromicroamperes.  This 
current  continuously  charges  the  capacitor  C  and.  unless  bled  off,  allows  e^  le  again  rise  to  saiura- 
liiMt. 

To  solve  these  problems  the  clreitii  shown  in  Figure  3. 33e  was  inventignied.  With  this  ctr« 
cull  we  ohuin: 


l«T 


£ 

i 


2 

1 


Zj  '  *  RjCj(S+l/RjCj)  IS+l/RjCj)  • 


Let  T  o  R.C,  >  RqC.  >  stabilixation  time  constant, 
c  11  2  3 


Then 


E 


R,CJS+1/T 

12  c 


and  for  a  unit  step  function  input 


O  R,C, 


t  -t/T 
-7;-  e  c 


Therefore  it  can  be  seen  that,  for  I  less  than  0, 1  t  .  this  circuit  will  integrate  with  less  than  10 

o 

percent  error.  In  practice  it  is  possible  to  oNain  values  of  T^.  ‘  M  seconds  or  mere  with  good 
long'term  stability.  This  allows  accurate  integrations  wer  about  S  seca€Uls.  For  most  measure 
meitfs  the  data  is  obtained  in  less  than  I  second.  The  slngle'iniegratien  circuit  employed  in  our 
tests  was  that  of  Figure  3.33c, 

Oo^le  Integrators 

The  practical  problems  associated  with  double  •integration  circuits  are  much  more  severe 
than  those  of  single  integrators.  The  most  common  method  of  performing  double  integration  is  to 
use  two  single  integrators  in  series.  This  yields  the  following  form  for  a  step  input: 


e 

o 


•5h/T 


c  . 


A  syalem  of  ibis  type  will  work;  however,  it  has  several  major  problems.  The  fttrm  of  the  es* 
poneotlal  re9ulreo  that  T  for  the  two  Iniegratora  to  aerleo  be  lartre  m  large  as  T  for  a  alngle 
ioMfrMor  to  yleW  the  aamo  acearacy .  Then  too.  tbe  oolpMl  drift  ibe  flrai  Mdogrator  U  in  tbe 
paaa  band  of  tb*  serond  and  la  tbereforo  imegraied. 


Much  of  the  series  amplifier  drift  problem  can  be  overcome  1^  using  the  circuit  shown  in 


Figure  3. 23d  for  which 


Eq  2(S  ^  l/2RgC^) 

’^i  (S+  2/R^Cj) 


2  2 

Let  R,C  =  4R„C„.  Then  for  a  step  input  e  =  t  /R  C,C  .  which  is  the  form  for  true  double  inte- 

1122  O  112 

gration.  However,  the  circuit  of  Figure  3. 23d  is  very  unstable.  Some  means  of  compensation 

must  be  supplied  to  eliminate  the  instability .  The  simplest  method  is  to  apply  a  leakage  path  from 

output  to  summing  point  and  to  block  long-term  drift  in  e.  as  shown  in  Figure  3. 23e.  This  circuit 

produces  an  oscillatory  form  in  the  output  voltage  for  all  values  of  R^  other  than  R^  =  R2.  If  R^  is 

made  equal  to  R„,  the  circuit  is  stable  and  we  have 
2 


fo  ,  2S  r _ ^  _ 

®i  L(s  +  +  z/R^Cj)  (s  +  i/2RjCj^) 


^  ^  _ 2S _ 

1 

.  ^b 

r  s  1 

n,'c,C2 

L(S+  l/RjC/J 

.(S  +  l/Rg'^g^^J 

For  a  step  voltage  input  we  obtain 


e 

o 


S' 


C* 


V 

2 

2C„T 
2  c 


„-t/T 

e  '  c 


This  is  the  form  desired,  performing  short-term  double  integration  but  having  long-term  stability. 

0 

T^  "  R2^2  made  large  to  obtain  accurate  data,  R^  is  limited  to  <  2x10  ohms  for 

stability;  therefore  Cj  is  determined  when  T^  is  specified.  This  means  the  gain  of  the  circuit  is 
controlled  by  the  value  of  Cj^.  For  the  actual  circuit  Cj^  is  calculated  to  be  several  hundred 


109 


nsicrofarads.  Low-leakage  capacitors  at  this  mafnitude  become  physically  impractical,  and  there¬ 
fore  this  circuit  has  limited  application.  To  attempt  to  eliminate  these  problems  the  circuit  in 
Figure  3. 23f  was  investigated.  If  >>  and  >>  R^  we  obtaiii 


fo  ^  ^  ^  2S  [ _ 

*^i  “i  Rj  CjCg  I  (S  +  l/RgCj)^  <S^  +  2S/RjC  J  +  l/RjCjCj^>. 


2S 


(s  +  1/2R2C3) 


(S  +  I/R3C3)  (S  +  2/RjCj)  (S  +  l/2RjCj^) 


and  if  =  ^RgC^  and  2RjCj^  =  *^3^3  '  again  obtain  the  proper  form 


E 

0  _  2S 

1 

s 

S  1 

■'l 

[(S+  l/RjC/l 

’  wl 

(S  +  l/RgCg)^] 

In  this  case  we  find  that  if  R3  and  Cg  are  specified  and  is  limited  to  a  reasonable  value,  gain 
can  still  be  adjusted  by  varying  R^.  This  circuit  can  be  built  with  reasonable  components  and 
double  integrates  nicely;  however,  it  is  very  unstable  because  the  voltage  divider  formed  by  Rg 
and  Rg  severely  limits  the  steady  state  stabilization  current  flowing  from  output  to  summing  point. 
Stabilization  can  be  restored  by  placing  a  capacitor  in  series  with  to  block  the  direct  current 
path.  The  circuit  shown  in  Figure  3. 23g  is  then  obtained.  Making  the  assumptions  that 
and  Rg  >  ^  Rg  see  that 


Z 


i 


[(S  +  2/RjCj)  (S  +  l/2RjCj^)] 


2  [(S+ l/2RgC3)  (S+2/R3C2)| 
^  ^3  l(S  +  l/RgCg)^  (S  +  l/RjCj)] 


no 


".'^l  ■  “jS  *“■  Vi  ■  "  “I'^b  ""  Vj  •  Vs  '  ’'b- 

Then 


^o  .  2f  .  2 

s 

^'b 

S  1 

[(S+1/R3C3)^1 

c  3 

(S  +  1/T^)®J 

which  is  again  the  form  desired. 

This  circuit  has  been  constructed  and  operated  with  many  values  of  T^  and  gain.  With  values 

Q 

of  Rg  ^  2x10  it  will  hold  zero  to  +100  mv  for  long  periods  of  time.  It  can  be  constructed  with  com¬ 
ponents  of  reasonable  value  and,  depending  upon  the  gain  desired,  accurate  double  integrations 
over  several  seconds  can  be  obtained.  Results  of  some  recent  experiments  using  the  displace¬ 
ment  test  sled  are  shown  in  Figures  3.21  and  3.22. 

Sample  Determination  of  Components 

For  purposes  of  example  the  sled  run  shown  in  Figure  3.22  will  be  used.  In  this  case  it  was 
desired  to  (a)  single  integrate  acceleration  to  obtain  velocity,  (b)  single  integrate  velocity  to  obtain 
displacement,  and  (c)  double  integrate  acceleration  to  obtain  displacement. 

Acceleration  (A)  =  50  g 
Velocity  (U)  =:  10  ft/sec 
Displacement  (D)  =  3  feet 

2 

Using  the  relationships  U  =  At,  D  =  Ut,  and  D  =  (1/2)  At  ,  the  values  of  the  integrator  gain  con¬ 
stants  can  be  determined.  For  case  (a)  above  we  obtain 

RjCg  =  t  =  U/A  =  10/(50)  (32. 2)  =  6.21  x  10  ^  seconds  where  all  components 
are  designated  as  shown  in  Figure  3.23c. 

Next,  values  for  R^  and  are  chosen  as  large  as  possible  consistent  with  the  desired 

9  *  6 

stability.  For  the  example,  R^  =  2  x  10  ohms  and  =  18x10  farads. 


If  RjCj  =  UgCg  =  T^.  then 


=  R,R„C,C„  =  (6.21  X  lO'^)  (2  x  10®)  (18  x  lo"®)  =  224 

C  1^14 


T  =14.95  seconds 
c 

R  =  T  /C,  =  14,95/18  X  lO'®  =  8.30  x  10®  ohms 
1  c  1 

C,  =  T  /R,  =  14.95/2  X  10®  =  7.5  x  lO'®  farads. 

4  C  4 


Using  the  same  method  for  case  (b)  we  obtain 


R  C  =  t  =  D/U  =  3/10  =  0,30  second. 

X  4 


Rj  =  2  X  10  ohms 


Cj  =  18  X  10  farads. 


=  (0.30)(2x  10®)  (18  X  lO'®)  =  10800, 
c 


T  =  104  seconds, 
c 


Rl  =  T^/Cj  =  104/18  X  lO'®  =  5.8  x  10®  ohms. 


Cj  =  lo'®  farads. 


For  case  (c)  the  circuit  of  Figure  3. 23(g)  is  used  and  we  obtain 

Cg/Cjj  =  t®  =  2D/A  =  6/(50)  (32.2)  =  3.73  x  10‘®. 

9  -6 

Let  Ro  =  2  X  10  ohms  and  C,  =  18  x  10  farads, 
o  D 


Cg  =  ^  10'®)/(4  X  10^®)  =  16.8  X  lO'^’, 


Cg  =  2.56  X  10  farads 


T  =  R,C  =  (2  X  10®)  (2.56  x  10  ®)  «  5.12  seconds. 

C  O  u 


112 


For  the  equations  to  hold,  must  be  much  greater  than  R^.  Therefore  let 


Then 

and 

Since 


Rg  =  100  Rg. 

7 

Rg  =  2  X  10  ohms 

Cg  =  T^/(2  X  lO')  =  2.56  X  10‘^  farads, 

R3C2  =  4R^C^ 

Rj  =  (2  X  10®)  (2.56  X  10'V(4)(18  x  lO'®)  =  7.1  x  10®  ohms 

Cl  =  4R2Cg/R^  =  (4)(2xlo’)  (2.56  x  10'®)/(7.1  x  10®)  =  2,88  x  lo'®  farads. 


These  circuits  were  constructed  and  the  data  shown  in  Figure  3.22  were  obtained.  It  can  be 
seen  that  in  all  cases  waveforms  from  integrated  channels  closely  resemble  the  output  of  trant  - 
ducers  measuring  the  functions  directly.  In  addition  all  peak  data  are  accurate  to  well  within 
+10  percent. 
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Section  4 


SPHERICAL  SHOCK  WAVES  FROM  UNDERGROUND  EXPLOSIONS 

R.  H.  Bishop 

Introduction 

fii  order  to  improve  understanding  of  cratering  and  seismic-coupling  phenomena  —  especially 
for  extrapolating  from  chemical  to  nuclear  explosions  —  an  analysis  of  shock  waves  in  solid  media  has 
been  made.  Only  spherically  symmetric  waves  are  assumed,  in  uniform  unbounded  media;  these  waves 
are  analyzed  throughout  the  entire  range  of  distances  from  close  to  the  explosion  center  to  distances 
where  linear  elastic  theory  is  applicable.  Principal  results  derived  from  the  analysis  are  pressure- 
distance  decay  laws  for  various  geologic  materials.  Time  history  profiles  are  not  obtained,  but  the 
"effective  shell  thickness*  parameter  (which  is  obtained  as  a  function  of  the  radial  distance)  indicates 
the  general  scale  of  the  shock-wave  radial  profile,  without  reference  to  its  shape. 

The  partial  differential  equations  expressing  consei-vation  of  mass,  momentum  and  energy,  the 
constitutive  relations  of  the  medium  in  which  an  explosion  occurs,  and  the  appropriate  boundary  tmd 
Initial  conditions  are  in  principle  sufficient  to  describe  completely  the  wave  disturbance  generated  by 
an  explosion.  In  practice,  however,  a  general  and  complete  solution  is  extremely  difficult  mathemati¬ 
cally.  Considerable  progress  toward  a  solution  has  been  made  In  the  last  few  years  by  use  of  electronic 
comi'uters,  particularly  in  the  Immediate  or  h>  drodyrismlc  vicinity  of  the  e3cplosion.  Nevertheless,  at 
present  there  appears  to  be  no  computer  calculation  capable  of  predicting  values  for  an  explosion  dis¬ 
turbance  which  agree  with  experiments  in  the  variety  of  media  in  which  data  have  been  obtained  and 
over  the  complete  pressure  range  from  hydrodynamic  to  elastic. 

In  addition  to  the  mathematical  difficulties,  lack  of  knowledge  of  media  constitutive  relstions  is  a 
major  shortcoming.  Since  methods  of  solution  by  computer  techniques  are  relatively  difficult  and  not 
readily  available  for  general  use,  s  simple,  analytical,  and  easily  applied  method  of  solution  is  highly 
desirable.  It  is  the  aim  of  this  report,  first,  to  maks  those  simplifying  approximations  and  reasonable 
assumptions  for  the  unknown  constitutive  relations  of  the  medium  so  that  an  analytic  solution  may  be 
obtained;  and,  second,  to  evaluate  the  resulting  solutions  by  comparison  with  experimental  data. 
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In  the  present  report,  solutions  are  obtained  which  are  analytic  in  each  of  the  several  different 
pressure  ranges  associated  with  explosive  disturbances.  Over  each  range,  an  approximate  relation 
between  stress  and  strain  is  described  by  a  set  of  constant  parameters,  and  a  corresponding  analytic 
expression  relating  pressure  and  distance  is  obtained  which  holds  only  within  that  range. 

It  is  known  experimentally  and  from  theoretical  considerations  that  an  arbitrary  initial  pressure- 

distance  profile  generated  by  an  intense  initial  disturbance  in  a  medium  is  generally  not  maintained  during 

subsequent  propagation;  it  rapidly  approaches  a  certain  "steady- state"  ultimate  shape.  The  "front"  of 

this  ultimate  hydrodynamic  profile  is  very  steep,  almost  discontinuous,  and  the  "tall"  has  a  nearly 

22 

unchanging  shape  during  propagation,  being  roughly  exponential.  Taylor  (1950)  has  shown  that  the  shape 
of  the  hydrodynamically  stable  profile  of  a  sufficiently  strong  spherical  shock  wave  depends  only  on  the 
adiabatic  exponent  of  the  medium.  Moreover,  the  profile  remains  always  "similar"  --  1.  e. .  a  given 
fraction  of  the  front  pressure  is  always  located  at  the  same  relative  location  between  the  explosion  center 
and  the  shock  front  position  at  each  instant  of  time,  provided  the  shock  remains  sufficiently  strong. 

During  subsequent  propagation  at  lower  pressures,  two  changes  take  place;  (1)  the  scale  of  the  profile 
varies  as  the  front  moves  along;  (2)  the  shape  of  the  profile  varies  with  R.  As  a  simplifying  assumption, 
we  ignore  the  change  in  profile  shape  (2),  This  change  is  small  as  long  as  the  transition  from  ambient 
pressure  to  front  pressure  is  irreversible.  Kven  assuming  an  Invariant  profile  shape,  one  must  still 
take  account  of  the  change  of  scale  (1),  since  this  is  a  first-order  effect  and  cannot  be  neglected  even  in 
an  approximate  analysis.  It  is  useful  to  define  a  separate  parameter  to  account  for  the  change  in  scale 
of  the  radial  profile  which  occurs  during  propagation  at  relatively  low  front  pressures  -•  i.  e. .  when  the 
shock  is  not  "strong*  in  the  sense  defined  by  Taylor.  This  parameter  is  related  to  the  radial  distance 
(at  any  instant)  from  the  wave  front  to  a  point  on  the  profile  where  the  pressure  has  decreased  to  a  given 
fraction  of  the  front  pressure.  This  distance  defines  a  "shell  thickness."  y.  which  contains  most  of  the 
wave  energy. 

Use  of  the  "shell-thickness*  parameter  simplifies  the  mathematical  description  of  shock  propaga¬ 
tion.  The  fundamental  partial  differential  equations,  which  contain  time  as  an  independent  variable, 
are  required  to  describe  the  transient  approach  to  a  steacty  radial  profile  shape.  Assuming  a  steady 
shape  to  have  been  reached  at  some  early  time,  we  derive  ordinary  differential  equations  in  which  the 
independent  variable  is  the  wave  front  radius.  R,  at  any  instant  of  time.  Integration  of  these  equations 
gives  the  curve  at  pressui^  as  a  ftmetion  of  radial  distance  in  any  given  medium.  This  latter  curve  is 
the  principal  object  of  concern  in  the  following  analyses,  since  most  of  the  available  data  from  «q>eriments 
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can  be  compared  with  it.  In  order  to  derive  the  curve,  however,  the  constitutive  (stress- strain)  relations 
of  the  medium  must  be  known  in  each  pressure  range.  The  initial  steady  shape  of  the  radial  profile  is 
determined  from  Taylor's  analysis,  but  the  change  of  shape  which  actually  occurs  during  subsequent 
propagation  is  ignored  in  our  analytic  approximate  solutions. 


Stress-Strain  Description  of  a  Solid  Medium 


The  Hugoniot  curve  refers  to  the  set  of  all  shock  front  conditions  (p^,  p^)  that  can  be  reached  front 

22 

a  given  initial  condition  (pj,  p^^).  In  a  gas,  one  may  use  the  following  Hugoniot  (Taylor,  1950), 

7  +  1  +  (7  -  1)  (Pg/Pj) 

'^s^'^1  ^l^^s  “  7  -  1  +  (y  +  1)  (Pg/Pj) 

provided  that  the  parameter  7  is  defined  as 

y  -  1  +(P3V^  -  p,Vj)/(E^  -  E^)  , 

in  which  p  is  pressure,  v  is  specific  volume  or  reciprocal  of  density,  and  E  is  specific  internal  energy. 
The  subscripts  1  and  s  refer  to  regions  ahead  of  and  behind  the  shock  wave,  respectively. 

The  parameter  y  is  equivalent  to  the  exponent  of  a  reversible  adiabat  only  for  "weak"  shocks,  so 

that  p  -  p,  is  small  compared  to  p,.  In  general,  y  is  a  function  of  p  . 

Si  1  s 

In  terms  of  the  relative  fractional  change  of  volume,  or  "strain,"  e  =  1  -  over¬ 

pressure,  <7  =  p^  -  Pj,  the  Hugoniot  may  be  reduced  to 

1/e  =  (1/2)(7  +  1)  +yPj/CT  .  (4.1) 

When  a  is  large  compared  to  p^.  Equation  4.1  approaches  the  well  known  limit  Pg/P^  =  (y  +  l)/(y  -  D 
for  "strong"  shocks;  and  when  a  is  small  compared  to  p^  Equation  4.1  corresponds  to  the  differential 
form  of  the  reversible  adiabat,  (da/dp)^  =  y^PjPj  ^  »  in  which  y^  is  the  value  of  7  corresponding  to 
the  ambient  condition  (p^,  p^).  That  is,  7  approaches  7^  as  a  approaches  zero. 

Equation  4.1  can  be  applied  to  a  solid  medium,  but  the  parameter  y  may  vary  rapidly  with  the 

front  pressure,  a  .  For  solid  media,  we  make  the  assumption  that  the  product  ypj,  in  Equation  4.1  may 

be  replaced  by  the  reversible  adiabatic  elastic  modulus  S^.  The  quantity  y  in  this  product  must  also 

be  replaced  by  y^  .  By  definition,  =  Pj(do/dp)j  ,  and  it  is  well  known  that  the  plane  wave  propa- 

1  /2 

gation  velocity,  Cj  ,  corresponding  to  infinitesimal  longitudinal  strains  is  (da/dp)j^  .  Thus, 
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°  **1^1  •  E^iwatlon  4.1  is  then  expressed  In  terms  of  for  any  solid  medium, 

°  =  ^1<1  °  ®s  '  <4.2) 

in  which  x  is  defined  as 

X  =  (7  +  l)a/(2Pj^C^)  .  (4.3) 

The  stress-strain  ratio  a/e  in  Equation  4.2  defines  a  finite-strain  modulus,  or  "shock  modulus," 

.  Thus,  we  interpret  x  as  the  relative  fractional  correction  which  converts  the  elastic-wave 

modulus,  S.,  to  the  shock  modulus,  S  . 

1  s 

The  modulus  correction,  x,  in  Equation  4.2,  has  a  simple  expression  in  terms  of  the  shock-front 
propagation  velocity,  U.  This  is  derived  from  the  Rankine-Hugoniot  equation  for  conservation  of  mass 
across  a  shock  front,  Ug/U  =  e  ,  in  which  u^  is  the  particle  velocity  at  the  front.  Eliminating  e 
between  this  equation  and  Equation  4.2  and  using  the  momentum  equation,  a  =  p^u^U,  we  obtain 

X  =  (U/Cj)^  -  1  ,  (4.4) 

It  should  be  recognized  that  the  strain,  c,  in  Equation  4.2,  takes  place  in  a  direction  parallel  to  the 
directio  1  of  plane-wave  propagation,  so  that  there  is  no  change  of  volume  in  any  transverse  direction. 
However,  there  is  always  a  transverse  stress  field  accompanying  the  longitudinal  stress,  a  . 

If  X  is  negligible  compared  to  1,  each  stress  component  is  directly  proportional  to  the  longi¬ 
tudinal  strain,  e  .  In  this  case  the  strain  at  the  front  is  infinitesimal,  and  the  corresponding  waves 
are  said  to  be  "sonic"  or  "acoustic."  In  a  solid  medium,  the  transverse  stress  field  depends  on 
Poisson's  ratio  as  well  as  on  the  elastic  modulus  S^.  When  a  is  large  compared  to  and  x  is 
large  compared  to  1,  stresses  generated  in  a  solid  by  shock  waves  are  termed  "hydrodynamic."  When 
a  is  comparable  to  and  x  is  comparable  to  1,  the  quantity  a  is  neither  acoustic  nor  hydrodynamic, 
nor  is  it  simply  related  to  the  longitudinal  strain,  e  . 

Reversible  Strain 

Let  us  consider  a  spherical  compressive  wave  in  a  small  volume  of  rock  subjected  to  a  uniform 
lithostatic  pressure,  p^ .  This  pressure  is  due  to  the  weight  of  air  and  rock  at  depth  h^  below  the 
surface  and  is  approximately  equal  to  B  +  Pj^ghg.  where  Pj^  is  average  density  of  rock,  B  is  barometric 
pressure,  and  g  is  gravitational  acceleration.  As  a  spherical  shock  moves  along,  the  front  pressure 
will  eventually  decrease  to  a  value  less  than  that  necessary  to  open  radial  (hoop-stress)  cracks.  If  the 


118 


shock  overpressure,  ”  Ps  “  Pi'  replaced  by  a  static  overpressure  of  the  same  magnitude  at  the  wall 

O 

of  a  spherical  cavity.  Lamb  (1960)  has  shown  that  the  static  tangential  stress,  a^,  due  to  the  excess 
radial  stress,  is  given  by  =  (1/2)  -  2pj),  provided  that  the  medium  is  linearly  elastic.  When 

exceeds  2pj,  the  cavity  wall  is  subjected  to  a  tensile  hoop  stress.  When  reaches  some  critical 
value,  T,  radial  cracks  will  be  formed.  If  we  assume  that  this  static  relationship  can  be  extended  to  a 
shock  wave,  the  parameter  T  is  equivalent  to  a  dynamic  tensile  strength  in  the  medium.  Then  the  shock 
front  stress,  aj^,  corresponding  to  the  onset  of  radial  cracking  is  defined  by  the  equation 
T=  (1/2)  (c^  -  2p^),  or  =  2(T  +  p^)  . 

The  actual  maximum  reversible  stress,  p^,  cannot  exceed  this  limit,  but  it  may  have  some  smaller  value. 
Thus,  in  general,  p^  is  less  than  or  equal  to  2(T  +  p^). 

2 

Since  =  ^i^l'  reversible  adiabatic  exponent  in  a  solid  medium  at  a  moderate 

3 

ambient  pressure,  p^,  may  be  relatively  large.  For  example,  if  =  2  gm/cm  ,  Cj^  =  4  km/sec  and 
Pj  =  GO  bars  (1000-foot  depth),  y  >  10^. 

When  X  is  small  compared  to  1,  the  strain  is  small  and  y  approaches  the  limit,  yj^  ,  Since  y^^ 
is  generally  large  compared  to  1  in  any  solid  medium,  and  =  7jPj»  Equation  4.3  is  practically 

equivalent  to 

X  =  o’/2p^  (a<p^)  .  (4.5) 

When  a  is  less  than  p^,  the  relation  between  stress,  0 ,  and  strain,  c,  is  given  by  Equations 
4.2  and  4.5.  Evidently,  the  stress  is  directly  proportional  to  the  strain  only  when  x  is  negligible  com¬ 
pared  to  1.  In  general,  the  relation  may  be  nonlinear  even  when  the  strain  is  reversible  (0  <  p^), 
because  x  is  not  necessarily  negligible  compared  to  1,  To  see  this,  we  observe  from  Equation  4.5 
that  X  may  have  any  value  less  than  and  p^,  in  turn,  may  have  any  value  less  than  2(p^  +  T). 

Thus,  X  may  have  any  value  less  than  1  +  T/p^. 

Cracking  Zone  and  Crushing  Zone 

When  the  radial  stress,  0j^,  exceeds  the  reversible  limit,  p^,  radial  cracks  will  appear  first, 
due  to  tangential  tensile  stresses  at  the  spherical  wave  front.  In  a  range  of  higher  values  of  peak  radial 
stress,  tangential  cracks  may  be  expected  to  form  as  a  result  of  shear  stress  failures.  Consequently, 
separate  rock  fragments  will  result  from  a  combination  of  processes  resulting  in  tangential  and  radial 
cracks.  As  the  radial  stress  increases,  the. fragments  may  break  up,  until  a  granular  composition  is 
left  as  the  "crushing- strength,"  p^.  Is  approached.  We  call  the  effective  dynamic  compressive 
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strength.  If  the  peak  radial  stress,  is  less  than  p^,  each  fragment  tends  to  retain  its  shape;  in 

this  i  ange  of  stress  there  is  no  excessive  transverse  bulging,  and  the  tangential  compressive  stress 

component  is  negligible  compared  to  the  radial  stress  component.  After  the  first  radial  cracks  have 

appeared,  the  tangential  stresses  become  zero,  and  the  radial  stress  component  becomes  dominant. 

All  these  conditions  apply  approximately  to  any  radial  stress,  which  is  less  than  the  dynamic 

compressive  strength,  p  .  The  region  of  a  medium  in  which  the  stress  a_  has  values  such  that 

C  R 

Pr  <  <  Pc  is  termed  the  "cracking  zone,"  and  in  this  zone  the  radial  stress  component  is  considered  as 

the  only  stress  of  importance  in  a  first  approximation. 


Because  the  radial  stress,  is  directed  perpendicularly  to  any  concentric  spherical  surface,  it 

will  diminish  inversely  as  the  square  of  any  radial  distance  in  the  cracking  zone.  This  is  readily  seen 
from  the  equilibrium  equation  in  spherical  coordinates. 
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in  which  the  tangential  stress,  a,,  is  set  equal  to  zero.  The  radial  stress,  a^,  which  is  a  solution  of 

t  R 

this  equation,  varies  inversely  as  the  square  of  the  radial  distance,  R. 


This  variation  is  strictly  true  only  for  a  static  radial  stress  field.  We  shall  assume  that  the 

inverse-square  law  is  also  obeyed  for  wave-front  stresses  in  the  cracking  zone,  (p^  <  <  p^).  It  should 

be  recognized  that  if  the  radial  stress,  Oj^,  exceeds  the  crushing  strength,  p^,  the  transverse  stress 

components  can  no  longer  be  neglected  even  in  a  first  approximation.  In  the  cracking  zone,  each  rock 

fragment  tends  to  resist  a  change  of  shape,  which  is  described  by  a  finite- strain  shear  modulus,  . 

As  the  radial  stress  increases  in  the  cracking  zone,  resistance  to  change  of  shape  begins  to  disappear  as 

fragments  are  broken  down  into  smellier  and  smaller  particles  of  more  regular  shape.  These  ultimate 

particles  tend  to  resist  change  of  volume  rather  than  change  of  shape;  hence,  the  shear  modulus,  p  , 

s 

becomes  negligible  in  that  range  of  stress  which  exceeds  the  crushing  strength,  p^;  this  range  will  be 

termed  the  "emshing  zone."  Now,  the  stress-strain  relation  can  be  described  by  a  finite-strain  shock- 

wave  bulk  modulus,  k  . 

s 

In  order  to  determine  a  stress-strain  relation  for  a  medium  which  applies  to  both  the  cracking  zone 

and  the  crushing  zone,  we  express  the  dynamic  shock  modulus,  S^,  in  Equation  4.2  in  terms  of  the 

shock-wave  bulk  modulus,  k  ,  and  the  shock- wave  shear  modulus,  p  .  By  analogy  with  the  acoustic 

s  s 

2 

modulus,  Sj,  representative  of  infinitesimal  strains  and  expressed  by  Sj  =  k  +  (4/3)p  =  ,  we 

define  finite  strain  or  shock-wave  moduli  by  the  corresponding  relation,  S  =  k  +  (4/3)p  .  Using 

s  s  s 


2 

the  well  known  relation  ji  =  PjCg  ,  where  Cg  Is  the  velocity  of  infinitesimal  shear  waves,  the  acoustic  ■ 
bulk  modulus,  k ,  can  be  expressed  in  terms  of  the  acoustic  plane  wave  modulus,  Sj^: 

k/S^  =  1  -  (4/3)(C2/C^)^  .  (4.6) 

In  the  cracking  zone,  it  is  expected  that  Poisson's  ratio  remains  nearly  the  same  as  its  "acoustic" 
value.  Thus,  k  /p  =  k/p,  k  /k  =  p  /p  =  f,  and  S  =  fS,,  where  f  is  a  constant  factor.  Solving 
Equation  4.2  for  x  gives 

x^  =  (4/3)(Pg/S^)  +  f(k/S^)  -  1  ,  (4.7) 

in  which  the  subscript  c  refers  to  the  crushing  zone  or  the  cracking  zone.  In  the  crushing  zone,  as 
already  explained,  the  finite-strain  shock-wave  shear  modulus,  p^,  becomes  negligible,  euid 
Equation  4.7  reduces  approximately  to 

x^  =  f  [^1  -  (4/3)(C2/C^)^j  -  1  .  (4.8) 

In  the  crushing  zone,  the  factor  f  is  not  strictly  constant,  as  in  the  cracking  zone,  but  depends  on 

the  amplitude  of  the  stress  a  .  In  the  absence  of  experimental  data,  we  shall  assume  Equation  4.8  in 

conjunction  with  Equation  4.2,  for  the  stress-strain  relation  in  the  crushing  zone.  In  Equation  4.8, 

the  factor  f  is  equal  to  the  ratio  5  IS.,  in  which  S  is  the  "rapid-strain"  shock  modulus,  and  S, 

SIS  1 

26 

is  the  "slow-strain"  acoustic  modulus.  It  has  been  observed  by  Wuerker  (1959)  for  "mine  rocks" 

24 

and  by  Watstein  (1953)  for  concrete  that  the  ratio  of  a  rapid-strain  modulus  to  its  corresponding 
slow-strain  modulus  is  about  2.  In  accordance  with  these  experiments,  we  use  f  =  2  in  Equation  4.8, 
so  that  the  modulus  correction  x^  in  the  crushing  zone  becomes 

x^  -  1  -  .  (4.9) 

Hydrodynamic  and  Plastic  Ranges 

When  the  front  pressure,  a,  exceeds  a  certain  minimum  value,  p^,  individual  particles  merge 

to  produce  a  continuous  fluid  with  high  viscosity,  which  exhibits  "plastic  flow."  When  the  pressure 

exceeds  a  still  higher  limit,  Pjj,  the  medium  displays  a  completely  hydrodynamic  character  as  a  low 

viscosity  fluid.  We  can  determine  an  approximate  value  for  p^  by  substituting  x  =  x^  into 

2 

Equation  4.3,  and  replacing  o  by  p^,  which  gives  p^  =  2x^p^C^Uy  +  1).  This  result  follows  from  the 
fact  that  we  have  identified  the  highest  pressure  in  the  crushing  zone  with  the  lowest  pressure  in  the 
plastic  zone.  The  expression  for  p^  involves  the  parameter  y.  From  experimental  Hugoniot  data, 
y  may  be  evaluated  using  Equations  4.2  and  4.3.  It  is  foxind  that  y  is  about  2  or  3  in  most  geologic 
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matex'ials.  These  values  for  y  are  accurate  in  the  hydrodynamic  range  of  pressure,  a  >  p„  ,  but 

approximate  in  the  plastic  range,  p  <  ct  <  p„  . 

m  H 

In  the  hydrodynamic  pressure  range,  a  >  p^^,  y  is  very  nearly  constant.  In  the  plastic  range,  y 
is  a  slowly  varying  function  of  a  ,  so  that  in  subsequent  calculations  we  shall  take  it  to  be  constant  for  ; 
all  pressures,  a  >  p^,  and  equal  to  its  value  when  o  >  Pj^  . 

Pressure  Zones  Near  an  Explosion 

The  schematic  diagram  of  Figure  4.1  illustrates  the  various  pressure  regions  surrounding  an 
luiderground  explosion.  The  modulus  correction,  x,  in  these  same  pressure  regions  is  shown  as  a 
function  of  ct  in  Figure  4.2  for  a  typical  earth  material.  The  information  of  Figure  4.2  has  been 
converted  to  a  stress-strain  diagram  in  Figure  4.3. 

Table  4.1  summarizes  the  modulus  corrections  and  the  stress-strain  relations  for  the  various 
pressure  ranges. 

The  only  certain  information  we  have  about  the  properties  of  solid  media  is  their  elastic  constants 
and  their  Hugoniot,  illustrated  in  Figure  4.2  by  the  two  regions  of  the  curve  where  the  slope  is  unity. 
Retween  these  extreme  limits  of  pressure  our  knowledge  of  medium  properties  is  imprecise,  if  not 
negligible.  The  two  straight-line  portions  of  the  curve  in  the  acoustic  and  fluid  pressure  ranges  must  be 
connected  in  some  fashion.  We  have  made  this  connection  by  assumption,  using  the  relations  listed  in 
Table  4,1. 
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Figuri’  4,1  Pressure  Zones  Near  eui  Underground  Explosion 
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MODULUS  CORRECTION  LOS 


Figure  4.2  Variation  of  Modulus  Correction  x  With  Shock  Front  Overpressure  a 


STRAIN  C 


Figure  4,3  Variation  of  Stress  a  With  Strain  e  at  a  Shock  Front 
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TABLE  4.1 


Modulus  Correction  and  Stress-Strain  Relation  in  Various  Pressure  Ranges 


Ccmdition  of  medium  and 
Zone  range  of  pressure _ 

I  Continuous  solid,  reversible 

strains,  o  <  p 


Modulus  correction  x 


X  = 


2P, 


Stress-strain  relation, 
_ g  versus  g _ 


n  Cracked  solid,  large 

fragments  or  radial  cracks. 

p  <  a  <  p 
r  c 


f  >2 
c 


1  1/e  = 


m  Crushed  solid,  small 
particles,  p^  <  a  <  p^ 


1  l/e  = 


f  =  2 


IV  Fluid,  hydrodynamic  or 


plastic,  a  >  p 


m 


{y  *  l)g 

2p,C^ 

I  A 


1/e  =  (^)  (Y  +  1)  +  p^C^/a 


Shell  Thickness  and  Wave  Energy 

Let  us  take  the  shock-wave  energy,  W^.  as  twice  the  kinetic  energy  in  a  spherical  volume  of 
radius  R,  the  distance  from  explosion  center  to  shock  front.  Near  the  front,  the  kinetic  energy  per 
unit  mass  is  known  to  be  practically  equal  to  the  internal  energy  per  unit  mass,  provided  the  front 
pressure  is  large  compared  to  ambient.  Thus 


W  =  4)r 
s 


/•R 


2  2  . 
r  dr 
r 


in  which  is  density  and  u^^  is  particle  velocity  at  radial  distance  r  from  the  center  of  explosion. 

We  now  define  the  effective  shell  volume,  V_.  to  be  equal  to  W_  /(p  ,  in  which  p.  is  density  and  u_ 

8  8/  \  S  8/  8  S 

is  particle  velocity  at  the  wave  front.  Then 


V  =  4irf  {p  jp  )  (u  /u  J^r^dr 
s  Jn  r  B  r  8 
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The  shell  volume.  V  .  can  be  expressed  in  terms  of  its  thickness,  y.  and  outer  radius  R,  so  that 
8 

Vg  =  4»R*yF.  in  which 

F  =  1  -  (y/R)  +  (1/3)  (y/R)*  .  (4.10) 

Now.  letting  n  =  r/R.  the  shell-thickness  ratio  y/R  is  found  from 

y/R  =  (1/F)j^  (p^/Pg)(u^/u^)*n*dn  .  (4.11) 

Equation  4.11  applies  generally  to  aqy  spherical  shock  wave,  provided  the  front  pressure  is  large 
compared  to  the  ambient  pressure,  p^  . 

R  follows  from  replacement  of  yp^  by  'in  Equatiem  4.1  that  Pj/P^  approaches  (y  -  l)/(y  -i-  1) 

when  o  is  large  compared  to  s^  .  From  Equation  4.3.  this  same  condition  results  in  a  large  value 

22 

of  X.  regardless  of  whether  a  gas  or  a  solid  is  being  considered.  Taylor  (1950)  in  his  analysis 
assumes  the  asymptotic  front  condition  P^/Pg  =  (7  -  l)/(7  +  D;  accordingly,  Taylor's  analysis  can  be 
applied  to  a  solid  medium  when  x  is  large  compared  to  1  and  the  shock  wave  is  strong. 

Taylor's  theory  provides  radial  profiles  for  inclusicm  in  Equation  4.11.  In  terms  of  Taylor's 
variables,  ^(n,  y)  and  ^n,  7),  the  result  is 

(y/R],  =  .  (4.12) 

in  which  the  subscript,  o,  designates  a  strong  shock  wave.  The  values  of  (y/K)^  calculated 
from  Equations  4.10  and  4,12  and  Taylor's  formulas  for  ^  and  ^  appear  in  Figure  4.4.  This  shows 
that  the  hydrodynamic  shell- thickness  ratio  is  a  function  only  of  the  medium  parameter,  >, 

and  does  not  depend  on  R. 


Figure  4.4  Relative  Shell  Thickness  Versus  Adiabatic  Exponent  for  Hydrodynamic  Shock  Waves 
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S  is  sem  (Figure  4.4)  that  the  shell-thickness  ratio  (y/R)^  is  invariant  as  R  increases,  provided 
y  remains  constant  and  the  shock  wave  remains  strong;  whoi  x  is  large  compared  to  1.  a  diock  is 
defined  to  be  "strong,"  regardless  of  the  type  of  medium.  Any  spherical  blast  wave  can  be  described 


adequately  by  two  parameters,  namely,  front  pressure,  o.  and  shock  energy,  W  .  The  shock  energy 

s 

2 

can  also  be  caressed  in  terms  of  the  relative  shell  thickness.  y/R,  and  the  energy  density,  .  at 

3  2  2 

the  fnmt.  From  preceding  definitions,  W  =  4  v  R  p  u  Fy/R.  Energy  density,  p  u  .  is  related  to 

S  S  8  8  S 

2 

front  pressure,  o.  since  the  change  of  kinetic  mergy  is  (l/2)u^  =  (l/2)o(v^  -  v^)  (e/2)(<r/p^)  . 

22 

Taylor  (19S0)  has  shown  that  there  is  a  simple  relati<m  between  the  front  pressure,  o.  and  the 
total  blast  energy,  .  This  quantity  is  defined  as  the  total  kinetic  and  internal  energy  associated 
with  all  the  material  inside  the  entire  sphere  of  radius  R.  R  is  sometimes  convenient  to  separate  the 


total  blast  energy,  W  ,  Into  two  parts.  The  shock-wave  energy,  W  ,  has  alreacfy  been  described.  It 
t  8 

is  associated  with  the  spherical  shell  whose  outer  radius  is  R  and  whose  inner  radius  is  R  -  y  .  We 
now  define  the  "core"  energy.  W^,  to  be  equal  to  ~  core  energy,  W^,  is  associated  with 

the  spherical  volume  whose  radius  is  R  -  y  . 

22 

Taylor  (1950)  defines  a  dimensionless  pressure,  f{n,  y),  analogous  to  the  relative  density,  tp, 
and  relative  particle  velocity,  which  we  have  already  used  in  Equation  4.12.  From  Taylor's  result 
in  our  notations,  the  front  pressure,  o ,  is  given  by 


where 


a  =  B(y)WjR 


B(y)  =  )r(7+  1)  L  +  (2/7)  (7-  1)' 


and  X  >>  1.  It  is  calculated  from  two  definite  integrals  through  the  use  of  Taylor's  formulas  for  ^,0. 


and  f . 


Ij(y)  "  ^  (('4*n^dn  , 
12(7)  •  ^  f(n,7)n*dn  . 


The  functicm  8(7)  It  plotted  in  Figure  4.5 
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ADIABATIC  EXPONENT/ 


Figure  4.5  Relative  Pressure  Versus  Adiabatic  Exponent  for  Strong  Shock  Waves 


From  the  strong- shock  energy  relations.  =  4irR^PgUgF(y/R)^.  u^  =eo/pj,  c  =  2/(y  +  1),  and 
Equations  4.12  and  4.13,  we  find  that  =  2  x(y  +  l)R^Ij(7)a .  W^/W^  =  2n(y  +  1)Ij(y)B(y).  and  also  that 
W^/W^  =  1  -  Wg/W^  .  The  ratio  W^/W^  designates  the  fraction  of  the  total  energy  which  is  partitioned  to 
shock  propagation  (i.e. ,  in  the  shell),  and  the  remainder  is  retained  by  the  core.  These  fractions  depend 
on  the  medium;  for  strong  shocks  (x  »  1)  they  are  functions  of  y  only.  It  is  also  evident  that  for  strong 
shocks  the  fraction  of  energy  partitioned  to  shock  propagation  is  independent  of  the  front  radius.  R, 
provided  y  is  constant.  For  example,  when  y  -  3,  W^/W^  -  0.8  .  Thus,  only  20  percent  of  the  total 
energy  is  retained  by  the  core  in  this  case;  whereas,  when  y  •  1.4,  W^/W^  ■  0,42  aiid  nearly  60  percent 
of  the  energy  is  retained  by  the  core. 


Heat  Energy  Loss  From  Shell 

As  shown  in  the  last  section,  there  is  no  significant  loss  of  energy  from  a  strong  spherical  shock¬ 
wave  shell,  since  y  is  nearly  constant.  As  the  spherical  front  radius.  R,  increases,  however,  the 
shock  necessarily  becomes  less  stnmg.  and  eventually  there  is  an  appreciable  net  loss  of  energy  from  the 
shell.  Then  becomes  a  variable,  depending  on  the  shock  front  radius  R  at  any  given  Instant.  Since 
the  adiabatic  e]q>anslon  curve  of  a  medium  lies  above  the  Hugonlot  curve  (Figure  4.6),  the  net  loss. 

(per  unit  mass),  is  always  less  than  the  heat,  Q,  which  corresponds  to  a  Hugonlot  expansion  curve,  in 
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which  the  final  specific  volume  is  the  same  as  the  initial  specific  volume.  Let  us  define 


Q  =  G  G  Q  ,  (4.1 

s  m 

where 

Q  =  heat  which  is  ultimately  lost  from  the  shell  and  retained  by  the  core; 
s 

Q  =  heat  associated  with  the  Hugoniot  expansion  curve  of  Figure  4.6; 

G^  =  factor  which  depends  on  shock  strength;  and 
G^  =  factor  which  depaids  on  type  of  medium. 

As  we  have  seen  for  strong  shock  waves,  the  shell  energy,  W  ,  remedns  constant  as  the  front 

s 

moves  along.  Therefore,  for  strong  shocks,  must  be  zero,  although  Q  is  never  zero.  For  shocks 
which  are  not  strong,  neither  or  Q  is  zero. 


Figure  4.6  Sioek  Energy  Relatione;  Pressure  Versus  Volume 


s 


The  heat  Q  is  the  approximate  energy  deposited  in  a  solid  medium  if  or  <  p^  ,  since  the  final 

specific  volume  is  nearly  equal  to  the  initial  specific  volume.  It  is  to  account  for  these  observations 

that  it  is  convenient  to  introduce  the  two  parameters  and  .  Also,  by  the  Introduction  of  parameters 

G^  and  G^  explicit  consideration  of  adiabats,  which  are  completely  unknown  for  geologic  media,  is 

avoided.  The  shock  strength  parameter  G^  must  go  to  zero  as  the  shock  becomes  strong  in  order  to 

account  for  Taylor's  result  that  the  shell  energy,  W^,  remains  constant.  The  factor  G^^  is  a  medium 

parameter  which  can  be  evaluated  approximately  from  an  experimental  determination  of  the  Hugoniot  in 

the  plastic  zone.  Explicit  forms  for  G  and  G  will  be  presented  in  the  following  sections. 

ni  X 

Referring  now  to  Figure  4.6,  the  heat  energy  Q  per  unit  mass  is  illustrated  graphically  on  a 
pressure-volume  diagram  as  the  area  between  the  Rayleigh  line,  AC,  and  the  Hugoniot  curve.  The  area 
enclosed  represents  the  approximate  net  hea^  energy,  Q,  deposited  in  a  unit  mass  of  condensed  material. 
Using  the  analytic  Hugoniot  curve  described  by  Equations  4.2  and  4.3,  the  heat  fraction  is  given  by 

h  =  Q/Ug  =  1/2  -  (l/x)^(l  +  1/x)  In  (1  +  x)  -  ij  ,  (4.15) 

2 

in  which  Ug  =  a(v^  -  v^)  *  oc/p^  represents  twice  the  kinetic  energy  at  the  shock  front,  and  is 
shown  in  Figure  4.6  as  the  rectangle  ABCD. 

When  X  is  small  compared  to  1,  Equation  4.15  reduces  approximately  to  h  °  x/6:  and  h 
approaches  1/2  when  x  is  large  compared  to  1. 


Energy  Propagation 

From  the  definition  of  the  effective  shell  volume  V^,  .  Differentiating  this  with 

respect  to  the  front  radius  R, 


dW^/dR  =■  ^J/dRj  +  uJ(dV^/dR)j  . 

2 

in  which  we  neglect  the  change  in  front  density  p^  compared  to  changes  in  u^  and  .  Consider  the 
energy  loss  UW^  from  the  shell  which  results  when  a  shock  front  advances  from  radius  R  to  radius 

9 

R  4  dR.  Then  dW^/dR  >  -4vR  ,  in  which  p^  is  the  ambient  density  oi  the  medium.  Eliminating 

dW^/dR  from  the  last  two  equations  provides  a  di^erential  equation  describing  variation  of  shock  energy. 
2 

u^,  with  dlstence,  R.  Using  Equation  4.14  and  introducing  s  •  R/y  we  obtain 


(4.16) 


lao 


(4.17) 


q  =  (l/F)zhG  G  (p,/p  )  . 

m  X  I  s 

in  which  q  is  a  parameter  which  denotes  the  heat  energy  loss  from  the  shock-wave  shell.  The  variable 
z  represents  the  ratio  of  the  front  radius  R  to  the  effective  shell  thickness  y. 

Shell  Thickness  In  Hydrodynamic  Zone 

In  order  to  integrate  Equation  4.16  in  the  hydrodynamic  zone  where  y  is  nearly  constant,  we  use 

2 

the  Hugoniot  Equations  4.2  and  4.3  and  eliminate  c  through  the  energy  relation  u^  =  oc/pj  .  This  gives 
2  -1  -1  -1 

Ug  =  Of  apj  (1  +  1/x)  in  which  a  =  (l/2)(7  +  1)  and  x  =  aalSy  From  these  results, 

2-2-12-1  2  -1 

Ug  =  «  Pji  SjX  (1  +  x)  ,  and  d  log  u^/d  log  x  =  (2  +  x)(l  +  x)  ,  when  y  is  constant.  Now  we  can 

obtain  a  differential  equation  from  Equation  4.16  which  relates  the  variation  of  x  with  R  to  the  variation 

of  z  with  R  in  the  hydrodynamic  zone. 


It  is  evident  that  this  equation  cannot  be  integrated  until  z  is  expressed  in  terms  of  x  and  R,  Let  us 
assume  that 

dlog  z/dlogR=  1/(1 +  x)  ,  (4,18 

and  also  that  this  relationship  is  valid  over  all  ranges  of  pressure.  It  is  apparent  that  Equation  4.18 
applies  to  strong  shock  waves,  where  x  is  large  compared  to  1,  since  we  have  already  seen  (Figure  4.4) 
that  Taylor's  theory  requires  z  to  be  constant  if  7  is  constant;  and  so  for  very  strong  shock  waves, 
d  log  z/d  log  R  =  0.  Very  weak  shock  waves,  or  acoustic  waves,  correspond  to  values  of  x  which  are 
very  small  compared  to  1.  These  also  are  described  by  Equation  4.18,  since  z  =  R/y  and  y  remains 
constant,  independent  of  R  so  that  d  log  z/d  log  R  =  1.  Each  "infinitesimal"  Fourier  component  of  a 
very  weak  disturbance  travels  at  the  same  "acoustic"  velocity,  Cj^,  as  is  evident  from  Equation  4.4, 
since  x  approaches  zero  for  acoustic  waves.  Therefore,  an  acoustic  wave  pulse  does  not  disperse;  it 
maintains  a  constant  radial  extension  whose  scale  is  measured  by  the  effective  shell  thickness,  y  . 

Using  Equation  4.18,  we  obtain  a  differential  equation  which  describes  the  variation  of  x  with  R 
in  the  hydrodynamic  zone. 


/d  log  x\  _  2  -t  3x  +  q(l  +  x) 
\d  log  Rj  ^  2  +  X 


where  q  is  defined  by  Equation  4.17,  and  F  has  been  set  equal  to  1.  Referring  to  Equation  4.10  and 
Figure  4.4,  we  see  that  y/R  is  relatively  small  for  strong  shock  waves,  and  F  is  nearly  equal  to  1  . 
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The  ratio  y/R  will  become  even  smaller  for  weak  shock  waves,  because  y  tends  to  remain  constant 
when  R  increases,  as  already  pointed  out.  Therefore,  we  can  assume  F  =  1  for  any  shock  wave,  with 
adequate  accuracy. 

By  eliminating  R,  a  differential  equation  for  z  is  obtained  from  Equations  4.18  and  4.19,  which 
is  valid  in  the  hydrodynamic  pressure  range  (7  =  constant), 

20) 

Now  we  refer  to  the  definition  of  q  in  Equation  4.17,  and  set  F  =  1.  Each  remaining  factor  in 
Equation  4.17  must  be  evaluated  in  order  to  integrate  Equation  4.20.  First  of  all,  it  is  evident  that  q 
is  directly  proportional  to  z,  and  the  differential  Equation  4.20  is  consequently  nonlinear  in  the 
dependent  variable  z.  Moreover,  it  will  become  apparent  that  the  remaining  factors  defining  q  in 
Equation  4,17  are  complicated  functions  of  the  independent  variable  x.  An  approximate  analytic  method 
is  desired  in  order  to  avoid  numerical  integration  of  Equation  4.20.  Fortunately,  one  can  derive  an 
analytical  approximation  which  is  sufficiently  accurate  for  most  practical  purposes. 


d  log  z  _  _  2  +  X _ 

d  log  X  (1  +  xU2  +  3x  +  q(l  +  x) 


Loss  Factor,  q,  in  Hydrodynamic  Zone 

It  can  be  shown  from  Equation  4.20  that  z  is  a  slowly  varying  function  of  x  when  x  is  greater 
than  about  5,  regardless  of  the  variation  of  q  with  x.  This  corresponds  generally  to  the  hydrodynamic 
pressure  range  described  in  Section  2.  Consider  now  each  of  the  remaining  factors  in  Equation  4. 1 7 
which  determine  the  shell  energy  loss  parameter  q.  In  Equation  4.17,  F  =  1  and  h  is  a  function  only 
of  X  given  by  Equation  4.15.  From  Equation  4.15  it  is  seen  that  h  is  a  slowly  varying  function  of  x 
when  X  is  greater  than  5.  From  Equations  4.2  and  4.3 


Pj/p^  =  1  -  2(1  +  7)’^1  +  l/x)"^ 


(4.21) 


showing  that  <Pj/Pg)  is  slso  a  slowly  varying  function  of  x  when  x  is  greater  than  5,  and  when  7 
is  between  2  and  3,  values  typical  of  hydrodynamic  pressures  in  solids. 


In  Equation  4.17  and  Equation  4.14,  is  a  shell  energy- loss  factor  which  depends  on  shock 
strength,  and  depends  on  the  type  of  medium.  We  have  seen  that  G^  approaches  zero  as  x 
becomes  large.  By  analogy  with  Equation  4.18,  let  us  assume  that  G^  =  1/(1  +  x),  so  that  Equation 
4.14  becomes 


Gj^  =  (1  +  x)(Q^/Q)  , 


(4.22) 
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and  Equation  4.17  becomes 


(1  +  x)q  =  G^zh(Pj/p^)  =  q^  .  (4.23) 

From  Equation  4.22  we  see  that  is  a  slowly  varying  function  of  x  in  the  hydrodynamic  zone. 

Note  that  G^  cannot  become  large  as  rapidly  as  x,  because  Qg/Q  must  approach  zero  as  x  increases. 

This  observation  indicates  that  qg  in  Equation  4.23  is  a  slowly  varying  function  of  x,  since  each  of  the 

other  factors  has  already  been  shown  to  be  a  slowly  varying  function  of  x.  As  a  result,  we  assume  that 

some  constant  average  value  ^2!  =  q  can  be  used  in  evaluating  Equation  4.20, 

'  'av 

/d  log  z\  _  _ 2  +  X _  , 

\d  log  xj  (1  +  x)  (2  +  q^  +  3x) 

which  can  be  integrated  to  give 


z/z.  = 


1  +  1 


1  + 


i,\= 


,3(b  -  1) 


h  +  b/x  V 

p^i/ 


b  -  2 
3b(b  -  1) 


(4.24) 


in  which 


b  =  (1/3)  (2  +  q^)  . 

In  Equation  4.24,  z  =  z^^  when  x  =  x.,  and  in  Equation  4.25, 


(4.25) 


(4,26) 


In  Equation  4.26,  we  first  examine  the  variation  with  respect  to  x  of  each  of  the  factors  over  the 
hydrodynamic  range  of  pressures. 

In  order  to  determine  z  in  Equation  4.2G,  let  us  find  an  average  value  for  z  over  the  hydro- 
dynamic  range  of  x.  To  do  this,  it  is  sufficient  to  examine  Equation  4.24  for  all  values  of  x  in  the 
hydrodynamic  range  which  are  smaller  than  x  =  100.  It  is  unnecessary  to  consider  Equation  4.24 
for  larger  values  of  x  since  z  remains  practically  the  same  as  the  strong  shock  limit,  (y/lD^  “  > 

given  in  Figure  4.4.  This  conclusion  also  follows  necessarily  from  Taylor's  theory,  because  the 
actual  density  ratio  Pj/Pg  given  by  Equation  4.21  differs  from  Taylor's  limit  (7  -  l)/(7  +  1)  by  less 
than  1  percent  when  x  is  greater  than  100. 

From  Figure  4.2  we  find  that  x^  is  the  lowest  values  of  x  in  the  plastic  zone,  and  from  Table  4,1 
it  is  seen  that  the  same  stress-strain  relation  is  used  in  both  the  plastic  zone  and  the  hydrodynamic 
zone.  In  Figure  4.2,  the  transition  between  hydrodynamic  and  plastic  pressures  occurs  where  x  =  Xj^  . 
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Below  this  point  y  is  not  constant,  and  the  hydrodynamic  equations  in  which  a  constant  value  of  y  is  not 
assumed  are  not  strictly  valid.  If  7  is  taken  to  be  constant,  as  a  simplifying  assumption  in  the  plastic 
range,  (x^  <  x  <  x^j),  we  can  make  an  additional  simplifying  assumption  in  order  to  avoid  calculating  a 
separate  lower  limit  x^  for  the  hydrodynamic-plastic  range  in  each  medium.  We  shall  assume  Xg  =  1/2 
for  the  "effective-average"  lower  limit  of  this  range  in  each  medium,  and 'the  corresponding  pressure  will 
be  called  the  "hydrodynamic  termination  pressure,"  Pg  (see  Figure  4.2).  As  determined  from  Equation 
4.3,  we  obtain  Pg  =  S^/(7  +  1).  This  may  not  be  the  "best"  choice  that  can  be  made.  Ideally,  one  should 
perform  a  separate  analysis  of  the  (curved)  plastic  zone  and  determine  the  variation  of  x  with  a  without 
restricting  it  to  follow  a  linear  relationship  (i.e. ,  constant  7)  as  we  have  just  done.  However,  the 
parameters  of  an  extensive  geologic  formation  are  not  sufficiently  uniform  to  make  such  an  effort  reward¬ 
ing,  and  our  objective  here  is  to  make  reasonable  and  simplifying  assumptions  so  that  the  problem  is 
tractable  in  analytic  form. 

From  Equation  4.25,  it  is  evident  that  b  is  always  greater  than  2/3,  since  q^  is  never  zero  in  the 
hydrodynamic  range  of  pressures.  It  follows  from  Equation  4.24  that  the  greatest  possible  values  of  z/z. 
occur  when  b  has  the  smallest  possible  value,  namely  b  =  2/3. 

Therefore,  evaluating  Equation  4.24  between  limits  x^  =  100  and  Xg  =  1/2,  and  assuming  b  =  2/3, 
we  obtain  the  greatest  possible  value  for  Zg,  namely  Zg  =  1.  8z^  .  Here  Zg  is  the  value  of  z  at  the  hydro- 
dynamic  termination  pressure,  pg,  and  1/z.  =  asymptotic  limit  (y/R)^,  shown  in  Figure  4.4, 

at  the  high-pressure  end  of  the  hydrodynamic  range,  where  x.  =  100.  An  average  value  of  z  can  be  taken 
to  be  half  way  between  l-8z^  and  z^,  so  that  z^  =  .  Because  of  the  assumption  b  =  2/3,  this 

particular  average  is  an  upper  limit  for  any  medium,  which  is  approached  only  by  a  nearly  loss-free 
medium  such  as  halite  or  granite.  From  experimental  data  we  find  that  x.^  =  3,5  and  b  =  3,28  in  tuff, 

ri 

which  differs  considerably  from  the  lower  limit,  b  =  2/3.  In  tuff,  it  is  found  that  z^  =  l.lGz^  .  For 

desert  alluvium,  b  =  12.8,  and  z  is  less  than  1.1 6z  .  It  turns  out  as  a  result  that  an  average  value  for 

a  o 

z  in  any  rock  medium,  including  alluvium,  is  given  by  z  =  1.2z  with  adequate  accuracy.  This  is  the 

H  O 

value  we  sliall  use  in  Equation  4,26  for  all  geologic  materials. 

Returning  again  to  Equations  4.15  and  4.21,  it  can  be  seen  that  both  h  and  (P^/Pg)  ^^re  quite  slowly 
varying  functions  of  x  when  x  exceeds  10,  since  7  in  a  solid  medium  is  generally  greater  than  2  and 
less  than  3.  Referring  to  Equation  4.24,  we  see  that  z  changes  very  slowly  when  x  exceeds  10,  for 
any  possible  value  of  b.  Therefore,  we  need  only  consider  how  those  quantities  vary  in  the  range  from 
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=  X.,  =  1/2  to  X  '  10.  The  same  is  also  true  of  G^,  given  by  Equation  4.22,  as  we  have  already 
shown  that  must  be  a  slowly  varying  function  of  x  when  x  exceeds  10;  then  P^^/Pg  approaches 

(y  -  l)/(7  +  1),  and  Taylor's  theory  becomes  applicable,  at  which  point  Qg/Q  approaches  zero. 


Referring  now  to  Equation  4.22,  we  recall  that  x  =  x„  at  the  upper  limit  of  the  plastic  zone.  It  is 
expected  that  the  heat  loss,  Q^,  should  be  approximately  equal  to  Q  for  any  value  of  x  less  than  x^j  in 
the  plastic  zone.  In  the  plastic  zone  where  x  <  x  <  x„.  Equation  4.22  is  simplified  by  this  assumption, 
giving  G^  =  1  +  X.  As  explained  previously,  the  general  plan  of  calculation  is  to  obtain  an  average  for 
h,  {p./p  ),  and  G  over  the  range  from  x  =  0.5  to  x  =  10.  It  is  not  clear  how  G  varies  when  x 

2  g  rrt  m 


m 


exceeds  x„.  Possibly  a  good  average  for  G  over  the  stated  range  is  the  value  corresponding  to  the 
n,  m 

point  where  x  =  x^^^.  Then  an  approximate  average  of  G^  to  be  used  in  Equation  4.26  is  =  1  +  x^j. 

In  principle,  one  should  now  take  the  average  of  h(Pj/p^)  over  the  entire  range  of  x  from  0.5  to  10. 
However,  it  was  not  possible  to  evaluate  G  for  any  value  of  x  greater  than  x„.  For  this  reason,  it  is 
expected  that  a  better  average  for  the  triple  product  might  result  from  emphasizing  the  upper 

end  of  the  range  of  x  near  x  =  10  when  an  average  for  the  double  product  (p^/p^)h  is  being  calculated. 
Since  p./p  is  a  decreasing  function  of  x  and  h  is  an  increasing  function  of  x,  the  product  varies  quite 
slowly  even  for  values  of  x  less  than  10.  Therefore,  it  is  sufficient  to  evaluate  (p./p  )h  at  x  =  10  in 
order  to  obtain  a  satisfactory  average  to  be  used  in  Equation  4.26.  From  Equations  4.15  and  4.21  we  find 
h,„  =  0.336,  and  (p./p  )  =  1  -  (20/ll)/(l  +  y)  .  Using  these  expressions  together  with  the  preceding 

equations,  ~  ^  ^a  °  obtain  a  general  expression  for  the  loss  term,  q^,  in 

Equations  4.25  and  4.26, 


q^  =  (0.4)z^  (1  +  Xjj)|^l  -  (20/ll)/(y+  l)j  .  (4.27) 

In  Equation  4.27,  z^  is  obtained  from  Figure  4.4  and  a  known  value  of  y  for  the  medium.  Values  of 

y  and  of  x„  for  the  medium  are  obtained  from  a  plot  of  x  versus  a  using  experimental  Hugoniot  data, 
li 

and  from  Equation  4.3,  As  in  Figure  4.2,  we  obtain  x„  from  the  point  at  which  the  plot  of  x  versus 

rl 

a  fir.st  departs  "significantly"  from  the  hydrodynamic  straight  line  where  y  is  constant.  Below  this 
point,  the  medium  is  plastic  and  y  is  not  constant.  The  exact  location  of  this  departure  is  more  or 
less  arbitrary.  As  a  rough  guide,  one  can  assume  that  the  departure  is  "significant "when  the  actual 
value  of  X  differs  by  about  30  percent  from  the  value  of  x  obtained  by  extending  the  hydrodynamic 
straight  line. 
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Variation  of  Pressure  with  Radial  Distance 
Hydrodynamic  Pressure  Range 

In  the  hydrodynamic  pressure  range  where  a  is  greater  than  p^,  the  value  of  x  at  any  distance 
R  is  found  by  integrating  Equation  4.19,  using  also  Equations  4.23  and  4.25, 


/ 

r.lr.  /  /  xl/3/1  +  b/x  \  3b 

where  R  =  R.  when  x  =  x^^,  and  from  Equation  4.3  we  obtain  the  pressure,  o/o.  =  x/x^  .  It  should 
be  obseiwed  that  y  is  assumed  to  be  constant. 


Crushing  Zone 

2 

When  the  wave-front  pressure,  a,  is  less  than  p  ,  we  take  x  =  1  -  (8/3){C./C,)  as  shown 

m  c 

by  Equation  4.9.  We  see  from  Figure  4.2  that,  in  general,  the  hydrodynamic  termination  pressure, 

P2  is  nearly  equal  to  p^.  Thus,  a  constant  value  of  x  determined  by  x^  from  Equation  4.9  can  be 

used  without  excessive  error  for  any  pressure  less  than  p„  and  greater  than  the  crushing  strength,  p  . 

A  c 

This  particular  modulus- correction  for  the  crushing  zone  is  merely  an  estimate,  made  necessary  by 
lack  of  experimental  data. 


Shock-wave  propagation  in  the  crushing  zone  is  not  described  by  Equation  4.28  since  y  in  the 
crushing  zone  varies  rapidly  with  pressure  a  .  In  the  crushing  zone  we  use  Equations  4.16  and  4,18, 
with  F  =  1, 


-  (d  log  Ug/d  log  r)  =  3  +  q  -  l/(l+x^)  ,  (4.29) 

in  which  q  =  G^(Pj/p„)zh(l  +  .  Here  the  heat  factor,  h,  is  x^/6  from  Equation  4.15,  using 

the  approximation  for  small  x  .  The  factor  G  is  obtained  from  Equation  4.22,  in  which  we  assume 

m 

that  Qg/Q  =  1.0  for  the  crushing  zone,  so  that  G^  =  1  +  x^.  This  combines  with  Equation  4.9  to  give 

G^  =  2  -  {8/3)(C_/C,)^  . 

m  21 


It  is  seen  from  Equations  4.21  and  4.9  that  Pj^/Pg  is  nearly  equal  to  1  in  the  crushing  zone.  Therefore, 


G  zx 
m  c 

6(1  +  X 
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and  Rquation  4.29  becomes 


-  (d  log  u^/d  log  r)  =  3  -  (1  -  zx  G  /6)/(l  +  x  )  . 

\  s  /  cm  c 

2  2  2  r  1 

Since  u”  =  <xe/p,  ,  it  follows  from  Equation  4.2  that  u  =  a  /  p,S, (1  +  x  )  and  Equation  4.30 

S  1  S[^1X  CJ 

becomes 

-d  logo/d  log  R  =  3/2  -  (1/2)(1  -  zx^G^/6>/(l  +  x^)  . 

Dividing  Equation  4.31  by  Equation  4.18.  with  subscript  c  in  crushing  zone, 

"(d  loga/d  log  z)  =  1  +  3x  /2  +  x  zG  /12  ,  ' 

c  c  m 

which  can  be  integrated  to  give 

2  +  3x 
_ c 

a/a.  =  (z./z)  ^  exp  f-(G  /12)x  (z  -  z.)]  . 

11  L  m  c  1  j 

In  Equation  4.32,  z  =  when  a  =  .  The  ratio  z/z.  is  obtained  by  integrating 

Equation  4.18  , 


1 


1  +  x„ 

z/z.  =  (R/R.) 


in  which  a  =  a.  and  z  =  z.  when  R  =  R.  . 

11  1 


(4.30) 


(4.31) 


(4.32) 


(4.33) 


Cracking  Zone  and  Reversible  Zone 

When  the  wave-front  overpressure,  a,  is  less  than  the  crushing  strength,  p^,  and  greater 
tlian  the  reversible  limit,  p^,  the  pressure  varies  inversely  as  the  square  of  the  radial  distance,  R, 
as  assumed  in  Section  2  for  the  cracking  zone.  Thus, 

a/a.  =  (R./R)^  ,  (4.34) 


in  which  a  =  a.  when  R  =  R. 

1  1 

When  a  is  less  than  p  ,  the  strains  are  reversible  and  q  =  0 
r 

Equation  4.25,  and  x  in  the  reversible  zone  is  given  by  Equation  4,5. 
from  Equation  4.28  with  b  =  2/3,  x  =  a/2pj,  so  that 


R/R^  =  (Oj/a) 


1/3 

\i+4p^7£ 


\2/3 


.  Thus,  b  =  2/3  according  to 
The  curve  of  a  versus  R  is  found 


.'X. 

(4.35) 
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when  o  is  less  than  p  .  fit  Equation  4.35,  p,  designates  the  ambient  pressure.  When  o  =  p  -  p.  is 
much  smaller  than  p^.  Equation  4.35  aj^ro^bhes  R  =  constant/o  ,  and  the  pressure,  a,  falls  off 
inversely  as  the  first  power  of  the  radius,  R.  This  range  corresponds  to  acoustic  overpressures,  for 
which  X  in  Equation  4.5  is  very  small  compared  to  1. 

Table  4.2  lists  the  stress-decay  formulae  associated  with  the  various  conditions  of  the  medium, 
and  the  corresponding  ranges  of  pressure. 


TABLE  4.2 

Stress  Decay  Laws  in  a  Solid  Medium 


Zone  Condition  of  medium  and  range  of  stress 

I  Continuous  solid,  reversible  strains, 

a  <  p 


Stress-decay  law 


/l  +  4p  /»  \ 
'l\  ■  (l44p;;3.j 


n  Cracked,  large  fragments  or  radial 

cracks,  p  <  o  <  p 
r  c 


a/o.  =  (R./R) 


III  Crushed,  small  particles, 

Pc<<^<P2 


o/o^  = 


(R./R) 


2  +  3x 
_ c 

2(1  +  X  ) 

c 


exp 


m 


1  +  X 


!(r/r.)  -  1 


IV  Fluid,  hydrodynamic  or  plastic, 

P2'  P2=  Pm 


j/o^  =  (R./R)'^ 


1  +  2p^c‘(b/o)/(T  +  1) 


1  +  2pjCj(b/crj)/(7  +  1) 


2  -  b 
b 


Initial  Conditions 

In  order  to  use  Equation  4.28  to  calculate  the  relation  between  a  and  R  in  the  hydrodynamic  zone, 

2 

one  must  first  obtain  the  "initial  value,"  x^,  at  some  given  radius,  R^,  Since  x^  =  (l/2)(7  +  Da^/pjC^ 
from  Equation  4.3,  it  Is  evidently  necessary  to  calculate  the  pressure,  o^,  at  some  given  initial  radius, 
Rj.  This  pressure  can  be  calculated  for  a  strong  shock  wave  (x  >>  1)  in  a  solid,  using  Taylor's 
formulae  given  above.  Strong  shock  wa'  .  x  >>  1 ,  are  always  produced  by  buried  nuclear 
explosions. 
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An  undergroond  chonical  ex|d<MEion  does  not  produce  n  strong  shock  wave,  md,  accordingly,  only 
a  fraction  d  the  incident  energjr.of  a  detonation  wave  is  transmitted  past  the  boundary  at  ttie  surrounding 
medium.  The  remaining  energy  is  reflected  back  into  the  detonation  products,  hi  order  to  calculate  the 
transmitted  energy,  one  must  first  obtain  the  energy  and  pressure  ot  the  incident  detonation  shock  wave. 

For  chemical  e]q>losions,  the  front  pressure,  of  a  detonation  shock  wave  is  obtained  from 

the  lUmkine-Hugoniot  mass  and  momentum  conservation  relations. 

u  =  AD  ,  (4.36) 

sD 


in  which  A  ==  (l  ~  P  /p  )  and  u  is  particle  velocity  behind  the  detonaticm  front.  The  ambient 
density  of  the  undetonated  explosive  is  designated  as  p^  and  the  detonation  front  velocity  is  D. 


Parameter  A  is  the  relative  fractional  volume  change  from  the  ambient  specific  volume,  p^ 


-1 


a  ,  to 


-1 


the  shock  f»*ont  specific  volume,  p^  .  The  factor  A  is  nearly  equal  to  0.25  for  several  types  of 


explosives  (Cook,  Keyes,  and  Ursenbach,  1961,  Lutzky,  1960), 


5,10 


In  order  to  find  the  pressure,  p^,  of  the  shock  wave  transmitted  into  the  surrounding  medium, 
one  must  first  obtain  the  reflection  pressure  as  a  function  of  the  reflection  particle  velocity  in  the 
explosive  detonation  products.  When  the  reflection  pressure,  p,  is  greater  than  the  incident  detonation 
front  pressure,  p^p,  the  reflection  particle  velocity,  u,  is  less  than  the  particle  velocity,  u^^,  of  the 
incident  detonation  wave,  and  the  curve  (p,  u)  is  described  as  a  "reflection  Hugoniot. "  When  the 
reflection  pressure,  p,  is  less  than  p^^,  the  reflection  particle  velocity,  u,  is  greater  than  u^^, 
and  the  curve  (p,  u)  is  termed  a  "reflection  isentrope. " 


These  reflection  curves  are  not  generally  known  from  direct  measurements,  but  formulae  for 

0 

calculating  an  explosive  reflection  Hugoniot  or  isentrope  have  been  given  by  W.  £.  Deal  (1958). 
Deal's  results  can  be  expressed  as  follows  in  our  notation.  When  the  reflection  pressure,  p,  is 
greater  than  the  incident  fi-ont  pressure,  p^^,  the  reflection  Hugoniot  is 


Ag  =  1  +  (1/4)(1  -  A)'Vl  -  u/u^jj)^  . 


Ag  -  1  -  (1/2  -  A)(l  -  A)'Vl  -  u/u^jj)^ 
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When  the  reflection  pressure,  p .  is  less  than  p^^^  the  reflection  isentr<^  is 


u/u  T.  =  1  + 
sD 


(4.30) 


0  =  (1/2  -  A)(l  -  A) 


-1 


The  transmitted  pressure,  p^.  and  the  transmitted  particle  velocity,  u^.  correspond  to  the  point 

of  intersection  (p  .  u  )  of  the  reflection  curve  (p.  u)  with  the  Hugoniot  curve  (a,  u  )  of  the  surrounding 
I  I  s 

medium.  The  Hugoniot  curve  (a.  u^)  of  the  surrounding  medium  is  generally  obtained  from  laboratoty 

experiments.  If  direct  measurements  are  not  available  for  the  medium  in  question,  one  can  use  a 

theoretical  Hugoniot.  derived  from  Equations  2.2  and  2.3  in  Section  2,  and  the  kinetic  energy  relation 

2 

at  the  front.  Ul2){ctlo^)  =  (l/2)Ug  .  which  is  independent  of  the  shock  strength.  It  is  found  that 


2  _  2o 

"s  ■  Pj(7  +  l)(l  +  l/x) 


(4.40) 


.  (7  +  1) g 
*  "  2 


(4.41) 


Here  is  the  longitudinal  wave  propagation  velocity  for  infinitesimal  strains  in  the  medium.  As 

already  mentioned  in  Section  2.  the  adiabatic  exponent  y  in  the  hydrodynamic  pressure  range  is  about 

2  for  hard  rock,  and  about  3  for  porous  wet  rock  and  porous  soils.  When  x  is  small  compared  to  1. 

Equations  4.40  and  4.41  reduce  to  the  correct  relation,  a  =  PjCjU^  .  for  small  acoustic  overpressures 

2 

(compared  to  the  modulus  regardless  of  the  particular  value  assumed  for  y  in  the  hydrodynamic 

pressure  range.  Therefore,  we  can  use  the  same  constant  value  of  7  in  Equations  4.40  and  4.41  over 

any  intermediate  pressure  range  with  reasonable  accuracy.  For  very  strong  shocks,  x  is  large  com- 

2 

pared  to  1,  and  Equation  4.40  reduces  to  the  simple  relation,  g  =  (l/2)(7  +  l)p,u  . 

z  s 

Shell  Thickness  Ratio 

In  order  to  calculate  the  shell  thickness  of  the  spherical  shock  wave  transmitted  into  the  surround¬ 
ing  medium  from  an  incident  explosive  detonation  wave,  it  is  first  necessary  to  find  the  shell  thickness, 
y^,  of  the  explosive  detcmation  wave,  using  Equation  4.11.  This  gives  the  (constant)  shell  thickness 
ratio,  z  »  y./K.  <  at  any  spherical  detonation  front  radius,  R  ,  in  th :i  explosive,  using  known 
relations  between  ipjpj,  (u_/u,),  and  n  «  r/R  for  an  explosive  which  has  a  given  undetonated  density,  p  , 
These  relations  are  termed  "radial  profiles."  Such  profiles  have  been  calculated  by  Lutzky  (1960)^^ 
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for  several  explosives,  including  TNT  at  a  density  p  =  1.58  gm/cm  .  Radial  profiles  for  pentolite. 

^  .  © 

RDX,  tetryl,  and  RDX/TNT  have  also  been  calculated. 

~  R  is  sometimes  necessary  to  obtain  the  shell  thickness  ratio,  z  ^  =  y  /R  .  for  an  explosive 
vliose  ambient  undetonated  density,  p^.  is  not  the  same  as  the  particular  density,  p^.  for  which  radial 
profiles  have  been  calculated.  An  approximate  method  for  finding  z  (at  density  p  )  from  a  known  value, 

&  Sl 

2 

z^,  calculated  at  density  p^,  is  given  by  Bishop  (1962a). 

z  Iz  =  (D  /D  .  (4.42) 

a  e  a  e 

Here  D  represents  the  detonation  front  velocity  which  corresponds  to  an  explosive  density,  p  ,  and 
d  & 

is  the  velocity  which  corresponds  to  a  density,  p^ 

In  order  to  determine  the  shell  thickness  ratio,  ^  shock  wave  transmitted  into 

the  surrounding  medium,  compared  to  the  shell  thickness  ratio,  z  ^  =  y  /R  ,  of  the  incident  detonation 

&  &  SI 

wave  at  the  medium  boundary,  we  observe  that,  in  a  time  during  which  the  detonation  front  moves 

a  distance  y^,  the  wave  front  in  the  medium  moves  a  distance  y^,  in  a  time  where  is  the 

initial  shock  front  propagation  velocity  in  the  medium  near  the  boundary.  Since  the  detonation  shock 

wave  can  supply  energy  "effectively*  only  during  the  time  the  effective  shell  thickness,  y^,  of 

the  shock  wave  in  the  medium  near  the  boundary  is  the  same  as  the  distance  that  its  front  moves  in  the 

same  time.  Thus,  y /D  =  y/U^,  or  y/y  =  U./D  .  In  terms  of  the  shell  thickness  ratio  z,  this 
a  a  t  t  t  a  t  a 

relation  is  equivalent  to  z  =  (D  /U  )(z  )(R./R  ).  Since  R  is  essentially  equal  to  R  near  the 
explosive- medium  boundary,  we  have,  finally, 

z,/z„  =  D /U,  .  (4.43) 

t  a  at 

The  shock  front  propagation  velocity,  U^,  in  the  medium  can  be  found  from 

since  p^  and  u^  have  already  been  determined. 

It  is  sometimes  necessary  to  adjust  the  transmitted  pressure,  p^,  and  the  transmitted  shell 
thickness  ratio,  z^ .  If  p^  is  equal  to  or  greater  than  the  plastic  flow  upper  limit,  p^,  for  the  medium, 
the  pressure  is  said  to  be  "hydrodynamic"  (refer  to  Figure  4,2).  In  the  hydrodynamic  pressure  range, 
the  hydrodynamioally  stable  shell  thickness  ratio,  Klven  by  the  graph  in  Figure  4.4 

as  a  function  of  the  acRabattc  exponent,  y  .  R  p^  is  greater  than  and  z^  la  greater  than  z^,  there 
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is  an  unstable  balance  between  the  pressure,  and  the  shell  thickness  ratio,  .  The  shell  thickness 
ratio  must  rapidly  shift  to  the  stable  value,  z  ,  and  the  front  pressure,  p  ,  simultaneously  shifts  to  a 

O  X 

corresponding  stable  value,  p  .  The  shift  is  completed  at  a  minimum  stable  radius,  R^,  given  (Bishop, 

2 

1962a)  in  terms  of  the  initial  radius,  R^,  at  the  explosive-medium  boundary, 

R  /K  =  (1  -  2/z  )'^  .  (4.45) 

o  t  o 

The  stable  pressure,  p  ,  is  found  in  terms  of  the  transmitted  pressure,  p,,  the  transmitted  modulus 
o  t 

correction,  x^,  and  the  stable  modulus  correction,  x^, 

p^  =  p^(x^/x^)  .  (4.46) 

In  Equation  4.46,  the  stable  modulus  correction,  x  ,  is  obtained  by  solving  the  following  quadratic 

2 

equation,  (Bishop.  1962a)  ; 


X  /x 
o 


t 


od  +  Z/x  )  -  1 
o 

o(l  +  l/x^)  -  1 


in  which  a  *  {ll2)(y  +  1)  and  x^  =  (U^/Cj)^  -  1 


(4.47) 


Comparison  of  Calculations  With  Experiments 

* 

Numerous  measurements  have  been  made  of  particle  velocities  resulting  from  explosions  in  various 
geologic  media.  Also,  from  new  developments  in  Instrumentation,  a  considerable  number  of  direct 
measurements  of  pressure  have  recently  been  obtained,  principally  in  the  hydrodynamic  zone.  Peak 
particle  velocity  data  may  be  converted  to  peak  wave  front  pressures  by  the  momentum  relation, 
o  =  PjUu^,  if  the  density,  p^,  and  wave  front  velocity,  U,  are  known. 

The  numerous  data  which  have  been  obtained  over  the  past  few  years  have  largely  defied  accurate 
theoretical  description.  Only  within  the  last  two  years  has  significant  progress  been  made  in  attempts 
to  calculate  peak  wave  pressure  as  a  function  of  distance  from  explosions  in  geologic  media.  Needless 
to  say.  very  little  progress  has  bem  made  in  calculating  the  particle  velocity  waveforms  which  are 
observed.  The  foremost  difficulty  in  the  theoretical  calculation  has  been  and  crntinues  to  be  our  great 
lack  of  knowledge  about  the  properties  of  these  media,  without  which  meaningful  and  accurate  calculations 
•  cannot  be  performed. 
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Extensive  data  have  been  obtained  from  both  chemical  and  nuclear  explosions  buried  in  four  media: 
granite,  halite,  volcanic  tuff,  and  desert  alluvium.  These  data  have  been  plotted  in  Figures  4.7  to  4.11, 
together  with  curves  calctilated  from  equations  derived  in  the  preceding  sections.  Figure  4.12  compares 
pressure-distance  curves  resulting  from  nuclear  explosions  in  various  media. 

Hugoniot  data  used  in  the  calculations  for  the  four  media  are  shown  in  Figures  4.13  through  4.16. 

In  these  figures,  the  modulus  correction,  x,  is  obtained  by  substituting  experimental  Hugoniot  data  into 
Equation  4.4.  Other  properties  of  the  media  required  in  these  calculations  are  shown  in  Table  4.3. 

Note  that  in  Table  4.3  the  infinitesimal  dilatational  velocities,  C^,  are  obtained  from  the  measured 

1  /2 

wave  velocities,  U,  using  Cj  =  U/(l  +  x^)  ,  from  Equation  4.4.  As  n»ay  be  observed  from  the  data 

of  Figures  4.7  through  4.11,  very  few,  if  any,  measurements  were  made  in  the  acoustic  regions  around 
the  explosions.  In  an  acoustic  region,  pressure  or  peak  particle  velocity  varies  inversely  as  the  distance 
from  the  explosion,  and  the  acoustic  wave  front  propagates  with  the  velocity.  C^.  Since  almost  all  of 
the  measurements  of  wave  arrival  were  made  in  the  crush-crack  zone  where  peak  particle  velocities 
decrease  with  distance  according  to  a  negative  power  greater  than  1,  we  interpret  the  measured  velocities 
as  being  those  of  finite  amplitude  waves,  and  determine  from  them  by  the  above  formula. 


TABLE  4.3 
Medium  Properties 


Mediiun 

type 

Explosion 
aerie*  or 
event 

Average 

density 

(gm/cm^) 

Measured 
wave  velocity 

U 

(km/sec) 

Infinitesimal 
wave  velocity 

(km/sec) 

X 

c 

_JL 

J!iL 

Minimum  hydro- 
dynamic  pre.^sure 

ph 

(kllobars) 

Granite 

Hardhat 

3.67 

5.5^ 

5.2 

0.111 

3.18 

0.72 

530 

Halite 

Cowboy 

2.15 

4.4* 

4.3 

0.094 

3.16 

0.23 

90 

Halite 

Gnome 

2.15 

4.8^ 

4.6 

0.094 

2.16 

0.23 

108 

Tuff 

Hobo 

1.85 

2.44< 

2.07 

0.385 

3.0 

3.5 

125 

Tuff 

Rainier 

1.85 

2.44® 

3.07 

0.385 

3.0 

3.5 

135 

Alluvium 

Scooter 

1.6 

1.20® 

1.04 

0.333 

3.0 

14 

135 

Alluvium 

Teapot  Ess 

1.6 

1.07^ 

0.93 

0.333 

3.0 

30 

135 

1.  Ferret,  1963a:  Chabai  and  Bass.  1963 

2.  Murphey,  1960,  Nicholls,  Hookei,  and  Duvall.  1960 

3.  Weart,  1960 

4.  Young,  1961 

5.  Ferret,  1961 

6.  Ferret,  1963b 

7.  Sacha  and  Swift,  1955 


CALCULATED  CURVES 


Lce  for  PeUetol  Figure  4.8  Pressure  Versus  Distance  for  a  Nuclear 

Explosion  in  Tuff 
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EXPCmMCNTAL 


To  illustrate  the  method  of  application  of  the  formulas  summarized  in  Table  4.2,  we  describe  briefly 
the  pi'oeedure  followed  for  obtaining  the  calculated  curves  of  Figures  4.7  to  4.11.  For  chemical  explosions, 
values  of  the  transmitted  quantities  p^,  u^,  and  in  the  medium  at  the  explosive  boimdary,  R  =  R^,  are 
determined  by  the  point  of  intersection  of  the  medium  Hugoniot  (a,  u^)  with  the  reflection  curve  (p,  u)  of 
the  explosive,  as  described  under  "Initial  Conditions"  and  shown  in  Figure  4.17  for  desert  alluvium.  The 


Figure  4.17  Desert  Alluvium  Hugoniot  and  TNT  Reflection  Isentrope 
(p^  =  1,375  gm/cm^) 
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shell  thickness  ratio,  z^,  of  the  transmitted  wave  at  the  boundary  is  determined  by  Equations  4.43,  4.42, 

and  4.11,  in  conjunction  with  the  radial  velocity  and  density  profiles  of  the  explosive,  as  explained  lander 

"Initial  Conditions."  The  initial  value  of  the  pressure,  in  Equation  4.28  is  p^;  similarly,  we  have 

X.  =  x^,  and  the  initial  radius  is  R^.  If  p^^  >  p^  >  Pg#  o  is  calculated  as  a  function  of  R  by 

Equation  4.28,  as  indicated  in  Table  4.2,  to  the  distance  Rj  at  which  a  =  p^  and  x  =  x^  =  1/2.  The 

2 

value  of  p-  is  obtained  directly  from  p.  =  p,(x„/x,),  where  x  =  (U,/C,)  -  1. 

4  ttl 

The  value  of  b  must  also  be  computed  for  use  in  Equation  4.28.  By  Equations  4.25  and  4.27,  b 
is  readily  evaluated  in  terms  of  the  medium  constants  y  and  x„,  since  the  shell  thickness  ratio,  1/z  , 

H  O 

equal  to  (y/R)^,  is  given  as  a  function  of  y  in  Figure  4.4. 

The  initial  value  of  the  shell  thickness  ratio,  z.,  in  Equation  4.24,  is  equal  to  the  transmitted 
shell  thickness  ratio,  z^,  at  the  boundary.  Since  x.  =  x^,  and  Xg  =  1/2,  we  can  evaluate  the  shell 
thickness  ratio,  •  from  Equation  4.24,  at  the  hydrodynamic  termination  point  where  the  pressure,  cr, 
is  equal  to  P2.  In  the  range  from  a  =  P2  to  a  =  p^,  the  pressure  as  a  function  of  distance  is  determined 
by  Equations  4.32  and  4.33,  in  which  the  initial  pressure,  p.,  is  equal  to  p„,  and  the  parameter  z. 
is  equal  to  z„.  In  the  crack  zone,  p  >  a  >  p  ,  pressure  decays  inversely  as  the  square  of  distance 

A  C  IT 

(Equation  4,34),  and  in  the  reversible  zone,  a  <  p^,  the  pressure  falloff  is  by  Equation  4.35. 

Since  there  are  no  data  available  regarding  the  values  of  p^  for  geologic  materials,  we  must 
estimate  p^  from  wave  pressure  and  particle  velocity  data  of  explosion  experiments.  In  Figures  4.7  to 
4.11,  we  note  the  minimum  pressure  level  at  which  o  varies  inversely  with  the  square  of  R;  for  pressures 
smaller  than  this  value  (which  we  take  to  be  p^),  the  falloff  rapidly  approaches  an  inverse  R  dependence. 
Values  of  p^  estimated  in  this  manner  are  listed  in  Table  4.4  . 

From  explosion  experiments  we  can  also  obtain  a  value  for  the  dynamic  compressive  strength,  p^, 

of  an  earth  material  by  observing  the  maximum  pressure  at  which  o  varies  inversely  as  the  square  of  R. 

Fortunately,  laboratory  determinations  of  the  static  value  of  p^  provide  a  helpful  guide  (Duvall  and 

7  23 

Atchison,  1957;  Warner  and  Violet,  1959).  ’  As  has  been  mentioned,  laboratory  tests  reveal  that  the 

dynamic  strength,  p^,  is  about  twice  the  static  value  for  several  materials  (Wuerker,  1959;  Watstein, 

26  24 

1953).  ’  Whether  or  not  these  results  from  tests  on  laboratory  samples  can  be  applied  to  gross 

regions  of  material  in  situ  is  unknown.  Nevertheless,  such  determinations  of  even  static  compressive 
strengths  give  valuable  information  for  order  of  magnitude  comparison. 
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TABLE  4.4 


Summary  of  Explosions  and  Corresponding  Numerical  Quantities  Used  in  Calculations 


TNT 

Ambient 

Calculated  Initial  values 

Hydrodynamic 
termination  values 
(x2  -  1/3) 

Deduced 

medium 

Medium  and 
type  of 
explosion 

radius 

(feet) 

Energy 

release 

Wt 

pressure 
Pi  =  Pghs 
(bars) 

<T. 

1 

(kilobars) 

P2  ^  ^ 

(kilobars)  ^t  *2 

^i 

(feet) 

H 

P  P 

(bars) 

Granite, 

nuclear, 

Hardhat 

33,7 

5  kt 

80 

45,400 

4.0 

10.5 

100 

227  0.75  15,6 

4000 

250 

Halite, 

TNT  (pelletol). 
Cowboy 

1  /I 

O.ISCW^^''’ 

Various, 
20  lb  to 
1000  lb 

60 

87 

TNT 

radius 

Rt 

10.4 

0.18 

Nonhydrodynamic 
initial  pressure 

2500 

10 

Halite, 

nuclear. 

Gnome 

28.4 

3  kt 

76 

28,900 

3.88 

10.5 

100 

144.5  0.97  16.3 

2500 

10 

Tuff,  240  ft., 
TNT  (pelletol), 
Hobo 

1,54 

968  lb 

11.4 

76 

1.54 

12.1 

2.57 

14.8  1.47  14,5 

1400 

5 

Tuff,  1000  ft,, 
TNT  (pelletol). 
Hobo 

1/3 

O.lSGWj.  ' 

103,  516, 
and 

973  lb 

64 

76 

TNT 

radius 

Rt 

12.1 

2.57 

14.8  1.47  14.5 

1400 

80  to  200 

Tuff. 

nuclear. 

Rainier 

23.5 

1.7  kt 

25  to  60 

3,960 

7.18 

8.0 

100 

14.8  1.38  10.5 

650 

Desert 

Alluvium,  TNT, 
Scooter 

14 

10®  lb 

7.4 

(131) 

29 

(14.0) 

18,7 

(22) 

8.0 

(15) 

3.4 

4.3  1.61  8.7 

6 

1 

Desert 

alluvium, 

nuclear. 

Teapot  ESS 

18.75 

1.2  kt 

4.4 

680 

11.5 

8.0 

100 

3.42  1.46  8.9 

It  should  be  noted  that  quantities  p^  and  p^  listed  in  Table  4.4  for  desert  alluvium  are  estimated 

from  the  Jangle  HE-2  data  points  plotted  in  Figure  4.11.  The  dynamic  compressive  strength  of  halite  has 

12 

been  measured  directly  by  Nicholls,  Hooker,  and  Duvall  (1960). 

With  regard  to  the  determination  of  in  Equation  4.24,  this  is  generally  equal  to  z^  (Equation  4.43) 
for  most  chemical  underground  explosions,  but  it  is  possible  for  z^  to  exceed  z^  with  certain  combinations 
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of  explosives  and  media.  Here,  1/z^  is  the  hydrodynamically  stahle  Aell  thickness  ratio  (y/R)^  obtained 
from  Figure  4.4.  If.  moreover,  the  initially  transmitted  pressure,  p^,  exceeds  the  hydixxilynamic  limit, 
as  noted  under  "Initial  Conditims,"  is  then  an  unstable  shell  thickness  ratio  and  rapidly  shifts  to 
the  stable  value,  z^.  Simultaneously,  the  pressure,  p^,  shifts  to  the  stable  value  p^  determined  from 
Equations  4.45,  4.46,  and  4.47.  In  Table  4.4  the  initially  transmitted  values  in  parentheses  for  the  Scooter 
explosion  represent  unstable  quantities  in  desert  alluvium. 

For  nuclear  explosions  where  the  starting  pressures  are  extremely  high  and  the  initial  disturbance 

is  fully  hydrodynamic  (o  >>  p„).  we  take  any  value  for  x.  which  is  greater  than  10,  for  example,  x,  =  100. 

n  &  1 

2 

Then  is  given  by  Oj  =  2pjCjXj/(y  +  1),  from  Equation  4.3.  Using  known  values  of  o..  and  y,  the 
corresponding  initial  radius,  R^,  is  obtained  from  Equation  4.13.  Now  in  the 

hydrodynamic  range,  a^>a>  p^,  the  pressure,  a.  varies  with  distance.  R.  according  to  Equation  4.28, 
and  the  subsequent  calculations!  procedure  is  exactly  as  described  above  for  chemical  explosions.  Note 
that  in  the  present  context,  =  z^,  and  z^^  is  given  in  Figure  4.4.  For  complete  details  of  each  of  the 
explosion  experiments  listed  in  Table  4.4,  reference  may  be  made  to  two  research  reports  by  the  author 
(Bishop.  1962a,  1962b).^ 

The  energy  release,  W^,  of  the  Hardhat  nuclear  explosion  in  granite  is  listed  as  5  kt  in  Table  4.4. 

It  should  be  pointed  out  that  there  is  an  uncertainty  in  this  value  (obtained  from  rudiochcmistry)  of  about 
30  percent.  The  error  associated  with  the  3  kt  value  for  Gnome  is  in  excess  of  30  percent,  but  the  errors 
for  the  remaining  nuclear  explosions  in  Table  4.4  are  believed  to  be  less  than  20  percent. 

Also  in  Table  4.4,  it  is  to  be  noted  that  the  initially  transmitted  pressures  ({^  =  (^)  from  the  pelletol 
(TNT,  ■  1)  chemical  explosions  in  halite  are  less  than  100  kilobars,  and  consequently  are  definitely 
not  hydrodynamic  (see  Figures  4.9  and  4.14). 

The  Initial  shell  thickness  ratios  (y^/R^  =  l/z^)  in  Table  4.4  are  of  the  order  of  0.1  for  all  the 
explosions  listed,  indicating  that  the  effective  length,  y^,  of  the  wave  transmitted  into  each  medium  is  of 
the  order  of  magnitude  of  10  percent  of  the  starting  radius,  R^,  . 

Values  of  Pj,  listed  in  Table  4.4,  refer  to  the  lithostatic  pressure  at  the  depth  of  the  associated 
instrumentation  from  which  data  were  obtained.  For  the  Rainier  explosion,  p^  varies  from  25  to  60  bars, 
corresponding  to  the  depths  of  the  various  velocity  gages  used  in  that  experiment  (Ferret,  1961).^^  As  a 
result,  one  might  expect  the  values  of  p^  and  p^  deduced  from  Rainier  to  differ  somewhat  from  those  for 
the  Hobo  explosions,  since  these  quantities  may  depend  upon  depth. 
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Note  that,  for  a  nuclear  eqdoaian,  flie  parameter  in  TaUe  4.4  refora  to  tiie  radiua  of  w 
"equivalent"  q>liere  of  TNT  that  would  releaae  the  name  ti^  eneror.  The  densl^.  p^,  at  the  equiva- 

3 

lent  TNT  qphere  ia  assumed  to  be  1.0  gm/cm  for  Gnome.  Hardhat,  mtd  Rainier,  wfaidi  are  compared  to 
TNT  (pelletol^  explosions  in  correqionfiing  media,  llier^ore. 

R^ffeet)  =  {19.7Hkiloton8)*^’  . 

3 

The  density,  P^,  is  taken  as  1.37S  gm/cm  for  the  Teapot  Ess  and  Nougat  nuclear  e^losions  in  alluvium, 
since  these  are  conqmred  with  the  Scooter  explosion  (Figures  4.11.  4.17).  for  iHiich  R^  =  14  feet. 

There  is  always  a  hy<hx>dynamic  xone  surrounding  aqy  tamped  nuclear  explosion  that  releases  more 
energy  per  unit  volume  than  TNT.  It  was  shown  from  Taylor's  strong-shock  theory  (Equation  4.13)  that 
initially  R.  =  ^B(Tf)W^/oJ  with  good  accuracy,  provided  x.  is  greater  than  1®.  Therefore,  the  particu¬ 
lar  distance.  R^,  at  vdiich  any  high  pressure,  a.,  occurs  (x^  >  10)  will  be  proportional  to  the  cube  root  of 

1/3 

the  energy.  .  It  is  then  evident  that  the  radial  distance.  R.  will  be  pnnwrtional  to  at  any  given 

pressure,  a.  in  the  crushing  xone.  the  crack  xone.  or  the  reversible  zone.  Cube-root  scaling  in  each 

low-pressure  zmie  results  from  the  fact  that  the  termination  radius.  R2,  at  the  beginning  of  the  crushing 

1/3 

zone  is  proportimial  to  .  This  last  proportionality  arises  from  the  fact  that  o  is  a  function  only  of 

(R/R^)  in  the  hydrodynamic  zone  of  a  given  medium  (see  Table  4.2).  and  R.,  in  turn,  is  proportional  to 


In  Figure  4. 7  are  shown  data  obtained  from  experiments  with  pelletol  (TNT)  explosions  in  volcanic 

tuff  (Hobo).  One  series  of  experiments  was  conducted  at  240-foot  depth  and  another  at  1000-foot  depth. 

Scatter  in  the  data  is  very  great.  The  source  of  scatter  is  unknown;  it  may  be  due  to  uncertainties  in  the 

values  of  W^,  gross  medium  inhomogeneities,  large  experimental  errors,  or  a  combination  of  these.  In 

order  to  give  an  idea  of  the  effect  of  an  inhomogeneity,  a  curve  has  been  drawn  in  Figure  4. 7  representing 

calculated  pressures  in  a  homogeneous  medium  with  acoustic  velocity  =  2.8  km/sec  and  density 
3 

-  2.2  gm/cm  .  Although  such  a  hypothetical  medium  does  not  represent  actual  tuff  characteristics,  it  is 
interesting  to  observe  that  only  a  slight  change  in  the  sonic  velocity,  results  in  a  considerable  change 
in  the  transmitted  pressures. 

Data  from  the  nuclear  exploslm  in  tuff  (Rainier)  are  seen  in  Figure  4.S  together  with  the  calculated 

13 

pressure- distance  curve.  Also  shown  is  a  curve  for  Rainier  calculated  by  Nuckolls  (1959),  and  data  of 
Figure  4. 7  from  the  Hobo  experimoits  at  240  feet. 

Figure  4.9  illustrates  data  from  vaulous  escplosions  in  salt  media  (halite).  The  data  are  very  good, 
and  are  well  described  by  the  calculated  curves  within  the  limits  of  experimental  error.  The  dotted  curve 
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is  from  a  computer  calculation  by  Seidl  (Letter  from  Lawrence  Radiation  Laboratory  to  Sandia  Corporation. 
Albuquerque.  New  Mexico.  Octc^r  10.  1961)  for  a  nuclear  explosion  in  halite. 

Data  from  the  5  kt  nuclear  eiqilosioa  in  granite  (Hardhat)  are  plotted  in  Figure  4.10.  Except  for  the 
two  valuable  and  independent  measurements  of  shock  pressure  near  500  kilobars.  the  remaining  data  are 
all  at  pressures  less  than  4  kilobars.  R  is  seen  that  not  enough  measurements  were  made  in  the  elastic 
region;  consequently,  it  is  difficult  to  estimate  the  reversible  limit  stress,  p^.  from  the  pressure  data  of 
this  figure.  The  com.<uter  curves  calculated  by  Seidl  (Lav/rence  Radiation  Laboratory  Memorandum. 
Livermore.  California.  UOPKA  62-6.  January  19.  1962)  agree  nicely  with  the  two  data  points  near  500 
kilobars.  but  below  10  kilobars  Seidl's  curves  lie  definitely  above  the  experimental  data  points. 

Numerous  data  from  explosions  in  desert  alluvium  (Figure  4.11)  cover  almost  completely  the  range 
of  pressures  from  500  kilobars  to  0.1  bar.  Most  of  the  data  are  results  of  velocity  measurements.  b>  the 
alluvium  hydrtx^amic  region,  shock  pressures  were  measured  directly  in  the  Nougat  series  of  nuclear 
explosions  using  a  new  "impedance-mismatch”  gage  described  in  Section  1  of  this  report.  Since  several 
of  the  jnelds.  W^,  of  the  Nougat  explosions  are  classified,  only  the  scaled  distances  (R/R^)  can  be 
presented  with  the  measured  pressures.  Considering  the  probability  that  the  desert  alluvium  medium  was 
not  identical  for  all  the  explosions  of  this  series,  it  is  believed  that  agreement  between  calculations  and 
experimental  data  is  excellent. 

The  theoretical  curve  for  the  Teapot  Ess  alluvium  (C^  =  0.92  km/sec)  has  been  plotted  in  Figure  4.11 
as  a  solid  line.  It  is  seen  that  the  Teapot  Ess  experimental  points  agree  quite  well  with  this  curve.  The 
same  is  generally  true  of  the  Nougat  experimental  points,  but  the  scatter  is  much  greater  than  that  of  the 
Teapot  Ess  points. 

It  Is  evident  that  all  but  one  of  the  Nougat  points,  ad  well  as  all  of  the  Teapot  Ess  points,  are  below 
the  corresponding  points  of  the  Scooter  explosion  (TNT).  This  is  probably  due  to  the  exceptionally  high 
value  of  the  Scooter  alluvium  dilatational  velocity  (C^  =  1.04  km/sec)  rather  than  to  any  great  difference 
between  nuclear  and  chemical  explosions  in  the  same  type  of  alluvium. 

It  must  be  emphasized  that  the  discontinuity  in  the  slope  of  each  calculated  curve  (Figure  4.12)  at  the 
hydrodynamic  termination  pressure,  p^.  is  not  realistic,  but  is  a  consequence  of  the  assumption  that  y 
remains  constant  throughout  the  hydrodynamic  and  plastic  zones  down  to  the  termination  pressure,  p^,  as 
shown  by  the  dashed  4S-degree  straight  line  in  Figure  4.2.  Experimental  Hugoniot  data  (Figures  4.13  to 
4.16)  reveal  that  y  is  not  strictly  constant  except  for  a  range  of  pressure  (o>  p^^)  which  is  much  greater 
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than  Pg.  In  the  region  of  pressure  near  P2.  Y  has  an  average  value  which  differs  from  that  of  the  higher 
pressure  regirai  which  was  used  in  the  calculations.  The  effect  of  including  a  variable  Y  in  the  calcula¬ 
tions  would  be  to  remove  the  discontinuity  and  to  produce  a  smooth  pressure-distance  curve  which  would 
lie  above  each  one  shown  in  Figures  4.7  to  4.11,  in  the  vicinity  of  the  pressure  P2«  The  difference  in 
pressures  would  be  comparable  to  that  between  the  experimental  Hugoniot  curve  {Figures  4.13  to  4.16) 
and  the  projection  of  the  4S-degree  straight  line  from  the  hydrodynamic  zone. 

Discussion  and  Conclusions 


Review  of  Assumptions 

Calculations  have  been  made,  using  a  constant  parameter  y  in  the  plastic  zone,  where  the  relation 
between  log  a  and  log  x  is  actually  a  curve.  Let  us  examine  the  resulting  error  in  halite.  Referring  to 
Figure  4.14,  we  obtain  the  "constant  y"  termination  pressure,  Pj.  from  the  intersection  of  the  x  =  l/2 
line  with  the  45-degree  hydrodynamic  straight  line,  which  gives  132  kilobars.  The  "exact"  value  of  Pj 
is  obtained  from  the  cuiwed  graph  at  x  »  1/2,  which  gives  139  kilobars. 

Another  assumption  of  consequence  is  the  use  of  f  =  2  in  Equation  4.8,  The  modulus  ratio, 
f  =  S  /S  ,  is  defined  in  terms  of  the  "rapid- strain"  shock  modulus,  S  ,  and  the  "slow-straln"  acoustic 

SI  3 

modulus,  Sj.  The  modulus  ratio,  f,  is  a  useful  parameter  only  in  the  crush-crack  zones,  where  it  is 

nearly  constant.  In  the  reversible  zone,  f  =  1  +  x,  which  is  not  constant;  nevertheless,  it  is  possible 
to  estimate  the  upper  limit  of  f  for  wave  front  pressures  near  the  reversible  limit,  where  a  =  p^. 

Since  p^,  <  o  <  p^,  in  xhe  crack  zone  (Figure  4,2),  the  upper  limit  of  f  in  the  reversible  zone  should  pro¬ 
vide  sn  estimate  for  the  nearly  constant  value  of  f  in  the  crack  zone  (and  also  in  the  crush  zone).  As 
shown  following  Equation  4.5,  x  may  have  any  value  less  than  or  equal  to  1  +  T/p^  in  the  reversible 
zone.  Here,  p^  is  the  ambient  pressure,  and  T  is  the  dynamic  tensile  strength.  Therefore,  f  may 
have  a  maximum  value  of  2  +  T/p^  in  the  reversible  zone.  The  value  of  f  in  the  crack  zone  is  probably 
comparable  to  2,  since  the  term  T/p^  does  not  apply  to  the  crack  zone.  The  constant  value  f  ==  2  is 
consistent  with  experimental  data  for  several  types  of  solid  media  (Wuerker,  1959;  Watsteln,  1953).  ' 

We  have  also  used  the  constant  value  <>  3  in  the  crushing  zone  (p^  <  a  <  P2)  when  making 

calculations.  As  shown  in  the  derivation  of  Equation  4.30,  the  relation  G  >  1  4  x  leads  to 

me 

2 

^m  '  ^  in  the  crushing  zone.  It  has  been  found  that  ■  2  gives  a  better  fit  to  the 

experimental  pressure  distance  points  (Figures  4,7  to  4.11),  In  the  derivation  of  Equation  4,30,  it  should 
be  noted  that  the  heat  loss  expression  h  ■  Q/u  ■  x  /6,  Is  only  approximate.  It  is  well  known  that 

8  C 
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the  energy,  Q,  is  not  the  exact  heat  loss  because  the  Hugoniot  e^qiansion  curve  (Figure  4.6)  only  approxi¬ 
mates  the  true  adiabat  curve.  Therefore,  the  better  numerical  value  of  2  for  may  be  interpreted  as 
an  empirical  factor  wliich  adjusts  the  product  hG^  in  Equation  4. 1 7  so  that  it  has  a  more  correct  average 
value  over  the  aitire  range  of  pressures  in  the  crushing  zone. 


hi  the  hydrodynamic  and  plastic  zones,  we  have  used  the  relation  =  1  +  Xjj  in  order  to 

evaluate  Equation  4.26.  Since  the  exact  value  of  Q^/Q  fn  Equation  4.22  is  unknown  when  x  is  greater 
than  Xjj,  we  have  used  the  equality  =  Q  only  when  x  is  less  than  Xjj. 

The  parameter  x„  is  evaluated  from  the  curve  of  x  versus  o  (Figures  4.13  to  4.16).  It  is  taken  as 
n 

the  point  where  x  on  the  curved  graph  is  about  30  percent  different  from  the  corresponding  intercept 
(a  =  constant)  on  the  projected  hydrodynamic  45-degree  straight  line.  This  point  is  defined  physically  in 
terms  of  the  minimum  hydrodynamic  pressure,  p^, 

Xjj  =  (1/2)(7+  l)Pjj/(PjCj)  . 


The  medium  is  said  to  be  plastic  at  lower  pressures  where  Pj  <  o  <  Pjj.  If  an  experimental  Hugoniot 
curve  is  not  available,  one  can  determine  x„  from  the  above  expression  if  the  minimum  hydrodynamic 
pressure,  p„,  can  be  estimated  --  for  example,  by  Interpolation  from  Table  4.3.  It  is  also  necessary 
to  have  a  numerical  value  for  the  acoustic  dilatational  velocity,  C^.  It  is  evident  from  Table  4.3  that  y 
is  about  3  for  porous  rock,  and  about  2  for  nonporous  rock.  The  acoustic  dilatational  velocity,  C^,  also 
determines  the  "termination-pressure"  p^,  defined  by  the  condition  x^  *  1/2,  so  that  p^  =  PjCj/(7+  1). 


It  must  be  noted  that  the  dilatational  acoustic  velocity,  Cj,  is  not  generally  identical  to  the  measured 

-1/2 

wave  front  velocity,  U,  but  is  calculated  from  the  equation  Cj  «  U(1  +  x)  .  The  modulus  correction, 

2 

X,  is  equal  to  x  =  1  -  (8/3)(C-/C,)  in  the  crush- crack  zones,  and  x  °  (7/2p.  in  the  reversible  zone. 

C  4  1  1 

Thus,  X  is  not  usually  negligible  compared  to  1,  and  the  wave  front  velocity,  U,  corresponds  to  that  of  a 
finite- amplitude  disturbance. 


Conclusions 

The  properties  of  a  solid  medium  which  determine  the  manner  of  spherical  wave  propagation  from 

buried  explosions  have  been  described  by  means  of  a  correction  factor  applied  to  the  known  elastic 
2 

modulus,  of  the  medium.  Presumably,  elastic  properties  are  well  known.  Also,  information  is 

available  (Hugoniot  data)  concerning  the  hydro<(ynamic  behavior  of  solid  media.  However,  at  the  inter¬ 
mediate  stress  levels  where  finite  amplitude  stress  waves  produce  such  effects  as  plastic  flow,  crushing, 
and  cracking,  no  experimental  Information  Is  available  for  geologic  materials.  Modulus  corrections  are 
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consequently  estimated  over  the  range  of  stress  levels  where  these  nonlinear  phenomena  occur.  From 
these  modulus  corrections,  stress- strain  relations  are  obtained  for  geologic  media  over  the  range  of 
stresses  from  acoustic  to  hydrodynamic. 

By  introducing  the  concept  of  wave  front  shell  thickness  and  by  considering  the  energy  associated 
with  an  explosive  wave  disturbance,  simple  ordinary  differential  equations  aie  derived  whose  solutions 
give  analytic  expressicms  for  peak  stress  as  a  function  of  distance  from  explosimis.  The  shell  thickness 
of  a  spherical  wave  front  r^resents  the  distance  behind  the  fr<xit  which  contains  essoitially  all  the 
energy  carried  along  by  the  wave.  Explicit  consideration  of  the  time  d^endence  of  the  wave  disturb¬ 
ance  is  avoided  by  use  of  the  shell  thickness  concept,  and,  because  of  this,  calculations  are  greatly 
simplified  and  result  in  analytic  solutions. 

Results  of  this  phoiomenological  description  provide  estimates  of  peak  stress  at  all  distances 
from  buried  explosions:  these  estimates,  when  compared  with  experimental  data,  are  as  accurate  as 
those  of  the  more  sophisticated  calculations  obtained  from  electronic  computers. 

The  shortcomings  of  the  description  are  evident  and  have  been  emphasized.  At  the  expense  of  mathe¬ 
matical  rigor  and  detailed  description,  simplifying  assumptlcms  are  introduced  in  order  that  an  analytic 
solution  may  be  obtained.  Nevertheless,  despite  the  approximate  nature  of  the  equations,  values  of 
stress  versus  distance  may  be  computed  for  buried  explosions  which  agree  with  measurements  to  better 
than  an  order  of  magnitude  over  the  complete  range  of  stresses  from  hydrodynamic  to  acoustic. 
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